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HAUPTAUFSÄTZE 


Spannungszustände und projektive Iransformationen. 
Von Robert Sauer in Aachen. 


ie bekannte „Dualität“ zwischen Lageplan und Kräfteplan eines Fachwerks läßt sich auf 
D die Differentialgeometrie übertragen und führt dort zu einem „dualen“ Entsprechen der 
Sätze über infinitesimale Verbiegungen von Flächenhäuten und über Spannungsverteilungen 
in Membranen. Diese Beziehungen wurden erstmals von W. Blaschke') deutlich erkannt 
und später von M. Lagally’) weiter untersucht und geklärt. Vor kurzem habe ich nun 
projektive Fragen der infinitesimalen Flächenverbiegung untersucht?) und eine einfache 
Methode gefunden, um aus den infinitesimalen Verbiegungen einer Fläche die infinitesimalen 
Verbiegungen der zu ihr projektiven Flächen’) zu berechnen. Die duale Übertragung dieser 
Methode liefert einen auch für die praktische Durchreehnung sehr einfachen Weg, um 
aus denSpannungsverteilungeneinerMembran dieSpannungsverteilungen 
der zu ihr projektiven Membranen zu gewinnen. Im I. Abschnitt «der vorliegenden 
Arbeit werde ich diese Untersuchungen für unebene und ebene Membranen durehführen und 
dureh Hinzunahme äußerer Kräfte noch verallgemeinern; der Il. Abschnitt behandelt das 
entsprechende Problem für räumliche Spannungsfelder. 

Die Arbeit ist so geschrieben, daß sie auch für flächentheoretisch wenig interessierte 
l.eser und ohne Kenntnis der zitierten Abhandlungen verständlich ist. An zweı einfachen 
Beispielen wird das Problem verdeutlicht. 


Il. Unebene und ebene Membranen. 


1. Bezeichnungen. Auf eine vorgegebene Membran wirken längs des Randes und in den 
einzelnen Flächenelementen vorgegebene Kräfte (Rand- bzw. Flächenkräfte): diese sollen 
sich im Gleichgewicht befinden mit den von ıhnen ın der Membran hervorgerufenen 
Tangentialspannungen. Ich setze zunächst die Membran als uneben voraus und werde 
die ebenen Membranen als Ausartung der unebenen in Nr. 6 behandeln. 


1), W. Blaschke: Reziproke Kräftepläne zu den Spannungen in einer biegsamen Haut, Int. Congr. of Math. 
Cambridge 1912, Proceedings 2 (1913), S. 291-214. 

°), M. Lagally: Über Spannung und elastische Deformation nnebener Membranen, ZAMM 4 (1924), S. 377 bis 38, 

3), R. Sauer: Wackelige Kurvennetze bei einer infinitesimalen Flächenverbiegung, Math. Ann. 1993, S. 673 bis 60; 
R. Sauer: Krümmungsfeste Kurven bei einer infinitesimalen Flächenverbiegung, Math. Zeitschrift 38 (10234), S. 468 bis 475. 

4) Vel. hierzu auch G. Darboux: Theorie des surfaces IV (18096). S. 78, 
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Wir denken uns die Membran als Fläche ®=1 in ein Gewebe dreier Flächenscharen 
ıt == const, # == const, w — const eingebettet; die Flächen «== const, » = const schneiden die 
Membran w=1 nach zwei Kurvenscharen. Das Flächengewebe soll so gewählt sein, daß in 
der Umgebung der Membran «, ®, w als krummlinige Koordinaten der Raumpunkte benutzt 
werden können. Dann ist der Ortsvektor eines Raumpunktes eine eindeutige Vektorfunktion 
vr (u, », m); insbesondere liefert x(w, r, 1) die Punkte der Membran. 

Die 6 Flächen wu, u + du, v, r+dev, w, w-+ dw schließen ein elementares „Parallel- 
lach“ ein, das im Falle »» = I mit der Seitenfläche (x,xXx.) du de der Membran anliegt (Abb. 1): 
die beiden Kanten x„dw, x„dwv liegen in der Membran, die dritte vom Punkt u/v/w aus- 
gehende Kante x. dw tritt aus der Membran heraus. Jedem Membranpunkt ordnen wir nun 
das durch X. X Xu bestimmte 1. a. schiefwinkelige Koordinatensystem zu, wobei die Einheits- 
punkte auf den Koordinatenachsen die Endpunkte der Vektoren x, X, Lu Sein sollen. Auf 
dieses von Punkt zu Punkt verschiedene Koordinatensystem beziehen wir die Flächenkräfte 
und die Spannungen durch den Ansatz 


(PXut 9X» + ri.) dudv — Flächenkraft auf das Flächenelement (r,„Xx,)dud ”) 


(!!x, +0" r,)du | | I; [ru du . i (N), 
Spannung am Linienelement 
| | Lu dr 


(0°! x,+ 0°? x.) dv 


welche von der umgebenden Membran auf das Flächenelement (x„Xr.)dud®e wirkt. Die 
Randkräfte können im Gleichgewiechtszustand als Spannungen angesetzt werden. Die Zeiger 
der oi® sind oben angebracht, um später nicht durch die Differentiationszeiger u, v» gestört 
zu werden. 

Die Spannungen o'’”’r„du und 0°"x„dve sind Schubspannungen; o'r,du und 
6” rY,„dr werde ich als Querspannungen bezeichnen. Wenn die Parameterkurven ein 
Orthogonalsystem bilden, greifen die Querspannungen senkrecht am Linienelement an. 











/U#i tl 
f 7 2° 14 Lu L ? 
u/v O,, Au urdu/v 
Abh. 1. Abb. 2. 


Aus den Spannungen an den Linienelementen x,du, x„dv erhält man die an einem 
beliebigen Linienelement x,du  xr.de® (vom Punkt «/»-+-dv zum Punkt «—+du/v; vel. 
Abb. 2) wirkende Spannung durch Vektoraddition (0'' x, + 0" x,) du +(0° x, +0” x,) dv. Die 
Spannungsverteilung in der Membran ıst also dureh die o°* vollständig festgelegt. 

Für die praktische Durchführung ist es nötig, die Flächenkräfte und die Spannungen 
Lu In 


. 


I In 


auf die Flächen- bzw. Längeneinheit zu reduzieren und auf die Einheitsvektoren 


Um ‚ en , 
zu beziehen. Man erhält als reduzierte, auf Einheitsvektoren bezogene Komponenten 
TE 


der Flächenkräfte bzw. Spannungen (alle Quadratwurzeln positiv!) 


( WER 4G. 
Ph ra =; Q=g \ In E; I r V Is: 5; 
| Yı1 Is» ir | Jıı Is» Is | YJıı Is» Ss (2) 
zu or NEF) & 22 o'? j N Pr. . BEZ 0” Ir 
ii Js» 
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2. Gleichgewichtsbedingungen. Das elementare „Parallelogramm* (1,%x x.) du dv befindet 
sich im Gleiehgewieht unter dem Einfluß der Flächenkraft (px. + 9X» + rx„) du dv und der 
an den Parallelogrammseiten wirkenden Spannungen. Der Flächensatz liefert 


ot oki... EITE ; ; j an 2.5 (>), 


und aus dem Schwerpunktsatz folgt dann die Vektorgleichung 


15) | 


. . 2 | „22 u Il » .I2 0 | 12 - 22 . . er 
a Pe 7: 1 Kun +0’ Luu + 00 £v Foy"LutIu Lot 0Iu LutrpPp lu tg tr mi, 1 


Die skalare Multiplikation’) von (4) mit Xu X X» bZW. Lu X Kuu bzw. Lu» X Xu» ergibt die drei 
Bedingungen 


akt 7 Yv, Lv») +20" (t4,Xv:Luv) + 0° (Ko, Ku u) tr Ku Ko Low) v, 


"kun kr t20"7 an Klum kun (rn +0u ) En ku tn ku dt rn Ku ud =O, I (D). 
0°’ lkKky, vs du u? 20" (2.,X du BD, (04°” t ou IKKvı ku Ko nv) PK». Xr NER IT) Zu m vo Kran Ku) oJ 


Für orthogonale Parameterkurven spezialisieren sich die Gl. (5) zu den bekannten Gleich- 
sewiehtsbedingungen von L. Lecornu und E. Beltramı?’)®). Durch Integration der drei 
partiellen Differentialgleichungen (5) für die drei gesuchten Funktionen 4 „, erhält man die 
mit den vorgegebenen Kräften im Gleichgewicht befindlichen Spannungsverteilungen der 


Membran. 


3. Projektive Transformationen. An Stelle der kartesischen Koordinaten .r, 9, z führen 
wir homogene Koordinaten ein 


& 7 Ü ; 
K Be = 3 z ; : ’ : (bb). 
ı) i ı) ı) 
i r ’ r i a 
Dadurch transformieren sich die Gl. (5) ın 
ıı I [ 5 R alte £& E c | 2212 © c C I. r |. - & 0 
0 SS seen -0 SS Sur 0 Ss, Sm Sm Sun] Sin Sim Sn . 
‘ ( 
Dr ı) 
11fe £ & & he E B = a 12 9,12 _% | - a 
0 Eu use ol 30 IE FRE m lo, 0 ) Oy 06 ,) S, Se Sm Sun ‚MM. 
I _ - | Sn . — 
1 ql: FO FRE PR T r|s Eukuntel 0 





Die nieht hingeschriebene dritte Gleichung folgt aus der zweiten durch Vertauschung von 
I und 2, « und », p und q. Die Ausdrücke |. . . .| bedeuten vierreihige Determinanten, in 
deren 2., 3. und 4. Zeile die zur angegebenen 1. Zeile analogen Glieder mit »,. Z und 9 stehen. 


Wir betrachten nun eine projektive Transformation: £, 7, £, 9 = Linearformen von £, 1. 
©, 9 mit konstanten Koeffizienten, deren Determinante nicht verschwindet. Ersetzt man £, ». 


A Rz ’ | Me: B\ 2 'd\? j B\ _ "\: 
©, 9 dureh & 9, 5, 9 und of, p, gr durch o’F- | ) ok,» 17 pP» q —| 9 gr | 9 r, 
so ändern sich die Determinanten |... .| um einen gemeinsamen nicht verschwindenden 
konstanten Faktor und die in (7) auftretenden Koeffizienten der Determinanten um den ge- 


j i J\° 
meinsamen Faktor | ,) i 
! 


Dies führt zu folgendem Ergebnis: 


Aus einem Gleichgewichtszustand zwischen den Flächenkräften p, q r. 
den Randkräften und den Tangentialspannungen o’*" einer vorgegebenen 
Membran (M) läßt sich ein Gleichgewiehtszustand für jede zu (M) projek- 


tive Membran (M) ableiten, indem man für (M) ansetzt’) 


s=wp, d=og4, Far, o'k—mok mit (S). 


5) Das skalare Produkt (x X v)3 wird mit Ga, d, 3» bezeichnet. 

6) J.. Leecornu: Sur V’equilibre des surfaces flexibles et inextensibles, Journ. de P’Eeole polyteehnique 29 (48), 1880; 
KE. Beltrami: Sull’equilibrio delle superfiecie flessibili ed inestendibili, Mein. della Ace. di Bologna 4 (0), 1881: vel. 
ferner B. Caldonazzo: Sulla meeeaniea delle superfieie, II Monitore teenico 26 (1920), S. 5 bis 20; hier wird insbesondere 
gezeigt, daß die Gleichgewiehtsbedingungen nieht nur für undehnbare Membranen gelten. 

7, Bei der Transformation der reduzierten Kräfte und Spannungen (vgl. (?2)) kommen neben © noch Längen 


und Flächenverzerrungen als Faktoren hinzu. 
13” 
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Li 


Geometrisch bedeutet den mit einer für alle Punkte gemeinsamen Zahl 


) 
multiplizierten Abstand der Punkte der Membran (M) von der Fluchtebene, 
die bei der projektiven Transformation in die unendlich ferne Ebene ab- 
gebildet wırd. 

Spezialisiert sich die projektive Transformation zu einer affinen, so wird » konstant, 
d.h. für jede zu (M) affine Membran (M) existiert ein Gleichgewichtszustand mit ungeänderten 


p. q. r, o’": die reduzierten Kräfte und Spannungen ändern sich nach (2). 


Unter der besonderen Annahme, daß keine äußeren Kräfte auftreten (p= gr), 
ist der Inhalt dieser Nummer das mechanische Gegenstück zum flächentheoretischen Problem. 
die infinitesimalen Verbiegungen der zu einer gegebenen Fläche projektiven Flächen zu 
bestimmen®)®). 


4. Beispiel. Die vorgegebene Membran (7) sei die Kugelfläche #°?+ ?—+z°+1: es 
bestehe Gleichgewicht zwischen einem konstanten, von innen her wirkenden Normaldruck 
1. > . . a . > n ” = A 

(Flächenkraft) und einer isotropen Spannungsverteilung in der Kugelfläche. Durch & 

Y | i 
„,» . 2°, wird die Kugel (JM) transformiert in das zweischalige Drehhyperboloid 
(MA) y’--2°-+-1=0. Nach Einführung der räumlichen Polarkoordinaten . = m cos u cos v, 
y mweosusinr, 2 wesin a hefert die Übertragung des vorgegebenen Gleichgewichtszustandes 


der Membran (M) auf die projektive Membran (M) nach Nr. 3 mit © sin?’« folgendes 
Ergebnis: 


Membran (M) 


P=-0=-0, R=2a: a Zt 212 — 0: 


p fi 0, r 2acosu: Ö R OT —ACOSUM, oe" 0: 


Membran (M) 


n n > . n jo sın? u - 3% 4 1 = ” 
»—=g=V, "= »>asinucosu: o'!—a ; 6” asıin’ an cost, ee | E 
COS MH 

‚2a nd 
i i : sin’ u = sın? u ER 5 12 
P=0-0, R=2ua —; Zrmg 2 =onauyil oe, 2-0, 
. ] l + COS” MH ] | —+ cos’ u 


Da die Durchmesser der Kugel (M) in die Parallelen zur Z-Achse transformiert werden, 
sind alle auf die Membran (M) wirkenden Flächenkräfte zueinander parallel und können als 
Schwerkräfte aufgefaßt werden, wenn man sıch (M) in passender Weise inhomogen mit Masse 
belegt denkt. Schneidet man von der Kugel (M) bzw. dem Drehhyperboloid (M) eine Haube 
ab, so besteht Gleichgewicht, wenn man die längs des Randes der Haube wirkenden Spannungs- 
kräfte als Randkräfte hinzufügt. 


5. Spezielle Parameterkurven. Ein «, »-Kurvennetz, bei dem nur Querspannungen (0'’ 0) 
bzw. nur Schubspannungen (6'' = 0° 0) auftreten*®), soll kurz als Querspannungs- bzw. 
Schubspannungsnetz bezeichnet werden. Indem man vermöge u" u" (u,v), «vu, r) 


neue Parameterkurven einführt, für diese die 6*’" berechnet und 0*"=-0 bzw. 0"! 0°" 0 
fordert, findet man folgendes Ergebnis: 


Zu jeder vorgegebenen Spannungsverteilung gibt es ein Schubspannungsnetz und unend- 
lich viele (uerspannungsnetze. Das Schubspannungsnetz ist bestimmt durch die Differential- 
eleichung 0" du’ +-20"dude +0”de”’- 0 und besteht aus zwei zueinander konjugierten 
Kurvenscharen; diese sind reell und voneinander verschieden für 0'' 06°” — (o'’)’<0®. In den 
(uerspannungsnetzen kann die eine Kurvenschar (unter Vermeidung der Kurven des Schub- 
spannungsnetzes) willkürlich gewählt werden, die andere Kurvenschar ist dann dadureh 
bestimmt, daß die Kurven des Schubspannungsnetzes und ‚der Querspannungsnetze sich 
harmonisch trennen, 


s) Vgl. hierzu die in 8) zitierte Arbeit von E. Beltrami. 
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Die Bedeutung der Schubspannungs- und (uerspannungsnetze liegt darin, daß sie bei 
der projektiven Übertragung nach Nr. 3 ihren Charakter beibehalten. Das flächentheoretische 
Gerenstück sind die krümmungsfesten®) und die wackeligen®?) Kurvennetze. 


Die Querspannungsnetze lassen sich durch Fadenmodelle mit zwei Scharen gespannter 
gegenseitig verknoteter Fäden darstellen. Die in den Fäden auftretenden Spannungen geben 
die Querspannungen und werden durch die in den Verknotungen angreifenden äußeren Kräfte 
im Gleichgewicht gehalten. Die von den Fäden überdeckte Fläche braucht im Falle ver- 
schwindender äußerer Kräfte nicht durch ein starres Material realisiert zu werden. 


6. Ebene Membranen. Die Membran (M) sei jetzt als eben vorausgesetzt (2= 0). Die 
Vektorgleichung (4) liefert dann nur zwei partielle Differentialgleichungen für die drei Spannungs- 
[unktionen ei" (u, r). Dies führt zu wesentlichen Unterschieden der Probleme der ebenen und 
der unebenen Membran; z. B. gibt es zu Jedem Kurvennetz einer ebenen Membran Spannungs- 
verteilungen, bei denen das Kurvennetz (uerspannungsnetz ist, während bei unebenen Mem- 
branen nur spezielle Netze als Querspannungsnetze genommen werden können. Die in Nr. 3 
behandelte projektive Übertragung läßt sich jedoch auch für ebene Membranen (MM) vornehmen: 
das in Nr. 3 abgeleitete Ergebnis bleibt dabei wörtlich gültig. 


II. Räumliche Spannungsfelder. 


Die Untersuchungen des I. Abschnitts sollen jetzt auf räumliche Spannungsfelder über- 
tragen werden. Die analog verlaufenden Rechnungen werde ich dabei nur kurz andeuten. 


7. Bezeichnungen. In einem vorgegebenen räumlichen Feld bestehe Gleichgewicht zwischen 
den an der Begrenzung und in den einzelnen Raumelementen angreifenden Kräften (Be- 
serenzungs- bzw. Volumenkräften) und den von ihnen hervorgerufenen Spannungen. 
Wir betrachten wieder das in Nr. I eingeführte Parallelflach vom Rauminhalt a... 1.2 du de dw, 


das jetzt aber nicht an die Fläche w= 1 angrenzen muß. Ferner setzen wir 

(Pu tg» + rr.) du dedw = Volumenkraft im Raumelement Q,. 1,1.) du de dw, 
(ut. ter.) dv dw | (tr Xi) dv dw 

o +» DH x ar ' x (9), 
(0°"' x, +0” x, +0” x.) dw du Spannung am Flächenelement 4 X. ı.) dw du 
(o®' » | 0°” E 1 0°? Bu d u d 1 | (Ku > vs du dv 





welche von der Umgebung auf das Raumelement %,, X, X.) du dv dw wirkt. Die Begrenzungs- 
kräfte können im Gleichgewichtsfall als Spannungskräfte angesetzt werden. Analog zu Nr. I 
ist durch die 0" die Spannungsverteilung vollständig festgelegt. 


Durch Reduktion auf die Raum- bzw. Flächeneinheit und durch Beziehung auf die Ein- 
‘ T In Lu . b . n u 
heitsvektoren ©", 0.” erhält man für die Volumenkräfte bzw. Spannungen (alle 
Hu Er „u 
(uadratwurzeln positiv!) 





( Js E 
P=p jl () g) MER R .. MER 
4 6 Y 
19 ) 9). 
li — gli Iii Sy; ai Iwii win wii 
Vor Den — In I Fr U 


wobei die g;, entsprechend Nr. 1 definiert sind und außerdem yg den positiven Wert von 
Ku Un, Lo) bedeutet. 


8. Gleichgewichtsbedingungen. Flächensatz und Schwerpunktsatz liefern 


EEE rin a ren Ye 
il » | 22 » 33 8 23 u ı O9 nsi .) „I? x 
oo Kuur O0" Ivo + 0 Kr u 20 Ivw 9 Ivwu r=0 Iuv | 
| 11 » 22 33 » - 28 m 23 „ - Mn | = 12 » | = 12 „ 
rOu Yu + Gy In + Ow I u + On Im tr Ow I. + Om Lu 64 kIwtrOu Io tr 0% Lu . . 1 ). 
+ pPXu + ET + "Lır 0 | 


Die skalare Multiplikation von (11) mit x. X x. bzw. X. X X. bzw. 1. x X, ergibt drei partielle 
Differentialgleichungen für die sechs Funktionen or.) - 
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9. Projektive Transformationen. Durch dieselben Überlegungen und Rechnungen wie in 
Nr. 3 findet man das analoge Ergebnis: 


Aus einem Gleichgewicehtszustand zwischen den Volumenkräften p,g, r, 
den Begrenzungskräften und den Spannungen o’* eines vorgegebenen räum- 
lichen Feldes (F) läßt sich ein Gleichgewiehtszustand für jedes zu (F) pro- 
jektive Feld (F) ableiten, indem man für (F) ansetzt’) 


ır 
_ 
F 
« 


’ I 
) vd», ( (vd, Y DT, oth aork mit | / 5 ö . . (12). 
Pop Img mo 
ı) 


Uber die geometrische Bedeutung von 9 
! 


vel. Nr. 3. 


10. Beispiel. Das vorgegebene Spannungsfeld (F) sei die Vollkugel #’+ 1? +2? 1; es 
bestehe Gleichgewicht zwischen einem von außen her wirkenden konstanten Normaldruck 
auf die Begrenzung (Begrenzungskraft) und einer isotropen Spannungsverteilung im Innern 
der Vollkugel. Volumenkräfte sollen nicht auftreten. Durch dieselbe projektive Trans- 
formation wie in Nr, 4 wird die Vollkugel (F) in zwei vom zweischaligen Drehhyperboloid 
“94° 2° 410 begrenzte, ins Unendliche sich erstreckende Raumgebiet abgebildet. Nach 
Nr. 9 erhält man mit © = w” sin’ u folgendes Ergebnis: 





Spannungsfeld (FF): 
ne EEE ER ni 
7 | N 7 
g*’ LCOSH, "rad Fa 0” am“ cos u; Volumen- 
che kräfte und 
" Schub- 
Spannungsfeld (F): 
RR spannungen 
> ! no „sin? u Kr 2 ae 0 
0 aa” sın? un cos, 6" —=amw? ’ 6" = a w*sın? u cos u, 
COS MH 
SIT 8 cin? ; ; ds ins Sega m’ sin?! u 
P am sınT u yw” + cos’ u, a mw sın U, ” Aue Ad z = 
y mw’ —+-cos” u 


Die an (FF) angreifenden Beerenzungskräfte müssen gleich den dort aus unserer Rechnung 
folgenden Spannungskräften I°® angenommen werden. Da sie parallel zur Z- Achse wirken, 
kann man sie als Schwerkräfte auffassen, wenn man sich die Begrenzungsfläche in passender 
Weise ınhomogen mit Masse belegt denkt. 


11. Spezielle Parameterflächen. Kin ı-, v-, w-Flächengewebe mit 0° = 0°' 00 bzw. 
oo" eo" 0 soll Querspannungs- bzw. Schubspannungsgewebe heißen, da in 
diesen Fällen nur Quer- bzw. Schubspannungen auftreten. Da für die 6 Funktionen o"(u,r,w) 
nach Nr. S nur drei Bedingungen vorgeschrieben sind, gilt hier folgender von Nr. 5 ver- 
schiedener Sachverhalt: 

/u einer vorgegebenen Spannungsverteilung gibt es nicht nur unendlich viele (Quer- 
sondern auch unendlich viele Schubspannungsgewebe. Für jedes Flächengewebe kann, man 
Spannungsverteilungen angeben, bei denen es (Quer- bzw. Schubspannungsgewebe ist. 

Natürlich sind die Quer- und Schubspannungsgewebe im nämlichen Sinne projektiv- 


- 


invarlant wie die Quer- und Schubspannungsnetze der Nr. 5. 


Die hier für Gleichgewichtszustände behandelte Methode läßt sich auf beliebige Be- 
wegungszustände verallgemeinern. Dabei kommen zu den Spannungen und Flächen- 
bzw. Volumenkräften noch Beschleunigungen hinzu, welche bei der projektiven Übertragung 
ebenso wie die Flächen- bzw. Volumenkräfte transformiert werden. Ferner gelten, wie ich beı 
späterer Gelegenheit genauer auseinandersetzen werde, analoge Sätze auch für Fachwerke. 
Man gewinnt dadurch aus den Spannungsverteilungen eines Fachwerks die Spannungsver- 
teilungen der zu ihm projektiven Fachwerke. Unter der besonderen Annahme verschwindender 
äußerer Kräfte ergibt sich der Satz von H. Liebmann'") über die projektive Invarianz der 
„Ausnahmefachwerke*. 421 

») Bei der Transformation der reduzierten Kräfte und Spannungen (vgl. (2)) kommen neben ® noch Längen-, 


Flächen- und Volumenverzerrungen als Faktoren hinzu. 


10), MH. Liebmann: Ausnahmefachwerke und ihre Determinante, Münchener Berichte, math.-phys. Klasse, 50 (1920), 
S. 197 bis 227 
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Ballistik '. 


Von ©. ». Eberhard in kEssen-Bredeney. 


m Handbuch der physikalischen und technischen Mechanik von Auerbach und Hort 

1928 Band II, Lieferung I, Seite 240 u. folg. konnte Verfasser zeigen, daß unter Zugrunde- 
legung eines Exponentialgesetzes für die Abnahme der Luftdiehte mit der Höhe über dem 
Erdboden alle senkrechten Flugbahnen der Geschosse mit beliebigen Anfangsbedingungen und 
Parametern sich in ein einfach unendliches System von im voraus berechneten Bahnen ein- 
ordnen lassen; des ferneren, wie man von vertikalen zu geneigten Bahnstücken übergehen 
kann, wofern diese Bahnstücke nur genügend wenig gekrümmt sind. Die Lösung setzte die 
klassische Formel für den Luftwiderstand 


dv Pi, ” 
Ir, 


nt fir) : Ö Io " K (vw) 
At I g I Yy 


voraus, wo » die Geschoßgeschwindigkeit, f die Zeit, r den Geschoßhalbmesser, P das Ge- 
schoßgewicht, ö den ballistischen Formwert, ö=og das Gewicht eines ebm Luft darstellt 
und g, zum Beispiel bei Fasella gleich 9,5047 m/s? gesetzt ıst. (Im Nachfolgenden sind noch 
an Bezeichnungen verwendet: s die Schallgeschwindigkeit bei der absoluten Temperatur T 
;, Schallgeschwindigkeit bei der absoluten Temperatur T,, für welche A «(r) ermittelt ist, 
I?, die Gaskonstante, p der Luftdruck, d der Neigungswinkel der Bahn mit der Horizontalen.) 


, 
N 


Inzwischen ist erkannt worden, daß die Verzögerung durch den Luftwiderstand richtiger 
eeschrieben wird 


dev rr ..., .(® S, 
dt P Ö Tin Ki, 
oder 
dv ri .9, 8 ( 8. „[V8, ri. 
7 re | KI—®): Ar, v”), 
s? T j %) ae s? pP. 
und da ö -T Ö TI also mit p = z I, Ö RT, ist, wird auch 


Fr PR, „flo : | n2 yole : “Pr ol: = 


Es soll nun gezeigt werden, daß auch für das neue Luftwiderstandsgesetz sich eine 
lösung für den vertikalen und für den gestreckten geneigten Schuß finden läßt, bei welcher 
man jede Bahn, also jede Anfangsbedingung und jeden Parameter in eine einfach unendliche 
Schar von Lösungen einordnen kann. 


Die Energiebilanz gibt für eine beliebige Bahn, wenn y nach abwärts positiv gerechnet wird, 


= ydy 2 y oder 2gy - 


2 sın d) dy sın d) 


Wir transformieren die Variablen in der von Langevin (Memorial de Yartillerie 
[rancaise I 1920, Seite 260 u. folg.) angegebenen Weise, indem wir setzen 


15} 65) 


s” Be: __ Ss s 
dy= ..dı; dr ME; dt=—drt; " ERS ER VE 
. Ss, ; a 5 5 
0 0 0 0 
dis? ww?) 2c,pP (mw) 
Dann wird obige Gleichung 5 IL a 
“ oO s’dy J sın ı) 
88 
D} 2} „ 2} Du} + „ . ( 
, \ d(sm) dw  m’ds’ dw? m” ds? dw? „Fr 
Nun ist z + — + — —— + mw? z 
s’dın dn s’ dy dn s. dy dn dy 


I) Ergänzungen zu dem Aufsatz: Benutzung der Fasella Tabellen usw. (Diese Zeitschrift Bd. XI (1031) S. 253 IT.) 





Ztsehr. f. angew. 
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N Ss” ' le) IAdT 
A PX N Ist, wıre dy Tdy 


Somit lautet unsere Gleiehune in den neuen Koordinaten 


dw’ w’dI__ 2e,pP im) 


dı 7 5 J sın #) 


f 0 


Ferner ıst bekanntlich 


’ 1 Ä N) RT dp gay ys 
‘ / une ’ I, somit — —(ö} 
”r u P 7 ” p IR, / I, m / 
] / ‘4 
S . «pp 4 KH Ta 
ume mit ’ - wird „dn und y ‚.e 
( N) S, Br p R, F; / / / 
Die Gleichung wırd also 
47 
lan“ n® AT ) pie X Dim) 
an ; r - 
dn I,dy J sın ı) 
| ’ AT ” BE 
Bisher ist noch keine Annahme über eemacht. Wır setzen nun ). nehmen also 
dy dy 
einen konstanten Temperaturgradienten an. Hiermit wird 
7 
dm? / Sc, Pp eı To pw) 
w—.g _— (2). 
dı) 4 sın ı) 


Für den vertikalen Schuß wird sin 9 1, für geneigten Schuß wird sin d = (sın V),. 
ein Mittelwert. 


Bei Anderune der Werte von e 
71 ‘] (] 


bzw. e In € ı do, kA 


„ und 7, ändert sich, wenn g konstant bleibt, nur c,, 
e also insgesamt nur der konstante Koeffizient des letzten 
(iliedes rechts. Mithin lassen sieh sämtliche Lösungen, genau so wie im Handbuch von 
Auerbach und Hort angegeben, finden, wenn man die obige Gleichung für alle mög- 
lichen e, Integriert. 


8” rn rn n . m . rn Y 
Aus dy = „dan N dı;, und aus 7 l,+49y, wobei T, gleich T an der Stelle y=-0, 
h I,dy 
findet man dıy), also 
I,„+ry 
T,,. Is tiy r 
7 IM -Iın„ — (3) 
/. „ti 
. . . » S . . 
somit , als Funktion von y und umgekehrt. Ebenso findet man » mW, und es ergibt sıch 
dy 
aus di "dureh Quadratur F. 
” 
Für den schrägen Schuß wird bekanntlich 
( (t] 
ydy | 4 -; / ydn 
tv VJdV > : oder hier te d) Ad) — — 
B Ei — 70° 
(Se 
Sn 
| | ydy 
also low” cos ı, loe'" eos ı) 0,13429 \ - (4). 
on 
0 2 


Um 9 zu finden, braucht man also auch nur nach Lösung der Hauptgleichung (2) obige 
Integrale für alle Parameter e, ein für allemal auszuwerten. 


FE 


. 





Lew, 
lech. 


‚lso 


h. 
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Das Verfahren bildet eine Ergänzung zu der stufenweisen Flugbahnberechnung des Ver- 
fassers (diese Zeitschrift 1931, S. 253 u. folg.). Statt des /-Wertes tritt hier (sin 9), auf. 
Dieser Mittelwert kann in gleicher Weise wie dort gewonnen werden, und der Unterschied 
zwischen dem Anfangs- und dem Endwert sowie dem Wert von sin (#,,) für die Mitte des 
betrachteten Bahnstückes bildet auch hier ein Kriterium für die zulässige Länge des be- 
trachteten Bahnstückes; dies kann man auch so ausdrücken, daß =; klein bleiben muß. 

in 

Während aber bei der stufenweisen Berechnung die rasche Änderung der dort benutzten 
»,(®) Kurve in einigen Geschwindigkeitsbereichen dazu zwang, die Stufen kürzer zu be- 
messen, fällt diese Beschränkung hier weg, wenn man nur die Gl. (2) für die richtige 
Div) Kurve integriert. Ferner entfällt die Beschränkung auf nicht zu große Höhenstufen, 
die dadurch bedingt wird, daß bei der stufenweisen Berechnung ein mittleres T und 
ein mittleres Luftgewicht eingesetzt werden muß. Man wird also die Methode bei Steil- 
bahnen in den steilen Bahnstücken, ferner bei sonstigen Bahnen in dem Teil, in dem die 
 (v)Kurve Maxima und Minima aufweist, mit Vorteil benutzen können, desgleichen für Über- 
schlagsrechnungen zur Beschleunigung der stückweisen Berechnung vom Jahre 1931. 

Noch eine Bemerkung ist nötig. Wir haben g als Konstante behandelt, während bei 
der Stufenbereehnung 9, insofern es das Geschoß beschleunigt oder verzögert (nicht aber, 
insofern es die Änderung der Luftdiehte mit der Höhe beeinflußt), als Variable behandelt 
wurde. Um also für ein anderes, mittleres 4 das richtige Resultat zu finden, wäre es deshalb 
am einfachsten, für ein etwa um 6°/,, kleineres 9, welches im ersten Glied der rechten Seite 
der Gl. (2) eingesetzt wird, ebenfalls die Integration durchzuführen und dann für das mittlere g, 
welches bei dem betreffenden Bahnstück gelten muß, zu interpolieren. 


Anmerkung: Die Langevinsche Transformation der allgemeinen Flugbahngleichungen 
für gekrümmte Bahnen führt an sich auf der Integration nicht leicht zugängliche Differential- 
oleichungen. Deshalb hat sie meines Wissens bisher keine praktische Anwendung gefunden. 
“ 2 h 2 aT . u: 
[rotzdem gibt sie ein bequemes Mittel, für die Annahme Fr /  konst. die Flugbahnen in 

. «N 
ein doppelt unendliches System von Bahnen einzuordnen. 


Die Transformation führt, wie oben dargelegt, zu der Differentialgleichung 


9 
dw mw”, \ 202 et To pw) 

u ä s ar ‚ 
dı, Er. 9 sın 


„ nur die Temperatur, bei welcher die Funktion 
P(w) als Funktion P (») empirisch ermittelt wurde, bedeutet. 


wobei, wie nochmals bemerkt werden möge, T 


Angenommen, die Gleichung wäre für beliebige ce, streng gelöst, so ergäbe sie die 
Beziehungen zwischen w, »; (und 9). 


Hätte man diese, so könnte man durch Quadraturen mittels 


ni) 


s” s Mi 
däy=-—;ädn; " EEE 
5 0 
für beliebige Temperaturen (bei konstantem Normaltemperaturgradienten) » als Funktion von 
y finden, d.h. für eine solche Beziehung zwischen T und y, daß einem bestimmten y ein 
- h dT . ' 
bestimmtes 7 entspricht, während 5 /) bleibt. 
( 

Nun ist es aber keineswegs nötig, vorgenannte Gleichung selbst zu integrieren. Hat 
man vielmehr die in y, 9 und » aufgestellten ballistischen Hauptgleichungen zum Beispiel 
in der früher?) angegebenen Weise integriert, und transformiert man die gefundenen Größen 9, ® 
(und ebenso x und #) mit Hilfe der Formeln 

So 5, . So Su 
dy: „day; w=—.v; dz „da: dr dit 
Ss r S Ss” S 
auf die Langevinschen Variablen, so kann man dann aus diesen die Werte von .r, 9, v und f 
für eine beliebige andere Temperatur im Gipfel bei normalen Temperaturgradienten mit Hilfe 
m u .„_dE . u 
der Transformationsformeln (5) ermitteln, solange man nur 4 7 = nicht variiert. 


?)O.v. Eberhard: Ds. Zs. Bd. XI (1931) S. 253. 
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Man braucht also nur für eine endliche Zahl von Gipfelgeschwindigkeiten und e-Werten 
Fluebahnen mittels des gewöhnlichen Verfahrens’) zu errechnen, auf &,n7,w und r zu trans- 
formieren und hat damit eine Flugbahnschar, in welcher man durch Interpolation jede beliebige 
Bahn in &,97, mw und r finden und diese auf ., 9,» und # transformieren kann. Dabei wird 
das bekannte e, im Abschußpunkt, an dem der Druck p, herrschen soll, 

‚Po, To ® 
0 ı) 


Pu 


(7 
r Ü 


wenn e, den Wert von e im Gipfelpunkt bedeutet, also kann man setzen 


f 
0 
j’ (+ ) R 
ey 


Man möchte zunächst glauben, daß der Übergang von .r, y,t auf &, 7,7 und umgekehrt 
der Übereane von £, y,r auf .r, 9, t Quadraturen erfordere. Das ist aber nicht der Fall, wie 
nachfolgende Betrachtung zeigt. 

Bezeichnet man mit T, die Temperatur im Flugbahngipfel, mit 7, diejenige am Ab- 
schußpunkt, und rechnet man y und 7 vom Gipfel anfangend positiv nach unten, so ergeben 
die Transformationsformeln durch Integration 


| Nnı 
„ l p , 
An |e l (n), 
E / RU 
also auch 
( | Y 
/ - / mn 
/ 1 | T, 


Für .r findet man zunächst 


Mit der vorhergehenden Gleichung wird hieraus 


"1+ Tif (eo) 
x \ J # ? ] de. 


l, 
also wenn £ als Funktion von » bekannt ist, 


„(1+Af,n) ds 


n Al : dı, also ‚, f» (N). 


p 


wie übrigens auch von Langevin gefunden wurde. Weiter ist 


ı/T,tiy 
yerira 
und hieraus folgt, wenn man ebenfalls den Wert von y einsetzt und r als Funktion von 7 
einführt: 


vr, 


Schließlich folgt aus f,(y) und aus T,+/iy=T, 


„ 
T, 


sy 4 
T, = T,1+Af, ml]; somit „= fie). 


J 


Bei den im ersten Teil dieser Arbeit angestellten Überlegungen, welche für den Steil- 


i 1 
schuß und für größere Teilbögen, die nur die Bedingung erfüllen müssen, daß to 9 klein 
bleibt, bestimmt sind, gelten jetzt auch die gleichen Betrachtungen hinsichtlich dessen, daß 
/ . . . 
r und vr von der Temperatur unabhängig sind. 426 
( ‘) 4 


3) OO. v. Eberharda.a.dO. 
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Fin neuer Ansatz zur Lösung räumlicher Probleme der 
Klastizitätstheorie. 


Der Hohlkegel unter Einzellast als Beispiel ). 
Von H. Neuber in München. 


1. Einleitung. Für die Lösung ebener und achssymmetrischer Probleme der Elastızitäts- 
theorie stehen bekanntlich Spannungsfunktionen zur Verfügung, aus denen Verschiebungen und 
Spannungen durch Differentiation von vornherein so gewonnen werden können, daß alle elasti- 
schen Gleichungen befriedigt sind. Dagegen fehlt es für die Lösung allgemeiner räumlicher Pro- 
hleme offenbar noch an einer solehen Methode. Der Grund hierfür ist in der Art der bisherigen 
Ansätze zu suchen, die sämtlich auf Integrale führen. Bereits Maxwell?) brachte Spannungen 
ınd Verschiebungen zu drei Funktionen in Beziehung und stellte durch Einsetzen in die räum- 
lichen Grundgleichungen die zwischen den Funktionen bestehende Differentialgleichung auf, 
ohne jedoch diese mit Hilfe eines weiteren Ansatzes in allgemeiner Form zu lösen. Bei 
späteren Ansätzen wurde einerseits derselbe Weg eingeschlagen, wobei es dann gelang, die 
jeweils auftretende Differentialgleichung in Integralform zu befriedigen. So enthält z. B. eine 
dieser Lösungstypen?) eine Integration nach « (cartesische Koordinaten), eine andere eine 
Integration nach r (Polarkoordinaten). Andererseits war man bestrebt, das Problem auf Rand- 
wertaufgaben der Potentialtheorie?) zurückzuführen, was bereits gelungen ist. Die hierbei 
auftretenden Integralgleichungen sind entweder nach der Fredholm schen Theorie oder nach 
lem Verfahren der sukzessiven Approximationen zu lösen. Obwohl jeder dieser Ansätze eine 
mathematisch einwandfreie Lösung des Problems darstellt, ist doch das Aufsuchen strenger 
lösungen in beliebigen Koordinaten mit verhältnismäßig großem Rechenaufwand verknüpft. 
lüs dürfte daher von Interesse sein, im folgenden einen sehr einfachen Ansatz kennenzulernen, 
mit dessen Hilfe es gelingt, Verschiebungen und Spannungen allein durch Differentiation aus 
einer allgemeinen räumlichen Spannungsfunktion zu gewinnen, die sich aus drei harmonischen 
Funktionen zusammensetzt. Das zugehörige Differentiationsschema wird vollständig sym- 
metrisch und läßt sich deshalb leicht auf beliebige Koordinaten übertragen. 


Als Beispiel wird das Problem des Hohlkegels unter Einzellast behandelt. Es werden 
bei drei verschiedenen Belastungsfällen die genauen Ausdrücke für die Spannungen angegeben. 


2. Die Lösung der elastischen Grundgleichungen mit Hilfe des neuen Ansatzes. In carte- 
sischen Koordinaten .r, y, 2 mit den Normalspannungen o,, 9,, 9-, den Schubspannungen 7, ,, 
Ty-, T.„ und den Verschiebungen £, 7, © lauten die Gleichgewichtsbedingungen °): 


00. , OT, OT. 
ET ED 50 ee a re | 
Our 0 02 


usw. mit zyklischer Vertauschung. Ferner bestehen zwischen Spannungen und Verschiebungen 
bei Zugrundelegung des Hookeschen Gesetzes und Voraussetzung kleiner Formänderungen 
bekanntlich die folgenden Beziehungen’): 


I(r OS r ) 
Or ul r usw... . (2), 
" oz m—?2 
wobei 
VS i ol? n OL r (3) 
Or 094 02 
gesetzt ıst, und 
‚{dS ,0n 
Txy (i N + USW. . Fe: a (4). 
er 0y 0x 


I, Vorliegende Abhandlung ist ein theoretisches Ergebnis der Forschungsarbeit, die vom Verfasser im Meeh. 


Teehn. Labor. der Teehn. Hochseh. München auf Veranlassung von Herrn Prof. Dr. L. Foöppl durehgeführt wurde. 


Hierbei möchte ieh mir gestatten, der Notgemeinsehaft der Deutschen Wissenschaft Für ihre liebenswiürdige Unter 
stützung meinen aufrichtigsten Dank auszusprechen. 

2) J.C. Maxwell: Seientifie Pap. of J. ©. Maxwell, Paris 1927, Bd. ?, S. 108 u. f. 

3) E. Trefftz: Mathematische Elastizitätstheorie, Handbueh d. Ph., Bd. VI, S, »., 

4) [,. Liehtenstein: Über die erste Randwertaufgabe der Elastizitätstheorie, Math. Zeitsehr., Bd. 20, S. 21, 
1924; ferner A. Korn: Über die Lösung des Grundproblems der Elastizitätstheorie, Math. Ann., Bd. 75, S. 497, 1914. 


5) A.u. L. Föppl: Drang und Zwang, I. Bd., 1. Aufl,, S. 16 u. f. 
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Durch Kombination mit den Gl. (1) erhält man die sogen. „elastischen Grundgleichungen“ 


1: m de 0 = 
+ BERN ZI a el 
m —20x 
wobei 
0” 0” 0” 6 
ne = z } - } ) £ )) 
0x oy 072 
bedeutet. 


Ks kommt nun darauf an, mit Hilfe eines neuen Ansatzes die Verschiebungen derart 
aus harmonischen Funktionen abzuleiten, daß die Gl. 3) und (5) erfüllt sind, ohne daß 


wie bei den bisherigen Methoden Integrale auftreten. Dies gelingt mit folgendem 
Ansatz: 
\F' 
( 
Gm -(@ 
E- N " f 
\F 
( _ 
2G HU», PP ar Te I A 
ne 7; j 
\F 
( ‘ 
6 ” Ri - r( P, 





Hierbei sind P,, P,, P, harmonische Funktionen, genügen also der Gleichung 
PR Te 5 
Um über den Zusammenhang dieser Funktionen mit der Spannungsfunktion F Aufschluß zu 


erhalten, setzen wir die Verschiebungen entsprechend Ansatz (7) in die Gl. (3) und (5) ein. 
ls ergibt sich mit Berücksichtigung von Gl. (8) 


/AD od N 
’ oP,, o®P. 
IF+0Cl-—+-—2+—)=26e....:....0 
X 07, 02 
une 
Ö om | 
\ | IF+ „26 e) 0 ee ee 
Or m—2 
Die (il. (10) verlangen 
i m 
IF-+ 52Ge=const, 
m—2 


oder, da die auf der rechten Seite stehende Konstante unwesentlich ist. 


.)\ 


8Be=ll- ZAPF. 2.2... rer 


[IE 


Setzt man dies ın Gl. (9) ein, so ergibt sich 


r & D A) PD, A D, 
Or 07] 02 


2l1-—)AF. ........ (2). 


[IE 


Während beı den bisherigen Methoden die entsprechende Gleichung nur 
ın Inteeralform lösbar war, läßt sıe sich hier in durchaus einfacher Weise 
befriedigen. Setzen wir nämlich 


F= PD, + ır D, + PD, n- z D, \ - 6 ; ; . R ri r (15). 
wober auch P, der Gl. (S) genügt, so wird 


Oo», A) PD, H OP», 


ä ») I 
Felt tie) re 


und 
| IF=V0 a : ; . r r j : k ; ® R j . , : r (15). 


Gl. (12) ıst schließlich vollständig befriedigt, wenn wir noch 


C=4[|1 n EEE. 


m 
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Mech. 





Prart 
[3 
lem 


(5). 


3 zu 
ein. 


(9) 


10). 


11). 


12). 


ISe 


14) 


I»). 


16) 
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»tzen. Diese Konstante trat schon beim achssymmetrischen Spannungszustand ®) in Erscheinung, 
ınd zwar wurde 


"0! 77 Er 


m 


sesetzt. Man erhält endgültig für die Verschiebungen 


OF 


2GE= ER BE ea a a 
0x 


Die Spannungen ergeben sich nach kurzer Zwischenreehnung aus den Gl. (2) und (4), wobei 
von den Gl. (11), (14) und (17) Gebrauch zu machen ist, zu folgenden Ausdrücken: 


0° F 0” F (‘ D, 0 D, O 2, 19 
Pe —— = = ya USW, 5 i : } ; . ; +], 
Toy? 02° 0x 0y 02 ei, 
"F « D, oO 2) " 
N — u USW. : ; ; . ; er F . .% zu), 
x u A) x A) Y A) Y A) r 


Man überzeugt sich leicht, daß tatsächlich die Gleichgewichtsbedingungen erfüllt sind. Ebenso 
sind die sechs Verträglichkeitsbedingungen erfüllt, die sich durch Kombinieren der Gl. (2), (3) 
und (4) noch aufstellen lassen. 


3. Deutung des Ansatzes. Der Ansatz soll zunächst an Hand der speziellen Klastizitäts- 
theorie näher diskutiert werden. 

Die Theorie der Torsion prismatischer Stäbe entspricht den Ausgangs- 
eleichungen 


(ı \ (\ N 
D, r D, j D, — D, (1. 2). D,: 2. PD, er, (20a). 
e x [Z ; a 3 
Es wird dann 
Y . a ( 2 ( \ 
F 0, s-= (i PD, n=Uaz, C% Day. 5 + + (a0. 


Die Theorie der Torsion runder Stäbe veränderlichen Querschnittes 
entspricht den Gleichungen 


(i ( 


D,d, 2 =D, = (a, vr), PD, y-o(a,r) . » . . (206), 
2 Ä . 


Qu 


A 
wobei mit «,r Zylinderkoordinaten zugrunde gelegt sind (r = yy’—+ 2”). Es wird 
BT ni fü _ - Pa .) 
F=0, E=b6, N=2W, [ Be ae 
Aus 1P,—0 folgt für g die bekannte Differentialgleiehung 


E hE 30 
Ne De BWuR; 
ee 5 li ah nn at an 
Ri x: ( Y Ri Y 


Der Fall F=O entspricht mithin der allgemeinen Torsıon. 


Die ebene Elastizitätstheorie wird erhalten mit 


Ri) PD’ r ’ 
—-, ®8,=0, ®%,=0, F=F+aß!’. . . (UN), 
Our s 


PD,=P,+ta PD’, D, 
wobei P, und P,’ zwei neue harmonische Funktionen sind, die nur von «= und y abhängen. 
F’ entspricht jetzt der Airyschen Spannungsfunktion. 


Dieachssymmetrische Elastizitätstheorie geht aus derselben Substitation her- 
vor, wenn wir voraussetzen, daß ®, und ®,’ nur von x und r (Zylinderkoordinaten, s. 0.) 


abhängen. F’ ist dann die achssymmetrische Spannungsfunktion '). 


Auf@Grund dieses Zusammenhanges wollen wir Fals „räumliche Span- 
nungsfunktion“ bezeichnen. 


6) H. Neuber: Beiträge für den achssymmetrischen Spannungszustand, Diss., München 1939, S. 3. 
") A. u. L. Föppl: Drang und Zwang, Il. Bd., ?. Aufl., S. 208, München und Berlin 1998. 
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Für den allgemeinen Spannungszustand ist bemerkenswert, daß immer eine der vier 
harmonischen Funktionen gleich Null gesetzt werden darf, ohne die Vollständigkeit zu be- 
einträchtieen. Dies geht aus der Substitution 


\cDb,' \cdb,’ \$b’ od,’ 
( ( ( / ( 
8, =, +2ad' —-I|2 -—- + +2), 8 =B)- : 
0 0 3% Q r I „ () Y A 2 l 1 + Rh r . 
(A )ır 
op», ap», 8) 
$D, m.’ 4. B__ dD, s F=F'’+2a®B,' 
‚ ‘ u 7, ’ 02 . 





hervor. welehe an der Allgemeinheit der vier Funktionen ® 
zu dem neuen Gleichungssystem 


P,, ®,, P, nichts ändert und 


os 


\F’ 
y u ’ , r) - ( 
I D, L. 7 D/' + PD, \ If, S ) e wi 7 PD’, 
(20h) 
\/ \ F’ 
p) ( 4  % ( 
2 \ 28. Ilrl 
7; ; 02 





führt, welches nur noch die drei harmonischen Funktionen P,’, PD’, P, enthält. Es werden 
also in Wirklichkeit nur drei harmonische Funktionen benötigt. Einer bestimmten har- 
monischen Funktion entsprechen mithin gerade drei voneinander ver- 
schiedene Spannungszustände. Die Mannigfaltigkeit des allgemeinen ela- 
stischen Zustandes ıst gleich der dreifachen Mannigfaltiekeit der harmo- 
nischen Funktionen. 


Für das Aufsuchen von Lösungen ist es zweckmäßig, jeweils ein Koordinatensystem zu- 
erunde zu legen, in welchem die Randflächen des Körpers enthalten sind. Aus diesem Grunde 
wollen wir auf krummlinige Koordinaten übergehen. 


4. Übergang auf krummlinige Koordinaten. Gehören zu dem orthogonalen Koordinaten- 
system u, r, 0, dessen Achsen die Riehtungskosinus eos (#, u) usw, mit den ursprünglichen 
Achsen bilden, die Spannungskomponenten 9,, 7,„ usw. und die Verschiebungskomponenten 
U’, V, W, so gilt, da sich letztere als Vektorkomponenten transformieren, 


"= eos (r,u)-£+cos(y,u)-y+cos(z,wW)-y usw... 2 2 20202020... (21). 


Die Riehtuneskosinus erhält man bekanntlich aus 


| ( c 


2. (22 
h,„ou 


eos (.r, ut) = USW, Zi )s 


wobei die A. ",. h,. der krummlinigen Verzerrung Rechnung tragen und aus den Gleichungen 


„ fjba\ , (öy’, (d2\ er 
hu [353 Ha) Han) WM 4 ae er 


zu bestimmen sind. (In Tensorschreibweise wird hy 91.) 


Setzt man für &, 7, g die Ausdrücke aus Gl. (15) ein und beachtet, daß sich auch die 
ersten Ableitungen einer Funktion als Vektorkomponenten transformieren, so erhält man 
schließlich 

| OH Va \y Oz 


N ( 
l | -P>aw H-2a®, -Da®D 
2Gh,\ du ou om "ou 


usw. . . . . (24). 
Für die Ermittlung der Spannungen verwendet man am einfachsten die zwischen Spannungen 
und Verschiebungen in krummlinigen Koordinaten bestehenden Gleichungen, die wohl erst- 
mals von Borehardt*‘) angegeben wurden. In der hier gewählten Schreibweise ergibt sich 


—r Ai. en. u h e usw (25) 
9 kn Be - . pP S a E i . ’ . _* 
“ hu \du HK, dv hr da Ba —1” 

h.o /V h.„o (U\\ 
ru» G = | et et 
me (noun, h 0» n,)) 


s) Borehardt: J. f. Math. (Crelle), Bd. 76, 1873. 









































INLeWw. 
Mech. 





vier 
ı be- 


208g) 


und 


20h) 


rden 
har- 
ver- 
ela- 

MO- 


I ZU- 
ınde 


ıten- 
chen 
nten 


(21). 


(22 
Igen 
(23) 


die 
man 


(24). 


Igen 
arst- 
sich 


(25) 


26). 
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Der Laplacesche Operator geht über in 
| (0 /h,hu ® 0 (*w h, eo | OQ (® h. 0 \ 


h,h,h"lWdu\ h, 0 = or\ h. or h. oO m/J 


Nachdem nunmehr alle mathematischen Vorbedineungen gegeben sind, soll das Ver- 
[ahren auf den Hohlkegel mit beliebige eerichteter Kınzellast Anwendung finden. 


ba) 
Om 


5. Der Hohlkegel mit beliebig gerichteter Einzellast. Ks werden zweckmäßig Polar- 
koordinaten mit 


e—ncost, yausınr cosm, ze esinesinm . Er - | 
‚uerunde gelegt. Aus den Gl. (25) erhält man 
h.=1, h,=Wu, 2 ee ee 


Wie aus Abb. I und 2 ersichtlich ist, fällt die Kegelachse mit der X-Achse zusammen: ferner 
ist die Kegelspitze zugleich Koordinatenursprung. Der Hohlkegel wird begrenzt von den 
Flächen »—=y und » =. Als erste Gruppe der Randbedingungen ergeben sich aus der Last- 
(reiheit der Oberfläche: 


Für ® Y\ 
0% 0, ‚ FpR vd, Tod "a De Ge Ge MI. 
und 6] 


also sechs Bedingungen. Weitere Bedingungen ergeben sich daraus, daß für einen beliebigen, 
durch den Hohlkegel gelegten Schnitt die an der Schnittfläche angreifenden Spannungen zu- 
sammen mit der Einzelkraft ein Gleiehgewichtssystem bilden müssen. Es muß also die 
.-Komponente der Einzelkraft entgegengesetzt gleich der Summe aller «-Komponenten der an 
der Schnittfläche angreifenden Spannungen sein. Das Entsprechende gilt für die y- und 
:-Komponenten. Wird der Schnitt längs einer Kugeltläche gelegt und sind P,, P,, P. die 
Komponenten der Einzelkraft, so folgen 


Y 


Pr \ run cos (x, u) dF DB er Er et tt 
F u=zu er,w 


Hierbei wurde der Einfachheit halber o,= r,,„ gesetzt. Weitere Beziehungen ergeben sich 
aus dem Gleichgewicht gegen Verdrehen. Es muß für irgendeinen Bezugspunkt das Moment 
der Einzelkraft dem Moment der an der Schnittfläche angreifenden Spannungen entgegen- 
vesetzt gleich sein. Wird die Spitze des Kegels als Bezugspunkt gewählt und sind M,, M,,, 
M. die Momente der Einzelkraft um die X-, Y- bzw. Z-Achse, so wird 


\ 


M,.= \ B. Tun (my Cos (2, 1) 2 cos (9, 1)) dF usw. Win. m,» 77 


F u=zu,v,w 


Indem nun verlangt wird, daß jeweils nur eines dieser sechs Integrale von 
Null verschieden sein soll, zerfällt die allgemeine Lösung des Problems 
zunächst in sechs voneinander verschiedene Partikularlösungen. Da die 
Lösungen für P, und P, und ebenso für M, und M, durch Vertauschung der Y- und Z-Achse 
ineinander übergehen, reduziert sich ihre Zahl auf vier. Da ferner die Lösung für reine Torsion 
(M.==0) bereits von A. FöppI®) mit Hilfe der "Theorie der 'Torsion runder Stäbe von ver- 
änderlichem Querschnitt gegeben wurde, bleiben nur mehr drei Partikularlösungen anzugeben. 





A. Die Einzelkraft wirkt in Richtung der A-Achse | 
und greift in der Kegelspitze an (s. Abb. 1). 


In diesem Falle müssen alle Integrale mit Ausnahme 
desjenigen für P, verschwinden. Der Ausdruck für P, 
lautet ausführlich: 

27nY 


P.=u”\ \(o,ecos® — T„„Sine)sine de dw . (32). | 
00 Abb. 1. ! 





», A. u. L. Föppl: Drang und Zwang, II. Bd., 2. Aufl., S. 108, München u. Berlin 198, 
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jeim Aufsuchen geeigneter harmonischer Funktionen, die diesen Bedingungen genügen, 


ge- 
langt man zu folgendem Ansatz: 
PD, Atln u —+1In (1 -+ eos r)) + CO (Inu + Insin r). | 
| er ee 
8 =B:.—, 5.0. PD, 0 | 
u 5 i 
Mit Anwendung der Gl. (13) erhält man als Spannungsfunktion 
F=(A+C)Inu+Aln(l+cose)+Cihnsine+Beosr .. 0, 5 


Hieraus sind zunächst entsprechend den Gl, (24) die Verschiebungen zu bilden. Aus den Ver- 
schiebungen erzeben sich schließlich mit Hilfe der Gl. (25) und (26) folgende Spannungen: 





| 
Ou [4 (2+@) Beost +(] 
u 
E | ‚_ cos NR | 
Or >| 4 ia Ib cöstv - eot’r 
ae l+-cosr r | 
| l 3 am 
Ow 2 | A ta IBecosr U | DE er >} 
ur | l+cosr sın” " 
sin ® | | cos? 
Tuv 9 | — la I)B+0-— = 
u” | | +-cos v sın® v| 
Tow —=U), T m 0) 


Man erkennt, daß zwischen o, und r,, folgende Beziehung besteht: 
Ö 


Ennadı EEE»: 4 2 ı : AUSPCHEEEEHEEEeE EEE | 


Die sechs Randbedingungen (29) reduzieren sich daher auf zwei. Aus beiden Gleichungen 
folgen: 

A 1 + cos y cos d B | 

bi (1 cos y) (1 COS od) E (a I) cosy)(l COS No) 


Zur Bestimmung von € steht noch Gl. (32) zur Verfügung. Es ergibt sich 


P.(a I) (1 cosy) (Il COS O) 
C hei A 


Ir (coso cos y) |eos” y : cos’) + (2 a) COS y COS 0] 


Die maxımale Spannung tritt an der Innenseite auf ("= 9, s. Abb. 1). Hier wird 





h coty z 
" a ‘. Wir erhalten 
sINnd cos od 
| P,.(eosy + cos0ö) |Bcosd (a I) cosy] 39) 
man ») 5 I? NG? “a, Ipne? ne? | ») OS v COS N] (. . 
ana ”) COS y |eos’ ; rCos + (2 a) cos y cos oO] 
Für den Vollkeeel ergibt sich 
P.(l-+cosy)[ (a Deosy] 40) 
Ö, “ Y r > . . 
max 2a” cos’ y 1 - (2 a) cos y + cos’ y) 











Die Lösung stimmt für den Vollkegel mit einer anderen Lösung für das achsial gedrückte 
IIyperboloid überein, wenn letzteres in großer Entfernung vom engsten Querschnitt betrachtet 
wird '®), 








Für die Kegelschale') ergibt sich dureh Grenzübergang y>6d und a>b eine ein- 
.) 
i . k ö „. 2ah 
fache Formel. Ist Ah die Wandstärke der Schale, so wird a?’ - b? en und wir erhalten 
OS v 














BD 


2rtahcosy 





(41). 





j) 
O max 





I, H. Neuber: Beiträge für den achssymmetrischen Spannungszustand. Diss., München 193, S. 37 u. f. 


11) Die Formeln für die Kegelschale gelten nur für nieht aHzu kleine Wandstärken, da dann die Voraussetzung 
kleiner Verschiebungen nieht mehr gilt. 
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B. Die Einzellast wirkt in Richtung der Y-Achse und greift in der Kegelspitze an 
(s. Abb. 2). 


In diesem Falle müssen alle Integrale mit Ausnahme desjenigen für P, verschwinden. 
Letzteres lautet ausführlich 
2% 7 
I 


P,= u” \ \[o„sin®cosw + 7,,c0osvcosw —- T,„Ssinw|lsinededw] . . . (42). 
0 ( 


Für die Erfüllung dieser Bedingangen ist der. folgende Ansatz geeignet: 





| 
| sın ® sin ® 
PP, =cosw|B H-D 
I-+ cost | cosr 

cos # | sin ® .-sın® 

D, ie; HE (43) 
« |, 1+cosr l -cose| 
A 
D, Pi 
» m 





Die Spannungsfunktion ergibt sich entsprechend Gl. (13) zu 
| B—C D+E] 
. . | a R > J | 
F=coswsmv|A+ ( k- r Es 
| lI+cosv 1-—cos®| 


Mit Hilfe der Gl. (24), (25) und (26) erhalten wir hieraus die Spannungen 





sinn cosm| I_CıH >al 2aE |] 
o s (ö-ra)l —-. +E&) | 
" u Pe ) I|+cosr 1 cost 
sin m cos | IL0_EH ZaUl 2a h--Ü D+E | 
Ö, 2 (a ‚(A—+ Y) r- t - 4 
u” | |-+-cosv ' 1 cos»  (l-+cosr)” (1 - cosr)) 
sin ? cos u ILC_E B—C D+E 
Fa - (a ‚A-+ 1) 7 m z, 
u | (I+cose)" (1 cos)? 
(dla). 
cos m | N B-(i+a)C D-+-(l+a E 
Fe — (a—1)(— A U’+-E)cost +2a(Ü-+kE)- 
u | l+cosv | vos 
sinesinw| af a E B—EÜ D-+E 
Tan 7 N m r . ol 
u° '1+cosr 1-—cosr (l-+cosr)’ (1 --coser)?’ 
sin m | er B+(l+aJC D-+(l-+a)E 
tu = — ka— 1) A+ll+u)(-C+E + | 
| lI+cosr | cos r 


Es zeigt sich, daß wiederum zwischen den an den Randbedingungen (29) beteiligten Spannungen 
eine Beziehung besteht, und zwar wird in diesem Falle 


cos tv -r,„simnv - n„cotw=0 . . .: .: 2 2 2 a 2 2». + Kb): 


Mithin entsprechen den Bedingungen (29) nur vier Gleichungen. Zusammen mit Gl. (42) 
stehen für die fünf noch unbekannten Konstanten daher fünf Gleichungen zur Verfügung, so 
daß letztere eindeutig festgelegt sind. 

Zur Abkürzung wollen wir setzen 


cosy=ec, cosd=d, N=(l-+e)(i+d|id eo +le Mil ed) (46). 


Dann erhalten wir: 


A 1 


(! (a 1) „Id e)” T (2 ayedi(] ed)] 


RB (1i+e)(1+d) 
3 | (a I) (d ey” (a l)(e + d)ıl ed) Jaecd(] ed)| 


e N 

D d—- 90 —A 

(' \ | (a- I)(d - eV’ +la VDle+M( ed) -2acd(l -cd)] (46a). 
E —- l— A 

(' V | (d c”’+(c+d) (1 ed)| 


P,(a—1l)N 
Iald e)A(d ce)” [(3 ed) (1 ed) (dd e)?] I (4 aed(l e d)°’, 
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Die maximale Biegungsspannung tritt an der Stelle w=180°, »—y (vgl. Abb. 2) 
a e a . > 

auf. Hier wird « .  . Es ergibt sich 

sın y 

P,sıny(l e)\t3ld ee? +d|3e - (a all cd) 


all e)itd | edyıl ed) (d e)’| +14 a)edil cd)”, ) 


Die maximale Schubspannung tritt an der Stelle » =%W*", »—=6 (vel. Abb. 1) auf. 
’ ; dc e. 
Hier wird :  . Wir erhalten . 
dsın y 
P,(a — 1)dı(ı1 — ce?) (1 — cd) 


’ . - ‘ S 
mas are td era ed) ed) (de?) +Ht—a)ed(l cd) vn 





Abh. 9. 


Beim Vollkegel müssen von vornherein die Konstanten D und E gleich Null gesetzt 
werden, da die zugehörigen Funktionen in der Achse Unendliehkeitsstellen besitzen. Man 
erhält so 

P,4 -a)siny(l-+e) 


( 


(49), 


’ 
MAX 


tal ec) [2 4. (3 a) e] 


P,(a In(Il-+e) - 
.) 2.19 5) (U). 
2nra?c2 +& -—-o)e] 

Setzt man andererseits in den Gl. (47) und (48) d 1, so erhält man die Formeln für 
den Vollkegel mit einer feinen axialen Bohrung. o,., ändert sich nicht, jedoch 
erhält man für 7,,.x: 


P,(a- 1) (I-+e) 


mar 7 a c|2+2 Je] ' 


(1), 


d.h. genau den doppelten Wert. Eine feine axiale Bohrung erhöht mithin die 
Schubspannung um 100°, 

Ferner lassen sich auch hier durch Grenzübergang die Spannungen für die Kegel- 
schale angeben, und zwar wird 


P, 


mr; ST 4 
zahsıny 


C. Die Einzellast wirkt in Richtung der Y-Achse. Ihr Angriffspunkt liegt auf der 
N-Achse im Unendlichen (reine Biegung). 
In diesem Falle müssen alle Integrale mit Ausnahme desjenigen für M. verschwinden. 
Der Ausdruck für M. lautet ausführlich: 
M,— u | | (tu» COoSw— T„weosvsinw)sne dvdw . . . ....(8). 
u ) 


Beim Aufsuchen geeigneter Funktionen gelangt man zu folgendem Ansatz: 








| 
cos m sın ® sın ® eh 
D, IB -D 5 
u“ | 1-+cost Il cost sın ® 
cosw|, . E " ” 
D, z L sın nr E(sin "Intg . -ecotr]|, re 5° 
u | u 
I | Ai " | 
D, ji Acost + E [cos elntg-, +1 ) - et 
| ur 





angew, 
Mech. 





‚bh. 2) 


(47). 


auf. 


(45). 


setzt 


Man 
(49). 
(0). 

für 


loch 


51). 


die 


ler 


len. 


ww‘ 

.. 

ww 
— 


4). 
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Mit Anwendung der Gl. (13) wird die Spannungsfunktion 


cos w sin ® | 1 | B I | 
(A+0)cosr +2Elcosteintg . +1)-+ 


f E r | ° (>), 
7 2 I+cosr 1 cost) 


Bei der Ausreehnung der Spannungen mit Hilfe der Gl. (24, (25) und (26) zeigt sich in 
diesem Falle, daß zwischen 6, und r,,. folgende Beziehung besteht: 


nt d . . :- 2: ss 2 sr 31: 8. 


Die sechs Randbedingungen (29) reduzieren sich daher wieder auf vier. Als weitere Be- 
dingung kommt noch Gl. (53) hinzu. Die Auflösung dieses Integrals bietet auch hier keine 
Schwierigkeiten. Es stehen mithin wieder für die fünf Konstanten gerade fünf Gleichungen 
zur Verfügung. Die Ausrechnung sei hier übergangen. Es sollen nun die Maximal- 
spannungen angegeben werden. Zur Abkürzung wollen wir setzen 

) 


u 
cosy=c, cosö=d, Intg>=g, Intgz 2 en ME 


ferner 
tall+ed)(l -ed),, TORE 2 » 
3 NT 13) 6(1 ed) + (4 a)ed(2 Hed)) 
| 
+7 (1 e@) bal+ed) + a) cd(IS+-Datacd) 
(I). 


4 ] 12 (4 a)ed(l ed) (2 at2ced) (4 a)(2-ta)ed(e dı? 
Ü ( «> ; 





+-8(4 -a)le - Mt -ed)fi+Hl+a)ed cd} 12ale dl eA)d d?)\ 

) 

Die maximale Biegungsspannung, die wieder an der Stelle = 150°, # = y (vgl. Abb. 2) 
auftritt, ergibt sich zu 

M.(A ce)” f | ce d’ 


l 


; - r “(12 -+5D5a+ta’)e (1 «?) 
mas ed ce) )Ni (A ce?) (1 dl?) Pe | 


-ra(V)  a)e+2a(l a) d +4 (48 Illa- )d(i+cd) 210 dd 
(9). 


+12 (1 ed)” = >(  a)ed? (-+DdDa)e all a)d 





+4(1 a)y(S-ta)d(e d)iq I, 


Die maximale Schubspannung (wieder an der Stelle vw. %°, » =& (vel. Abh. 2). wird 


M.siny (I ce) dA’ na 2 (1 a) —+ (4 a)ed 6. 
tal e dr) ie m .- 


T max 7 a’ (d ’ e) e? N \ I | d? 


+4a9+/ll+a)ed(d - c)+4 (4 a)ea+bed (6-+ta)ed 


(GO). 


-12( ed) (1 de’) all ale Sad-+-2(t > a)e’d 


d 





+4(1 -a)(S+ta)el(e —- d)(g h), 


Beim Vollkegel müssen von vornherein die Konstanten D und E gleich Null gesetzt 
werden, da die zugehörigen Funktionen in der Achse unendlich werden. Es werden 


3M, d+o?[d0 -a)e+2( a] 


) | 6 
Oman ar |6 (Al e)+(4—a)e(2+e))] (sh), 
9 M.(1 + e)’siny ai 
gi; rn — - (62). 
na “zarcel6(l ce) +(4  a)e(2+e)] 


Für den Vollkegel mit einer feinen axialen Bohrung wird (aus Gl. (60) mit 
d 1) 


IS M.(1 + e)’sın y 
IIAX 7 a” ce" [6 (1 e?) +(4 a) e (2 ce)” 


(d.h. Erhöhung um 100 "/,). 


14" 
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Für die Kegelschale wird durch Grenzübergang y—6: 


‘ 


ha 8. | 63 
N e=> 2 2 s sin? y A 
und wir erhalten als Spannungen 
M. M,tg y 64 
"max" 7a? hcos y _. zarh (04) 


Durch Überlagerung mit dem Spannungszustand B können wir jetzt auch die Formeln 
für den Kegelstumpf angeben (s. Abb. 3). Zu diesem Zweck zerlegen wir die Kraft P, 


On 








F% 1 ul “ 
r-?! ’ on 
h; 
Abb. 3 


die senkrecht zur Achse im Abstand a-coty LI (l seı die Länge des Kegelstumpfes) von der 
Spitze angreifen soll, in eine in der Spitze angreifende Kraft gleicher Größe, die mit P,, zu 
identifizieren ist, und in ein Biegungsmoment vom Betrage  Pfa-coty 1), welches gleich M. 
zu setzen Ist, 


Für die Kegelstumpfschale wird dann 


P-I P 


I 
nn | y 
a” h cos y i ah m 


1 
’ max 


(62). 


Schließlich seien hier vor allem die Formeln für sich schwach verjüngende 


Wellen angegeben (Vollkegelstumpf für kleines y), die sich mit e = 1- — tg? y bei Vernach- 
lässigung höherer Potenzen von tg y ergeben: 
| 2 I 
/ 24 P(3 4 ey) 
Pi I m da m Be (66) 
oO = "ur ER . R ID), 
: 2 er I\ I Sm Er. 
31 2 alt ! — a’ 
m) m 


Andererseits lassen sich die genauen Werte für den Hohlzylinder, ermitteln, indem wir 





a” bh? 
e=1-,.,,d=1-,,, und «==» setzen. Die maximale Schubspannung wird 
2 u‘ Zu’ 
| 9\ .) | 
er (3+ \a+li \d>| 
| mi m | Br 
Tax = | a a ie klar FR 20 4 Sn Gr 


(1 — „.)(a’ h*) 


Für den Kegel bei zusammengesetzter Bıiegung ergeben sich mithin 
erößere Biegungsspannungen, aber kleinere Schubspannungen als beim 
Zylinder. Dies erklärt sich in einfacher Weise dadurch, daß beim Kegel die Biegungs- 
spannungen zugleich Querkraftkomponenten besitzen. Die maximale Schubspannung erhöht 
sich durch eine feine axiale Bohrung um 100%, 


6. Zusammenfassung. Ein neuer Ansatz ermöglicht es, die drei Komponenten des 
elastischen Verschiebungsvektors allein durch Differentiieren aus vier harmonischen Funk- 
tionen abzuleiten Das zugehörige Gleichungssvstem läßt sich in einfacher Weise auf 
krummlinige Koordinatensysteme übertragen. Es stellt auch dann noch die allgemeine Lösung 
des elastischen Zustandes dar, wenn eine der vier Funktionen gleich Null gesetzt wird. 
Als Beispiel wird die Spannungsverteilung im Hohlkegel unter Einzellast angegeben. 365 











(63) 


(64). 
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Theory of Vibration of Buildings During Farthquake. 
Bv M. Biot in Louvain (Belgien). 
Gebäudeschwingungen bei Erdbeben. 


Zusammenfassung. Die Bewegungen eines Gebäudes, das von einem Erdbeben getroffen 
wird, haben den Charakter vorübergehender Schwingungen. Die Untersuchung derartiger 
Vorgänge nach dem Heaviside-Verfahren ist den Elektro-Ingenieuren geläufig. Für den 
Fall der Erdbebenschwingungen wird indessen die Anwendung dieser Methode wegen der 
Unregelmäßigkeit der Erdbewegung zu mühsam. Außerdem braucht man nicht die Bewegung 
eines Gebäudes genau zu kennen, sondern nur die Scheitelwerte der Beanspruchung. 

Die folgende Arbeit enthält eine Theorie der Schwingungen eines Systems unter einem 
vorübergehenden Anstoß. Der erste Paragraph enthält die theoretischen Grundlagen. Der 
Hauptgedanke ist die spektrale Zerleguug des Seismogramms, d.h. es wird die Frequenz- 
verteilung des Erdstoßes untersucht. Dann ist es ziemlich einfach, die Wirkung eines solchen 
Stoßes auf ein Gebäude zu ermitteln. 

Im zweiten Paragraphen wird diese Methode auf die Untersuchung der Erdbeben- 
sicherheit eines Gebäudes angewendet, das einen elastischen ersten Stock hat. Eine Tabelle 
ermöglicht die Berechnung der Schwingungsperioden des Gebäudes. Als Endergebnis erhält 
man die durch die Elastizität des ersten Stockwerks hervorgerufene Verringerung der Be- 
anspruchungen., 


1. General Outlines of the Theory ').. We shall confine ourselves to general outlines of the 
mathematical theory and we will neglect the damping. A more detailed account of the 
theory is given in the appendix. 

The analysis starts from the fact that any vibration of an elastie undamped system 
may always be considered as a superposition of harmonies. Although the method is valıd 
for any type of motion, shear, bending, or torsion, let us only consider for the present the 
horizontal vibration of a building. This building, like any elastie system, has a certain 
number of so called „normal modes“ of vibration, and to each of them corresponds a certain 
frequeney. When a building vibrates in a normal mode, all the displacements have the 
same phase, i. e. they all reach the maximum at the same moment. The „shape“ of the 
oscillation is that of standing waves, and the higher the frequeney the higher the number of 
these waves. We will show that any motion can be caleulated when we know these 
modes of vibration. 

Consider the ith floor of mass m;. The amplitude of osecillation of that Nloor in the mode 
of order k is called A, y;;. These quantities contain an arbitrary constant A,, corresponding 
to the fact that a free oscillation has an arbitrary amplitude. 

In view of the application of the theory we have to give to these amplitudes certain 
values and determine these constants so as to satisfy the condition 

i 


Im; A,” Y'ik =]. 


The amplitudes become then 


M; 1. Y;r = - 
ik An Nik | hr m; 
The amplitudes «;,; are then said to be „normalized“. We suppose that we also know the 
frequencies ”"% corresponding ee ee 
The next step is to ealeulate the statical defleetion of the building under a constant 
horizontal acceleration j, that has been applied very slowly. In other words, we have to 
ealeulate the deformation that the building would take if the gravity had an horizontal com- 
ponent j,. Im general, this deformation will be very eomplieated and will be a combined 
ee ee developed ıs valid for any general case, but 
we will restriet our considerations to the study of horizontal defleetions. 
By applying the so ealled property of orthogonalıtv »f the defleetions Hin, namely, 
i 


Y 


m; U; U;; UV h k 


I) See also: M. Biot, Theory of elastie systems vibrating under transient impulse with an application to earth 
quake-proof buildings, Proe. Nat. Acad. Se. Vol. 19, No. 2, pp. 262-268, Feb. 1933. — This work was done during the 
author’s stay as Research Fellow at the California Institute of Technology. 
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it can be proved that the statical defleetion @; produced by the acceleration j, ean be 
expressed as a sum of amplitudes «;; of the normalized free osecillations by the series 


l; 
\'C 
MH; „U;k; 
Z vd] 
where 
(]; In ] 
j 
GC Jo > m; tt;, 


If, instead of a gradual applieation of the acceleration j,, it is suddenly applied at time 0, 
the building is going to vibrate around a new position of equilibrium which is given by the 
previously considered statical defleetion. Ilenee, the motion of each floor is given by 


I; 
’ C, 
M; (f) \ or HM;]; | | COS 7; a . 


This shows that when an horizontal acceleration is suddenly applied to a building the 
maximum defleetion is twice the statical one. 


From this last formula, we may deduce the motion due to an arbitrary acceleration 
(see Appendix D. The results are the following: 


let an horizontal acceleration j, y' (f) act upon the building between the instants O and T. 
Consider then the funetions of the frequeney vr, 


. 
f, eo) =\y(T)cos2rnvrdı 
u 
r 
f. dv) =\y (Msin2avrdr 
ir 
and F(v)=yf’e)+fr(r). 


After the ınstant T the motion of the building is composed of a series of free osecillations 
each of which gives to the different floors an amplitude®) 

(',. u, 

1 tr 
- önv: Flvr). 

ID 
The funetion Fr) is nothing but the frequeney distribution of the impulse, or its 
speetral distribution. The function 2» Fr) is dimensionless and may be called the 
reduced speetral distribution of the impulse. 

Ilence we have the following theorem: When an arbitrary horizontal accele- 

ration j,y'(#) of finite duration acts upon a building, the resulting vibration 
at theend ofthe impulse is composed ofa series of free osecillations each 


Ä | i kUkR _%. 
of whieh has an amplitude equal to the corresponding term s» ofthe 
VD 


statiecal deformation due to j, and multiplied by the value 2rv,F(r,) of 


the reduced speetral distribution of(yfH for the eorresponding frequeney. 
IE we want to know the amplitudes during the earthquake, we must naturally use a 
funetion Fr) corresponding to an impulse ending at the moment we consider. 
It will be notieed that we consider only the amplitudes of the free oscillations and 
not their phases, This is justified by the fact that we are not interested in the motion 


itself of the building, but merelvy in its maximum amplitude. This maximum is the sum of 


the amplitudes of each separate free oseillation. It will not always be reached because it 
supposes that an instant exists for which all the free oscillations have their maximum de- 
Nleetion simultaneously. However, this maximum will many times be nearly reached in a 
short time, and in anv case it Is the highest possible value. 

Remarks. 1. The method may be used for any type of deformation. For instance 
in the case of torsional deformations. Call M; the torque produced around the axis, at the 
ı'" floor, by a unit horizontal acceleration; then 7, M; will be the torque due to an 
acceleration j,. Ihe normalized angles of rotation of each floor in the torsional normal 


2), See Appendix 1. 
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modes being 9;, and I; being the moments of inertia of the floors, the rotation due to the 
statiecal deformation is 


l; 
\ > 
J; Vi GI; k 


with 
i 


C=jı & M; 9x. 


The eondition of normalization of the ©;, becomes here, 
i 


PERS 


In order to get the amplitudes of rotatıion due to an acceleration we must ealeulate as 
before the expressions 
(j; \ 
‚ I, 2rıv; F (v,.). 
DR 
II. The steadvy state oscillation due to an harmonie horizontal acceleration j, et is 
eiven for the defleetion deformation by 
h 
\ Ch; U;]; 1 
DI Fig ' | (07) | 


DR 


serrvit 


Mi 


where o=2nv. This gives an experimental method of measuring the «;, by produeing 
resonance in a model. 


III. In order to show the main properties of a speetral ceurve F(») we ealeulate its 
value for a very simple case; we will surpose that y(f) is a sine eurve sin? N t of total 
length T. "The corresponding function F(v) has a sharp peak in the vieinity of » N and 
its approximate value is 

I sınz(v -N)T 
2 a(w—N) 


This shows that in general a speetral distribution eurve will show a certain number of peaks 
if there exists certain periods in the seismogram (fig. 2). 





[27] 














Fig. 1. 


According to recent observations?) there seems to exist characteristic frequencies of 
the ground at given locations. These frequencies would be given by the peaks in the 
speetral curve. 

If we possessed a great number of seismogram speetra we could use their envelope 
as a standard speetral curve for the evaluation of the probable maximum effeet on buildings. 
The speetral eurves would be of interest as much to seismologists as to civil engineers. 

The direct computation of such speetra might be tedious but automatie electrical methods 
can be easily imagined such as a photographie record passing in front of a photo-electrie cell 
acting upon a tuned eireuit. 


2. The Eifectiveness of an Elastic First Floor. As an application of the method we will 
caleulate the reduction of stress due to a so called elastie first floor in an earthquake-proof 
building. 


3) Suyehiro, Tokyo. 
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let it be of reetangular shape. We assume 

that the deformation is an horizontal shear as 
| shown in fig. 3. Furthermore, the shearing 
rigidity and the mass of each story are supposed 

FF to be constant from the second floor to the 
| top. The first floor only will have a different 
rigidity. In order to simplify the analysis we 
consider the upper part of the building, i. e. 
N) all the stories except the first one as a conti- 
nuous beam having shearing deformations only. 




















We use the following notations: 


Abb 


x.  coordinate measured downwards from the top as origin 
horizontal defleetion of the building 

h height of the upper part of the building 

M mass of the upper part of the building 

» number of stories of the upper part 


hı height of one story 


N force necessary to displace two conseeutive floors of a unit length relatively 
to each other 


MV 
7 mass density of the equivalent beam 
n 
„= Kh, rigidity coefficient of the equivalent beam 


(ı force necessary to displace the second floor of a unit length relatively to the 
eround 


7 | 
E „, Propagation speed of a shear wave in the beam 
h | | | 
A P time Necessarv for such a wave starting from the hottom of the beam to reach 
the top 
f 
M r 
” 
li s 
" 
7 B* y 


According to the general theory we have to consider the normal modes or free oscillations 
of the equivalent beam, which are given by the equation, 


u u 


tl = nm 0 
, Q “2 0 f* 


This may be written with dimensionless variables, 


0: „ 0" „ 
er. ur 


By putting „== z(£)ei’' we get, 
d? z ” 
= ren. 
de 
The general solution 1s 
z=Acos/i&+BsinÄAE. 
Consider the boundarv conditions, 
AN 0 
er 
Nu 
7) 
3 (ru »—=h. 
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Ssume (r ' a » . 
r as Put R= ratio of the rigidity of the first story to that ofthe others, and a Rn 
. Ä .,. 
ırıng n being the number of these other stories; the boundary eonditions become 
osed 
dz e 
the : 0 E—0 N. 
rent as 
Ai dz, i 
5 e, / ne 1 zZ v0) = | (2). 
de 
INtI1- 
nly. From equation (1) 3 —=0 3 —=AcosiE. 
"rom condition (2) igi=a. 
The roots A, of this equation correspond to the free oseillation frequeneies of the building. 
We ehoose eertain values of a corresponding to certain simple values of AR and n as follows: 
| 
) 
N 
15 1.66 25 5 
10 1.11 1.66 10 
ely \ 
D 0.556 0,834 
The values of /, as a funetion of a are given in the following table: 
the a Io I bes bez by hr 
0 0 Ey 2 37 I an 
0.596 0.68 3.31 6.31 9.48 12.60 19.75 
0,834 0.80 3.38 6.41 9.51 12,02 15,75 
C -  4® - „em 
ch 1.33 0.89 ).4) 6.45 9.54 12,65 15.77 
1.66 1.03 3.08 6,53 3,59 12.69 15.80 
2,50 1.15 1.19 6.69 0.67 12,76 15.85 
3.39 1,23 3.86 6,74 I 15° 12,82 15,01 
5.0 1.32 1.04 6.91 9.90 12.93 16.0 
10,0 1.14 1.30 1,22 10,18 13.20 16.24 
x 7 In In In In Ila 
2 2 2 2 2 2 
ns In fig. 4 are plotted the values of 4,’ where /, = ka - The period T, corresponding 
to 7% IS, 
7 an a 
T.=— 
Al; 


which would oceur ıf the 
the only elastieity being 


It ıs interesting to compare the fundamental period T, to that 7, 
building would be perfeetly rigid from the second floor to the top, 
due to the first floor. We get, 








Au vom M _2Znt, 
M „rGu- dd / In ’ 
dF (i “a 6, RR | Ze 
ann iu | Be, in ie | 
"he ratio of frequencies fy 7 - is a funetion of a. A, | 
i a 12) 
10, | 
z Io r fu 08 X, 
A | 
06 
| 2 
0 2.50 0.725 rn ; | 
dei 3 | 
0,556 0,910 3.35 0,674 
02 ei | 
0.834 0.875 5.0 0,590 / x, %s | 
l ı & 
1,11 0,840 10, 0.455 a IT IETUTTEH 
1.66 0.800 x o | Fig. 4 
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The shape of the normal modes is given by 
tk) = A, cos in 2 
We have to determine the eonstants A, by the condition of normalızation which reduces 
here to the integral 
h 
m\ A) cos’/, dr =] or M A, \cos’i,EdE=1. 


IE u" 


By integrating and using the equation A,tg/, «a 


A, 


3 | | 7ı ’ 


M| «a C08? 7] 
The statical deformation of the beam due to an horizontal acceleration g is then 


h 


\'C 


na) „urlr), 
A Vi 
h qmh ZN 
where U; egm\u,(e)dr —— (|0084,. Ed} 
. | M f a cos? 4;|ı 
) ' 1,2 
= | Ak 
232... 81nd;  @acos/A; 
\ecos 4, 2dE w 
N A] Al; 
er 4, ’ 
Since on | the statical deformation takes the form 
h 
27 he COS Al, a in 
üle)=2gt,. a \ Ä — 4 c08ÄArLE. 
u AaACOS" /ı | 
Ak | I +- 9 
| Al 
2a cos 4, 
Put — + =Bı 
a u OS” /7;| 
Ar |1-r | 


in? 


These coeffieients are dimensionless and funetions of a only. Their values are given in the 
following table: 








a B, B, B, B; 
0 .* () () () 
aA 
0,996 2,336 V.OOSOO 0.U00692 V.OOOLBS 
0.834 1.734 — 0.0116 0.UHO09TO 0.BOO2OS 
1.11 1.125 0.0137 0.00122 0.000272 
1.66 1.001 0.0166 0.00177 V.OOOBSH 
>50 V.I02 0.0191 0.000234 0.000544 
3.33 0.802 0.0203 0.00274 0.VOOHS2 
2, 0.710 0,0210 00334 — 0.000912 
10,0 0,612 0.0212 000404 — (,00125 
N 01,520 0.0192 0.00417 000151 


The statieal deformation is then given by: 
h 
ülx) = gt, 2 BrcosirE. 
By using this result we are able to find the amplitudes of oscillation produced by an 
acceleration gy(#). 
We build the reduced speetral intensity 27» F(r) of the function of time y(#). The 
amplitude of each free oscillation of the beam is then 


(a) gt, Br; cos hr & 2rv,;F(r,). 
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We are also interested in the total shear produced in this beam. The amplitude of 
hear produced by each mode is 


OU; mol, 


S;:=gM-B;4, sin 4, &2rzv,F(v,). 
A r hı ( ze I; 4 ’] I; I; I ( l, 


I: u 


Replaeing 2», by , we get 


FE F(v,,) 
S;; qM B, 4," sın A, & / Re. 


For the fundamental oseillation the shear is maximum at the bottom (E== 1) 


F(v,) 


S=gqgMB,4’sin4, 


0 

For the osecillations of higher order, the maxımum value is 
.„Fo;,) 

Sı: —( M B; Alk } 


Put B.i21’sni =UC.fe), B; 72° 0:(a). 


The maximum shear corresponding to each free oseillation ıs then given by 
F(v),) 
 ° 
The coeffieients C),(a) are dimensionless and functions of a only. Their values are 
plotted in fig. 5 and given in the following table: | 


S; = 9M C,(a) 




















a C, C, C, C, 
0) 0 0 v 0 07 i 

0.556 - 0,340 0,0486 — 0,0139 0.00622 
0.834 0,397 0.0666 — 0,0199 0.000942 
1.11 - 0,438 V.OS1S — 0,0254 0.0123 
1.66 0.496 0.107 0.0377 0.0177 
2.50 -0,544 0.132 0.0517 0.0254 
3.33 — 0,572 0.150 0.0623 0.0324 
5.0 0.600 0.172 0.0812 0.0451 
10.0 0.629 0.196 0.105 0.071 

X 0.642 0,212 0.129 0,091 





We see that the fundamental oseillation is by far the most dangerous, and that, for 
siven values ofthe »’,s, the influence ofan elastie first floor is important 
only for values of a smaller than 3. For example, it would not be in general of great 
advantage for a seven story building to build an elastie first floor having only 50°, of the 
rigidity of the other floors. 


0 


Appendix 1. 

Vibrating systems under transient impulse. The following theory gives a method of 
evaluating the action of very random impulses on vibrating systems (1. e. effeet of staties on 
radio-cireuits or earthquakes on buildings). In the following text, we will use the language 
of meechanıies. 

Consider a one-dimensional continuous elastie system without damping. The free 
oseillations are given by the solutions of the homogeneous integral equation, 

h 

yow\o(S)alrg)y(S)dE. 

da 
Due to the nature of the kern there exists an infinite number of characteristic values 
; of» and of characteristie funetions ; solutions of this equation. These functions give 
the shape of the free oseillations of the system. "They are orthogonal and have an arbitrary 
amplitude. This amplitude may be fixed by the condition of normalization, 

h 


KOrMORFEEE 


da 
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We now suppose that certain external forces f(x) are acting on the system, these forces 
being expressed in such a way that the produet of the displacement y by f(x) represents 
the work done by this force. For example, if f(=) is a moment y will be the angle of 
rotation around the same axis as the moment at that point. 
It can be easily proved that the statical deflection of the system is, 
y Hi 
ee te Eee ie a 
(); 
ö fir) 
| 


! 
in a series ofthe orthogonal 
ol) 


where €; ıs the Fourier eoefficient of the development o 


funetions ;, 

f(x) 2 
‚=20;y, 
Ol) 


hence 
h 


C; \ f(& ; (&) d & . 
u 

If the applied forces are variable with time and harmonie of the form f (x) e ® ! the defleetion 
is expressed by the expansion, 

CH | 


(N); (ı) 


"ı-(@) 


The amplitude is composed of each of the terms of the statical deformation (l) 
multiplied by a resonance factor. 

The motion due to a sudden application of the forces is of the same type, and can be 
deduced immediately from the preceding harmonie solution. 

By using Heaviside's expansion theorem we get, 


\' ( i UF 
Ja = 2 


„; 


The amplitude due to a sudden applied force f(r) is composed ofa series 
ofoseillationseach ofwhich has an amplitude equalto twice the correspond:- 
ınz term of the statical deformation (Il). 

We will now investigate the action of varving forces of the type fr) y (t); these forces 
are supposed to start their action at the origin of time and to keep on during a finite time T'. 

Using the Heaviside method, and the indieial admittance (3), the motion after the 
ımpulse has disappeared IS eiven by, 


1-08 will - . » : 2 2 m 2 2 2 2 a + MB 


r 
'd | 
HM \ 79a (t r)y(T dr. 


mn pe ; e 1 
Mn ".. sin; flo; \y(r)coso;rdr] coso;t|o;\ v(r)sino;rdrl}. 
} 0 u 


This motion is the superposition of free oseillations. Their respective amplitudes can be 
physicallv interpreted as follows: 


. 
Put fo )=\y(N)cos2rnrtrdr, 
u 
r 
f. (2) = \y(r)sin2avrtrdr, 
u 
i ; 2) 
where v ıs the. frequeney =, 


The amplitudes of composing free oscillations may then be written, 


('; 9; ri . 
2a yo )+tfron). 


(Dj 


Now, according to the Fourier Integral, 


v(d=2\f,(w)cos2rnvtdvr+2\f,w)sin2nvidv . . 2.2... 0). 
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Irces This shows that the expression 
ents v .o \ ’»= 
e of Fo)=yW +’) 
nay be considered as the „speetral intensity“ curve of the impulse. 
The amplitude of each free oseillation due to the transient impulse is, 
(1), CO; 9; 
FE ee ee a A 
VD; 
onal i i 
ıe expression 2» F(r) is a dimensionless qu: ıat we will call „reducec 
TI N 2 F(») is a dimensionless quantıty that Il call | l 
speetral intensity“. Wethen have the following theorem: When a transient impulse 
ıcets upon an undamped elastie system, the final motion results from the 
superposition of free oscillations each of which has an amplitude equal to 
er 
the eorresponding term of the statical deformation (Il) multiplied by 
(VD; e 
the value ofthe reduced speetral intensity for the eorresponding frequeney. 
tion The advantage of this theorem is that for the ealeulation of the motion it replaces a 
complieated impulse by a speectral distribution which is always an analytical funetion of the 
(2) [requency. 
ä This theorem could also have been established by starting from (2) and using direetly 
the Fourier integral. We will apply this last method in order to generalize the theorem to 
ı (1) the case of an elastie system with viscous damping whose motion is defined by the equation 
“urayterymiluvl) . 3 er eriir s 
) hi 
The impulse is supposed to be given as before by the speetral distribution (4). Introdueing 
a complex speectral distribution, 
(3). ?6b)=hle) if.) 
we may write, 
ies Pen 
n.d- vd)=/jy Bea rn ZN, 
FL 
= where according to the Fourier integral 
F \ y 
the nis 
q (v) == | y' (r)e -RFT AT. 
" 
The funetion 9 (r) is holomorphie; its expansion in a power series is 
] gr 
N 
Ypir) iu 
n! 
where 
e 
14,6 2a)" ty dat. 
be 0" 
Calling M the largest value of |y (f) and I (vr) the coefficient ofi in the variable » considered 
from now on in the complex plane, we have 
M 2 
( (v) er p> zzIt(v)7 l i 
> 1n)| | 
This shows that for > T, g(») e®='”t has an upper limit. 
Consider now the elastic system (6) under an harmonie impulse A e2=irt, the correspond- 
Ing motion Is 
Aerivt| | | ] 
"JU 2er TE TREE DEGETETE 
smTpm vr —rv, vv v.| 
Ihe quantities »,, and v», are complex frequencies 
», =ai+ß, „=ai—ß. 
The free oseillation of the system is damped and given by 
4). 


E znatcos2nßt. 
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Aceording to (7) and (5) the motion due to the impulse Ay (f) will be, 


F 


dv. 
; v,| 


.] ne 008 
\y()ezirt 


. p ıy N 
JL /) nt KL v, 


y(#): 5 


F 


We have seen that (vr) e2’"‘vt has an upper limit, and bv using then the method of eontour 
integrals and residues, we find 


Ari 


y(t)= lq (v)zrirt 7 (v,) ee Tivat], 


t pm 


At the time T when the impulse has ceased, the amplitude is 
y(h) == — - end Q (ai | pP) p 2naT 


A 
I’ m(p—+ ar) 
This last result generalizes formula (5) to the ease of damping. We have 
to consider a complex frequeney ai+p and the analvtical prolongation 
y(at+p) ofthe speetral distribution ge) ofy(d). 


The quantıty is the defleetion for the statie deformation due to a force A. 


Appendix Il. 


A method for computing free oscillations in case of non uniform mass and rigidity. 
Whenever we know the shape 7; of the free oseillations of the building, we can caleulate 
the corresponding frequencies or 7; by the following method: 

Consider the building as a eontinuous shear-beam and oscillating freely with the fre- 
queney of order k. When the amplitude is maximum the kinetie energy is equal to zero 
and the potential energy is 

Ti Iyn\® 
et 1] 1;\" 
via, 
3 \. - | dr | 
0 
On the other hand, the potential energy passes through zero when the strain disappears; 
at that moment the kinetie energy is maximum and has the value 
h 
I ( | 
3 \ me) dr Yr ar, 


0 
Kuquating those two expressions we get the value 


hi 
| ur) | 
M 


7 
- 


d un | ” pr 


d sE 
mx h 


Im(a)y da 


which is independent of any arbitrarv eonstant multiplying 7. 
The orthogonal funetions 77, representing the shape ofthe free oseillations may be found 
bv the following semi-empiriecal method. We note that those funetions are defined by equation. 


d d Un 


I 
udl.r) 
; d.r 


4 m (ar) en” 97 0, 
dr! ' 


which by a change of the independent variable 


(da 
u Jul) 
becomes 

dÜ 2 

\|d 2? lv. (2)] + m [er (Z)| eo’ y—=d. 


A l. (2 
This ıs the equation of buckling under a load P of a beam of moment of inertia 
pp 


E nm 1ı " 


ig) 
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If we eonsider an elastie strip of uniform thickness h, its moment of inertia will have the 
value I under the condition that the variable width A satisfies the equation 


bh h® ] P A 


12 Emuo um’ 


where A is an arbitrary constant. We may choose for z any arbitrary scale eonvenient for 

















our 
we are interested only in the shape of the functions. 
In order to realize the given boundary eonditions, the strip will be repeated symmetrically 
around a point representing the top of the building (fig. 6). The deformation of the half 
> 
ne 
| | 
| yo | | 
u ei u‘ 
A. | | 
on | | 
Raseıza] | | 
Fig. 6. 
ity. ' rn 
ate strip under various loads will give the corresponding functions Y7;,. The first deformation 
only is stable, so that the others will have to be stabilizel by a special but very simple 
bon. device. Knowing the shape of the „7, we then compute the attached frequencies by the 
method above. 
ero 
Appendix II. 

Method for caleulating the acceleration spectrum from a displacement spectrum. Suppose 
that we know the displacement diagram of the ground, in the form = (f), v being a certain 
arbitrarv length., and T an arbitrary length of time. The acceleration of the ground is 

rs! ® 
‚Nr den MEY 
from there, 
23 mo . " 
"HP hei: 
We can easily deduce the speetral distribution of the unknown acceleration diagram (f) 
[rom the speetral distribution of the displacement diagram by the following considerations: 
Consider the speetral funetion G (vr) of y(d 
Go)=yg’o)+g’o), 
m 
y,()=\gy(T)eos2avrdr, 
u 
. 
nd Jr) \gp(r) sm2rrrdr. 
N. u 
From the Fourier integral we get, 
= x 
(= 2\g,() cos? rvrtdvr+2\g,(r)sin2avrtdr, 
0" 0 
a FT 


ve’ (MT 2\(27vrT)g,(v)cos2art dv +2\(2 ar T)g,(v)sin2avtdev. 
0 ( 


This shows that the speetral distribution F(v) of the acceleration diagram y(f) is 
Fo) = (2avrT)G(v). 
This work was undertaken at the suggestion of Prof. Th. von Kärmän. The author 
wishes to express his appreeiation of the continual interest Prof. von Kärmän and 


Prof. R. R. Martel have taken in its progress. 363 
(California Institute of Technology, 1932.) 
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Fehlerabschätzung für das Störmersche Integrationsverfahren. 
Von Günther Schulz in Berlin-Wıilmersdorf. 
(Aus dem Institut für angewandte Mathematik an der Universität Berlin.) 


ur Durehführune der numerischen Integration einer gewöhnlichen Differentialzeleichung 
zweiter Ordnung 


d? Ic v 3 / 
De Ela NM. . or een. 
2. dx . R 
deren rechte Seite die erste Ableitung dt nieht enthält, hat sich ein von ©. Störmer (Oslo) 


angegebenes Verfahren‘), das sich ohne weiteres auch auf Systeme solcher Differential- 
eleiehunzen übertragen läßt, als besonders geeignet erwiesen. Mit ihm haben Prof. Störmer 
und seine Mitarbeiter zur Bestätigung seiner Theorie des Nordlichts in langjähriger Arbeit 
eine eroße Anzahl Bahnkurven von Elektronen im magnetischen Feld der Erde berechnet’). 
Dabei hat sich das Verfahren immer wieder von neuem durch seine leichte, bequeme An- 
wendbarkeit und die Genauigkeit der erhaltenen Ergebnisse bewährt. Seine Bedeutung für 
die praktische Integration der Bewegungsgleichungen der Mechanik, falls die wirkenden Kräfte 
nur von Ort und Zeit abhängen, liegt daher auf der Hand. Die vorliegende Arbeit gibt erst- 
malig eine exakte Abschätzung des bei Verwendung der Störmerschen Integrationsformel 
begangenen Fehlers sowohl für die einzelne Differentialgleiehung (1) als auch für Systeme 
soleher Gleiehungen. Bei Störmer findet sich in den angegebenen Arbeiten nur eine Ab- 
schätzung der Größenordnung des Fehlers für einen Punkt der die Integralkurve approximieren- 
den Punktfolge unter der Voraussetzung, daß die vorhergehenden Punkte fehlerfrei sind’). 

Im Abschnitt 1 wird zunächst Störmers Hauptformel mit Rücksicht auf die folgende 
Fehlerabschätzung auf eine neue Weise abgeleitet und ihre praktische Anwendung kurz ge- 
schildert. Es zeigt sich, daß dieses Verfahren, das ebenso wie das von J. C. Adams ange- 
gebene*), mit dem es im Grundgedanken weitgehend übereinstimmt, sukzessive Extrapolationen 
benutzt. Während aber bei Benutzung des Adamsschen Verfahrens die obige Differential- 
gleichung zweiter Ordnung durch Einführung von z=y in ein System von zwei Differential- 
eleichungen erster Ordnung für die Funktionen # (N und „(= x (f) verwandelt werden muß, 
ist dies bei Verwendung der Störmerschen Formel unnötig, die Berechnung von # (f) wird 
vielmehr erspart, und man erhält unmittelbar Werte für das gesuchte Integral der Differential- 
oleichung (1). Ist ein System von Differentialgleichungen, etwa 


S=fla a, H, BIER . 5: ea en 


zegeben, so fällt diese Ersparnis für die Rechnung noch stärker ins Gewicht. Die in 2 und 3 


zegebene Fehlerabschätzung verwendet mit den erforderlichen Verallgemeinerungen Methoden, 
die R. v. Mises’) bei der Herleitung einer Fehlerformel für das Adamssche Verfahren be- 
nutzt hat (Trennung von (uadratur- und Interpolationsfehler, Angabe einer Differenzengleichung 
für die Fehlerschranken, Majorisierung der Fehlerschranken dureh eine Partikularlösung dieser 
Differenzengleiehung). In 4 werden die erhaltenen Formeln auf Systeme von Differential- 
eleiehungen ausgedehnt, wobei sich eine interessante Anwendung einiger Sätze von G. Fro- 
benius über Matrizen mit positiven Elementen ergibt. Schließlich sind als Beispiel in 5 für 
eine dem Institut von Herrn Prof. Störmer freundlicherweise zur Verfügung gestellte Be- 
rechnung der Bahnkurve eines Elektrons die Ergebnisse der Fehlerabsehätzung nach den hier 
eerebenen Formeln angegeben). 


I1),(C', Stürmer: Archives des seieneces physiques et naturelles, 4. per., tome 24 ıGeneve 1907), p. 5—18s, 113- 
IDS8, 221-947, 317-364, insbesondere p. 221—235. Das Verfahren ist ferner in folgenden Arbeiten von Ü. Störmer 
dargestellt: Skrifter ute. av Videnskaps-Selskapet i Kristiania (Mat.-naturvid. Klasse) 1913, Nr. #4, insbesondere S. |] 
bis 11: Nr. 14, insbesondere S. 41 bis »7. Norsk matem. Tidskr. 3, 1921, S. 121 bis 134. — Comptes rendus du Congres 
intern. des mathem. Strasbourg 1921 (Toulouse 1921), p. 2227. 

Die reehnerisch gefundenen Bahnen sind übrigens kürzlich mit sehr guter Übereinstimmung von E. Brüche 
auch experimentell hergestellt worden. Vel. C. Störmer: Ztsehr. f. Astrophys. I, 1930, S. 237 bis 274, wo auch weitere 
literatur zu finden ist. 

3) (', Störmer: Skrifter utge. av Videnskaps-Selskapet i Kristiania (Mat.-naturvid. Klasse) 1913, Nr. 14, S. 57. 

ı F., Bashforth. ). C. Adams: Anattempt to test the theories of eapillary action, Cambridge IS, p. 15 
1, 43-46, 58-62. Das Adamssche Verfahren ist auch in der in Fußnote °) angeführten Arbeit von R. v. Mises ge 
sehildert. 

5), R.v. Mises: Zeitsehr. f. ang. Math. u. Mech. 10, 1950, S. 81 bis 92. 

6) Die Arbeit von J.Tamarkine: Sur la methode de ©, Störmer pour l’integration approchee des equations 
differentielles ordinaires, Math. Ztsechr. 16, 1923, S. 24 bis 219, befaßt sieh nieht mit dem eigentlichen Störmerschen, 


sondern mit dem Adamsschen Verfahren. 
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°n. 1. Ableitung der Störmerschen Formel. Von der rechten Seite der gegebenen Differential- 


rleichung (1) sei vorausgesetzt, daß sie der Lipschitzbedingung 


fix, t) fr BISRIRB |. . 2 2 2 Hr 2 » 1 5.5 
ung für beliebige x, .*, t eines gewissen Bereiches der «-t-Ebene genüge; dabei ist x eine ge- 


wisse Konstante. Alle Verfahren der numerischen Integration bestehen darin, daß anstatt 
des zu gewissen Anfangswerten gehörigen Integrals .r (f) der Differentialgleichung eine Folge 

(1), von Näherungswerten «,, #3, #3, .. ., die den Werten #,, t,, #f,, ... der unabhängigen Veränder- 
lichen # entsprechen, berechnet werden. Die dabei gemachten Fehler seien 


slo) 

Et 2 de A re ae 
1al- [3 . * “ * * 
Ier Ks sei ferner f,,, — t„=At=h als Schrittweite angenommen. Die Integration sei bereits bis 
mr zur Stelle t=#, durchgeführt, d. h. es seien schon Näherungswerte ..., Won, Ina Un für die 
t?). wahren Funktionswerte ..., (fu...) “(tu _,)» “-(#f„) bekannt, ebenso die begangenen Fehler 
An- 5 Eyoos En 15 &n oder wenigstens obere Schranken ihrer Absolutbeträge '), Die Störmer- 
für sche Formel, die wir jetzt ableiten, ist eine Vorschrift zur Berechnung von =, ; , aus den be- 
fte reits bekannten Werten. Einmalige Integration von (1) liefert 
"st- tı 
nel FU yeee ROT 10; .% © | 4. 5 
me In-ı 
ıb- eine zweimalige Integration ergıbt 
en- 
| lu+ı Y 
ıde (tt, +) —- el) lt Ile  )=)do)fla,ddt. Eee  ° 
(ir tı I 
=c 
B Ersetzen wir in Gl. (6) na durch » — 1 und subtrahieren die so entstehende Gleichung von 


u. Gl. (6) unter Berücksichtigung von Gl. (5), so erhalten wir 
al- 
al- 27 furı G 27 Y 
ıß, # all,ı,) It-\fiae.ddt+ \do\fir,tydt \doe \ fir,Bdt. ...(N. 
rd In-ı ty 27 in-ı In-ı 
Hierbei sind — wie überall in dieser Arbeit die I rückwärtsgenommene Differenzen, 
die dureh 
>) lu, Un Un, ?u,„=Au, an rn Krk 6 
definiert sind und bisweilen auch mit V bezeichnet werden*). Nun ersetzen wir in Gl. (7) 
13 f(x, t) durch das Interpolationspolynom r-ten Grades in # (in der Lagrangeschen Form), das 
D, an den r—+1 Stellen #,, fu 1» ==, fur die Funktionswerte f(x (t,),t,.), Fir tu Eu)» 
no Ä f(x (,_ tur) annimmt, und rechnen die Integrale, die im Integranden Polynome in f ent- 
ng | halten, aus, was gewisse Zahlenkoeffizienten a ergibt. Ersetzen wir ferner auf beiden Seiten 
er alle x (#,) durch «,. so entsteht 
al- 
Tr r r r 
VO - “ Y . u 5 . Y R 9 - ö 
D x. 4 ‚==h' > Ar ofn o+ > A fr ofn 0 R} d: Oo fn oj= h? > Ar ofn 0. . . (9) 
ul D==% Oo 0 o—\ e=V 
e- 
er mit d«ro+aro- a"ro=Aaro. Zur Abkürzung ist hier f(r,,t,) = f, gesetzt worden. Die a’ 
rühren von dem ersten Integral in Gl. (7) her, die a” vom ersten Doppelintegral, die a’”’ vom 
zweiten Doppelintegral. Die erhaltene Näherungsformel ist bereits die Störmersche 
® Formel für den „Näherungsgrad“ r in der Gestalt, wie sie unserer Fehlerabschätzung zugrunde 
er . x . u.4 j ; r R 
.; gelegt wird. Für die praktische Rechnung ordnen wir nicht nach den Funktionswerten f,_r; 
res fun, fu, sondern nach den Differenzen f,, Afu, .--, J”f„ und erhalten so ın 
he u) D) 75 gr 73 7 > ' 7, P » 
re u My ri = h? (3, fi. u. B, If, + Bs 5 A 1 ..0. 4 Dr j’ Fu . . . . . . (10) 
dc die Hauptformel, nach der gerechnet wird, in der Gestalt, wie sie von Störmer angegeben 
er worden ist. Die auftretenden Koeffizienten a,. und 30 hängen durch die Beziehungen 
7) Diese Anfangswerte und die Absehätzung ihrer Fehler muß man sieh auf andere Weise verschaffen, vgl. 
. dazu R. v. Mises: a.2.0.8.89 und Günther Schulz: Ztschr. f. ang. Math. u. Mech. 12, 10932, S. 44 bis 59, ins 
- besondere S. 21T. 


s) Man beachte, daß in den Arbeiten von C. Störmer die Jvorwärtsgenommene, durch Jun = Uyzı Hy 
definierte Differenzen bedeuten. 
15 
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PAR Bd, + B,Tr B, t- Bst... + BD, 
1 Fr. 3 ,\_ 
a,, (1 JB. (7 )3. = ("3 
a, (3)3: HS)%+...+[2)3, ee er 


3\_ Pi. 
üys 2 A 3 )Br 





zusammen. ihre Zahlenwerte sind 


Ko L. G;, 0 

13 | | 
1 R A, ur Ayo 
doy [> >| ; >. 1» 

1 „ | | 
Ayo (; OAzı 12°’ Uyg — 29 Ay; — [> 

2) 11 37 2 a“ 19 

A u. Oro . Urn — =. A 
240 “ 15 a 3 5 “240 
> ) ) l > l 7 19 73 , 

Po | , Pı ”, ’» R 1» 3199 I’ 4) Pam Zn 


Die Angabe der Werte von a’, o, a 0, & ro Ist überflüssig. Als Probe hat man 


Die Verwendung der Störmerschen Formel soll der Kürze halber gleich für das 
System (2) unter Annahme des Näherungsgrades r =3 gezeigt werden. 








In | In 1 Finn Ua \ 1° Fa en YIn-a Pyn In-ı 8’ On 
Ir, 3 j IFn- 3 ange EP fa-r Yn-3 A In-3 2 0 
In 3 Laer 3 1" Lu? 2 In 3 : 1’ In 2 a In 3 . Pe In 3 Ps , 
law;  Afa! „. |A’fn-: lu: A 9n-:| .. ’n-ı 
7 2 Una I" ni In 2 ! P In- 1 ut HUn-2 Yu ı In-> 1’, 1 . 
lau. Aal a5 It Iyn-ı| Im-ı I" gq, 
e. ni 1’ In-ı i Ef, —— | 5% Ya IAn-ı Zi Ten 
rn — / NEED | I — | On mm 
bi LT ec Pius; In RE Sen Mn un Pas, In EN 
ee 


Ist die Integration bis zur Stelle f,, zu der die Funktionsnäherungswerte x,„, %„ gehören, be- 
reits vorgeschritten, so können beide Differenzenschemata bis zur stark gezeichneten Linie 
berechnet werden. Formel (10) und die gleichgebaute für I? 9, ,, die auf der rechten Seite 
die g-Diflerenzen statt der f-Differenzen enthält, gestatten dann aus f,„., Ifu. If, P’fns In: 
Ian. "gu. .Pg„ die Werte Pa,4,, 1 9y„., zu berechnen. Die Differenzenschemata können 
sodann durch Berechnung von Aa. +,, An+ı, AYn+ı Yn+ı fn+ Aln+o» Pfan+zs Afnt+ı 
Ins Ins Yun gu, in der angegebenen Reihenfolge vervollständigt werden. Für 
die praktische Reehnung empfiehlt es sich, statt der f,, 9, die Größen &,—=h?’f,, 9,=h’g, 
einzuführen und für diese Werte die Differenzenschemata zu bilden. 


2. Rekursive Abschätzung des Fehlers. Zur Durchführung der Fehlerabschätzung sub- 
traleren wir Gl. (7) von Gl. (9) und addieren die Identität 


Fr r 
v0 Jr? B} (a', 077 4 O0 a; 0) f(x (Fu 0). In. 0) h’ > Ar O fr (th, . 0), I 0) . » (12). 


v (=) 





rew, 
ech. 


N—9 
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e- 
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te 


»n 
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Nach Einführung der Fehler (4) erhalten wir 





r 
3 Hi h?’ >'ar olf o,In—o) f(x (in - o), bi - o)] 
o=V 
r a7 
+h|n Narofl& (ne), tn) fa, Dat 
u, In-ı 
r Ins ß 48). 
+ N rofl (no), tn-)— \do\fle,dat 
= () ty In 
r I Ä 
Ih? > a” ro f(& (tn-o), n-o) -\do\fla,d)dt 
0 0 27 ı I 1 


Die Maxima der Absolutbeträge der letzten drei eckigen Klammern für alle in Betracht 
kommenden Intervalle seien durch e,', e,’, ce,” abgeschätzt; die erste eckige Klammer ın 
Gl. (13) schätzen wir mittels der Lipschitzbedingung ab und erhalten 


r 
1 > ‚ 0) Yı— . 
E l’ en+tılr h’ | Arofn—o + Cr mıt Cr er + 6, t ie . . . (l 1). 


p 0 


[23 [224 


Um einen Ausdruck für e,', e,’, ec, zu gewinnen, beachte man, daß jede der drei letzten 
eekigen Klammern in Gl. (13) den Unterschied zwischen einem Integral über f(« (t),t) und 
dem entsprechenden Integral über das Interpolationspolynom r-ten Grades darstellt, das an 
den Stellen #,, t,_1+ ::., fur dieselben Funktionswerte annimmt. Das Restglied der Inter- 
polationsformel lautet 


TER Dr 1 ABEL) Ba lasse 77 A a DEETR EEE ZZ DEE 


wobei f"+D (r) die (r + 1)-te Ableitung von f (x (ft), $) nach ft an einer Zwischenstelle 7 zwischen 
t, tus Anis =, Eur Ist. Wir können daher setzen 


cr =h Ta Di PH N Bade. . : 
oder, wenn t—t,„—=hu gesetzt wird, 
: hr+3  _ | , i 
er = FED! fer+D» max © | +1l)...(u+r)du| .. 2.2.0.0. (16). 
Analog ergibt sich 
hr+> \ 
u da EEE \ “) vu+1l)...(u+r)du ...2...(6%), 
hr+3 s A | ’ 
= 725 Ch Aue u a  , 


Die Anwendung des Mittelwertsatzes ist statthaft, da die Integranden in den Integrations- 
intervallen ihre Vorzeichen nicht wechseln. Die in den Gl. (16), (16°), (16°) in den Absolut. 
strichen stehenden Integrale bezeichnen wir mit 7,', 1,”, 1” (I’<0, 1," >0, 1,”<V0) und 
geben nachstehend die Zahlenwerte. Wir haben also 














Per; r+1]) 
fl r i s nn ’ | " . _ 
C, BED! 'f Be } mit fi I,|\+\1.,. | +1. |... AM. 
A 1 r 1 2 l l 
r=]: ‚eg, 54 z I, 7, h=7 
l 19 1 
_ Oo. [4 ni. [7 N; [220 an. 
- ,=-- = , = „h 
ie a A ef be: 19 
=: IL: 30 Fu re I, 60 I, 6 
9) 363 41 27 
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Zur Vereinfachune setzen wir in Gl. (14) 


zh? :/Groi = Ao r } . : . . . . : . . i ; j (15) 

und erhalten 
ent den tn-, A,|n|+ 4A, len -ı|+:.::+4,I|u-r|)tcr (n=r,r-+1l,...) (19. 
Die Koeffizienten A,. A I,. e, auf der rechten Seite sind sämtlich nichtnegativ. Als 


EEE 
bekannt müssen wir die Fehler &,, &,, .... &, oder wenigstens obere Schranken ihrer Ab- 
solutbeträge annehmen: es seı 


£ N £ N 


0 u» l 


>. 05 Be 
Mit Hilfe der Ungleichung (19) soll nun eine Gleichung zur rekursiven Berechnung oberer 
Schranken X, ,, Arts, »-. hergeleitet werden, so daß |e,; | eis De 


eilt. Für das Foleende führen wır noch zur Abkürzung ein 


O, r |; =... +18 N \ 


m I | ] . y “Ay, . X, - T N\ 


In «der Uneleiehune (19) setzen wir an Stelle von n die Werte r, r-+-I1l. r-+-2. ...n In 


u >85,  . 0). 


und addieren diese » +1 Ungleiehungen und dazu noch die Ungleichung 


€, €, N i X, r ro N . . . . . . . . . . . (21). 
In der so entstandenen Ungleichung wird die linke Seite verkleinert, wenn wir dort  &,4,—&, 
schreiben, was wir noch weıter zu &,+, &, verkleinern. Wir erhalten schließlich 
lEn+ı Ey + 4, (0, 0,4) 7 A, (0, ’ Ör_s) ... +4, ‚6, +3 O,) 
(22). 
A,0, rt (n r+1l1)e,+X,+4X,_, (n r, r-+1, 9 


Rechts stehen hier ebenfalls nur positive Glieder. Denken wir uns diese Ungleichung zu- 
nächst für »=r hingeschrieben und rechts 2, |, |&, + »:., |&r durch X,, X,, .... X,, also 
Os 012 22. 0, durch S,, S,, ..., 8, ersetzt, «dann ist die rechte Seite eine obere Schranke 
von |#,4,,, die wir mit X, +, bezeiehnen wollen. Führen wir die gleiche Überlegung für 
» vr +1, r+2,... dureh, dann erkennen wir, daß aus dem Gleichungssystem 


A,d(S, ) S-)+:.:+4,r. (Su 3 5) | 


ur tr — r+ De, +, + N, _, n=r, r+l, 9) 
rekursiv eine Folge oberer Schranken X,+,. Ar+s, --- für |e,+,|, ler+s|, -.. gewonnen 
wird. Für die praktische Berechnung ist eine andere Formel günstiger, die man erhält, wenn 


Man erhält so 


(29) 


‚ ) V 4 g f | | q Ip 
Ant: 2A, X, 1 A, Ay r4, Ay MER A, X, r TC, 


| (24) 
(n=r+lr+2,..)J ne Ze; ;; 


Man beachte, daß diese Gleichung auf Grund ihrer Herleitung nicht für n=r gilt. Um 
\,+, zu berechnen, muß man auf Gl. (23) zurückgehen. Diese lautet im Spezialfall n=r 


, ) v r ’ 4 | r i Ir 
Ars 32 5 2. a 2. HE E04, 5a 


unterscheidet sich also von Gl. (24) durch das Vorzeichen des zweiten Gliedes auf der rechten 
Seite. Damit ist die Fehlerabschätzung grundsätzlich erledigt. 

Der Umweg über die durch Aufsummierung aus (19) entstehende Ungleichung (22) und 
über die darauffolgende Differenzenbildung, die zu (24) führt, ist nicht zu vermeiden, da man 
unmittelbar aus (19) nur eine Rekursionsformel X,+,=2X,+X,„-,+:... (mit dem Plus- 
zeichen auf der rechten Seite, während (24) dort ein Minuszeichen hat!) folgern könnte, die 
wesentlich größere, praktisch nieht brauchbare Schranken liefern würde. 

Setzen wir beispielsweise (für den Fallr= 4) N =-X, =-X\,=A, A, 0, so folgt aus 
Gl. (25) und (17) 

) 


u 2 ’(5) 
elf 


max * 


Alleemein hat also die Fehlerschranke für einen Punkt die Größenordnung Ah’”+°, wenn die 

vorangehenden Punkte fehlerfrei sind. €. Störmer kommt an der in der Fußnote °’) ge- 

nannten Stelle auf anderem Wege unter Vernachlässigung höherer Potenzen von h zu dem 
+.) ‚e) 


etwas günstigeren Zahlenwert 80 


0 an Stelle des obigen 





ech. 


IS) 


19). 


Als 
Ah- 


2). 


'n 
In 


d. 


.). 
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3. Independente Abschätzung des Fehlers. Während die Benutzung der Gl. (24) als 
Rekursionsformel erfordert, daß für die Bereehnung einer Schranke X, für den Fehler |e, 
sämtliche vorangehenden Schranken bis X,_, schon bekannt sein müssen, werden wir jetzt 
eine andere Formel angeben, die \e,) abzuschätzen gestattet unabhängig von der Kenntnis 
ler Schranken für die früheren Fehler. Dazu gehen wir auf die Differenzengleichung (24) 
zurück. Sobald aus Gl. (25) X,;, berechnet ist, können jetzt X,,X,....,X,.;, als die 
segebenen Anfangswerte angesehen werden, die zusammen mit der Gl. (24) eine Lösung X, 
eindeutig bestimmen. Wir werden eine andere einfach gebaute Lösung X, der Differenzen- 
»leiehung (24) angeben (die natürlich zu anderen Anfangswerten gehört), so daß stets X, X, 
(in 1.2,3,...) Ist. Damit ist ebenfalls eine Abschätzung für die Fehler |e,' gegeben, die 
allerdings im allgemeinen ungünstiger sein wird. Die zu benutzende Lösung X, der Differenzen- 
rleichung (24) setzen wir zusammen aus einer partikulären Lösung € z," der zugehörigen 
homogenen Gleichung und einer partikulären Lösung —e,: (J4,+ 1, +...-+ A,) der inhomo- 
venen Gleichung. Dabei ist z, diejenige Wurzel der Gleichung 


27 de 


sr 22’ et A 2" + A,2 +... +4,-,2 +, : .: 2:0. . 26), 


die reell und >1 ist. (Um zu zeigen, daß es stets genau eine solche gibt, dividiere man 


(11. (26) durch z’, setze — Z und berücksichtige, daß A,, A,,..., A, nichtnegativ sind.) Nun 


eilt der leicht zu beweisende Satz: Sind X, und X, zwei Lösungen der Differenzengleichung (24) 
und setzen wir X, — X,=0o,, dann ist 


Ö 0 h) 


i , p) . a 1 


0 Q, 0). Ör+; U, +4; N) u , , 


/ 


hinreichend dafür, daß d, 0, also X,’ X, für alle » ist. Wir wählen daher in 


Pi _ = r IN 
Kr up TE Zu) A A ea 


”. A,+4,+...+4, 


die Konstante Ü derart, daß die Bedingungen (27) erfüllt sind, daß also 
ry z r u 7 r A. r rs r rn ‚  ry r 
BE: Baer re ur Brass ha Ars er =, 


eilt, wobei die X,, X,, ..., X7,+, die bekannten Anfangswerte sind. Diese Wahl ıst stets 
möglich. Man kann beispielsweise wählen 


Cr 


| (29), 


(! X+ 


wobei X eine gemeinsame obere Schranke von X,. X, ..., Nr+, Ist. Bei dieser Wahl von € 
erhalten wir unter Berücksiechtieunge von (18) die Fehlerformel 








N : \ + r u N ce 0020 
— j al Ba xh? 4” Pen 
mit r (80). 
m \' Br 
UV 00r 
Für die A’ hat man folgende Zahlenwerte: Wo 
{ 19 0.008 
2) __ (3) .) ? 
ru, AP—=—, Fu Aw 15° 
+.) .») U 
Zur Berechnung der Wurzel z, ist die Gleichung . 
0,00. 
j ar + 1 2 z' de a’ 1 j . 2000 er | 
W,(z) =. . . > 2 > 22 + % , BD wo I m I 10 78 
[A424 1 Z 
\' l-r [7] i , 
Pa Aro = N Abh. P, 
0= 0 
aufzulösen. Um bequem einen Näherungswert für z, zu finden, der nötigenfalls zu verbessern 


ist, sind in Abb. I die Kurven W,(z) für r -1,2,5,4 für das Intervall z==1,00,...1.25 
eezeichnet. 











‚ . .. e en . Ztschr. f. angew 
230 Schulz, Fehlerabschätzung für das Störmersche Integrationsverfahren Math. und a 





Auf eine gewisse Lücke in der Fehlerabschätzung muß noch hingewiesen werden. Die 
Er 

Berechnung des Quadraturfehlers e, verlangt eine Abschätzung von dgr+i auf der exakten 
Integralkurve. Diese (r-+1)-te Ableitung läßt sich ausdrücken durch f, fx, ft; fxx, -.. (ab- 
hängig von « und f) sowie durch x und = selbst. Höhere Ableitungen &, x" ... treten wegen 
“= f nicht auf. Für »(f) kann man nach Durchführung der näherungsweisen Integration 
(etwa aus der Zeichnung) plausible Schranken annehmen, sodann die Fehlerabschätzung durch- 
führen und nachsehen, ob die Schranken weit genug waren. Für & (t) kann man ebenfalls 
(aus der Zeichnung oder aus dem Differenzenschema) leicht plausible Schranken annehmen; 
aber man kann nachträglich nicht feststellen, ob die Schranken weit genug waren, da das 
Störmersche Verfahren keine Werte von .r liefert. 


Praktisch ist diese Lücke ohne Belang, da man bei genügend enger Intervalleinteilung 
dem Differenzenschema sehr sicher Schranken für x und & entnehmen kann, die man ziemlich 
roh vergrößern mag, ohne daß dies auf die Abschätzung von |f”+"”,., wesentlichen Einfluß 
hat. Eine praktisch noch bequemere Abschätzung ist in dem in 5 durchgeführten Beispiel 
angewendet worden. Bei der Anwendung des Adamsschen Verfahrens auf die Differential- 
gleichung = — f(x, f) tritt die geschilderte Schwierigkeit in der Fehlerabschätzung nicht auf. 


4. Systeme von Differentialgleichungen. Die in 2 und 3 gegebenen Fehlerabschätzungen 
lassen sich ohne besondere Schwierigkeiten auf Systeme von Differentialgleichungen ver- 
allgemeinern. Im folgenden werden die Formeln für das System (2) von zwei Differential- 
gleichungen angegeben; wie sie für mehr als zwei Gleichungen lauten, ist daraus ohne weiteres 
zu erkennen. Die Beweisführung wird nur an den Stellen ausführlicher gegeben, wo grund- 
sätzlich neue Überlegungen gegenüber dem Fall einer einzelnen Gleichung hinzukommen. 


Gereehnet wird mit den Formeln 


5 Ki ri - h’ (d, fn+ 3: Eu + +2, 4" fn) 


E u (32), 
I? An rs h? (8, AntB: Igun+t- . + Pr 1" Yu) | 


wobei zur Abkürzung f(#,,9,,H) = f, und 9(#,, 4», = 9, gesetzt worden ist, nach dem am 
Ende von 1 gerebenen Schema. 


Wir setzen die Gültigkeit der Lipschitzbedingungen 


fu d- Far dere try) 
a ER x | 5 A 
74H. —- ar, dir +, ly— el] 


für die rechten Seiten der Gl. (2) voraus. 2x. #y, /x, /, Sind vier gegebene nichtnegative 
Konstanten. Abzuschätzen sind die Fehler 


„us Er); weh Wh): . » x 3 2» 5 ır - Ta. 
An die Stelle der Ungleichung (14) bzw. (19) treten die Ungleichungen 


1 
j? € vun + . \’ (Axo Ex, n #) + Ayo ey, n 7) ) -r- C, mit C, — ng + ee + Fe 
. pe ” = > - 
/ ud 
ö . (BP). 
1 eyınHı Y (Bxo exX,n— Oo r Bo Ey, n— 0) 7 d, mit d, d, + d, +d,” 


(n=r, r+1l,r+2,...) 


Vz \) 
\ 





Dabeı ıst zur Abkürzung gesetzt 


„ 


Ax h? laro‘ Axo, Ay h° Adro 4, DO» Fi h? Adro Bxo, /y h? laro!| Bo ö (36). 
Ferner sind e,', e,", ce, , d,, d,', d’" so definiert wie er, €, , cr ın 2 und 3, nur daß an 
Stelle von f(x, #) jetzt fe, y,M) und 9(#,9, N tritt; und zwar gilt auch hier 


h" ui “r+9 I l h" er r+»D I 37) 
r (r +1)! / max r’ (ir (r +1)! 4 ® ’ ° ‘ . . . . (4 . 
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(Gegeben sind obere Schranken X,, X,.:::: Ay, Yar Yır-:-, Y. für die Fehler der 
ersten r + 1 Punkte, so daß also 
wi 4 ü we . d | 4 
ER ä., Is; uses fü, ER | 
| (38). 
Eyo : Fu; Ey; > Y BAG ey og ; } r | 
Die gesuchte Folge der Fehlerschranken X, ,, AXr+3s+: Yr+ı1+ Yrtas... berechnet sich aus 
r r 
r fr r ' r v r v r 
AÄr+: 2 Er A, 1 +2 Axo “ir 0 +> Ayo Y; 07 C, 
Q 0 v 4 
(39). 
r F 
Yr+ı . 74 Y,+ Y, 1 Ze: Bxo Är QO +3 Byo Yr 0 4 d, 
0=V 0 — U 
Ä r 
r fa r r y r y r u 
rw 2X, Bu ‚+2 Axo u o+ > Ayo Ya O0 T C; 
0 0 [% 0 
! (40). 


. r 
Y, + 1 2 5. ge 1 u \ . Bxo An _ LU r \ y RB, (2) Ei #) H d 
— ” —— - ” 


0o=—=V 0 IH 
S N 





(n=r+1l,r+2,...) 


Um eine Abschätzung von |&x,„|; £&y„| unabhängig von der Kenntnis der Schranken 
der vorangehenden Fehler zu erhalten, berechnen wir zunächst aus Gl. (39) X, ,. Yr+, und 
sehen X,, X3,..:, Artı, Yır Yas- +» Yr+ı als die gegebenen Anfangswerte an. Diese be- 
stimmen mit den beiden Differenzengleichungen (40) eindeutig eine Lösung X,„. Wir werden 
auch hier wieder noch eine andere Lösung X,', Y„' der Differenzengleichung angeben, so daß 
stets X, = Xu, Ya Ya n=1,2,3,...) ist. Ein Lösungssystem für das zu (40) gehörige 
homogene Gleichungssystem ist Cz,”", Dz,", wobei z, eine Wurzel der Gleichung 

#.h?” —W,(z) %, h? 
= () ee de 5 
A. R® A,h? — W,(2) 


ist. (Die Funktion W,(z) wurde in Gl. (31) definiert.) Setzt man nämlich Cz,” und Dz,” 
in das zu (40) gehörige homogene Gleichungssystem ein, so erhält man für € und D das 
homogene lineare Gleichungssystem 


r r 
(zr+ 2 ER „Tr op ® \’ Axo zı 0\, C+ \’ Ayo zr 0.D —{) 
( 0 d u 


(42). 





r r 

N’ Bxo272.C+ a’t! 22” +27") + V’ Byoz'-2\-D=0O 
year a u 

Damit dieses System eine nichttriviale Lösung €, D zuläßt, muß notwendig seine Determinante 
verschwinden. Das gibt gerade Gl. (41) zur Bereehnung von z,. Statt ihrer betrachten wir 
zur genaueren Diskussion nacheinander die Gleichung für 1 


x. h? — A x, h? 
i 0 Be ne . (445) 
Ah? i,h?— A 
und die Gleichung für z 
Wie... ir anne 


Hat man nicht ein System von zwei, sondern von m Differentialgleichungen, dann tritt in der 
„Säkulargleichung“* (43) eine analog gebaute m-reihige Determinante auf. Zu jedem die Gl. (41) 
befriedigenden z2-Wert können dann aus Gl. (42) die zugehörigen Werte von € und D bis auf 
einen gemeinsamen Faktor bestimmt werden. Für die majorante Lösung X,, Y, Kommt 

ebenso wie im Fall einer einzelnen Differentialgleichung — nur eine solche Wurzel z, ın 
Frage, die reell und >1 ist. Eine solche existiert aber immer. 
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Nach einem Satz von G. Frobenius hat nämlich eine Matrix mit nichtnegativen Ele- 
menten, wie es die Matrix der mit A? multiplizierten Lipsehitzkonstanten ist, mindestens L; 
eine charakteristische Wurzel, die reell und positiv ist®). Im Falle » =2 läßt sich der Be- 
weis dafür auch unmittelbar führen. Setzen wir nämlich e) 


“ . er V 
Kx . Ay nz 2 p . Ax Ay ur Ay Ax | . 


so folgt für die beiden Wurzein A,, A, von Gl. (43) 


A,.,=W(» typ’ ). 
Da »?> I, so ist 1, sicher reell und positiv. 


Die Gleichung W,(z)= A hat ferner für 1>0 immer genau eine reelle Wurzel, die 
I ıst. Analog zu Gl. (25) setzen wir also 


Aycı Ay d, zr AxCr Rx d, 


An Ü2," pP D A + h? A» | ° . . . . (49). 


ii. h’ATA 
Die (von » unabhängigen) Konstanten auf den rechten Seiten der Gl. (45) sind partikuläre 
Lösungen des inhomogenen Gleichungssystems (40). X, und Y,„ befriedigen also dieses \ 
System (40). Der obige Ansatz für X, und Y, ist ersichtlich nur statthaft, wenn 4=F0 ist, 
was durch Vergrößerung einer Lipschitzkonstanten stets erreicht werden kann. 4A” ist 
in (30) definiert. Daß es nur einen brauchbaren Wert von z, gibt, wird im folgenden gezeigt. 


o>5 


Die Lösung (45) von (40) ist nur dann als majorante Lösung brauchbar, wenn C und D 
einerlei Vorzeichen haben, so daß C und D als positiv gewählt werden können. Wir be- 
schränken uns hier der Einfachheit halber auf den Fall, daß die Matrix 8 der Lipschitz- 
konstanten keine verschwindenden, sondern nur positive Elemente hat. In diesem Fall 
lautet der bereits erwähnte Satz von Frobenius: 


Sind die Elemente einer Matrix alle reell und positiv, so hat ihre charakteristische 
Gleichung eine Wurzel, die reell, positiv, einfach und absolut größer ist als jede andere 
Wurzel (sog. „Maximalwurzel* der Matrix). 


"robenıus zeigte ferner: Ist 4” die Maximalwurzel einer Matrix WU mit positiven, 
reellen Elementen, dann hat das homogene lineare Gleichungssystem mit der Matrix M — 4*€ 
nur eine (bis auf einen gemeinsamen Faktor bestimmte) Lösung, deren Elemente alle als 
positiv angenommen werden können. Das homogene lineare Gleichungssystem mit der Matrix 
M - IE kann ferner nur dann eine Lösung haben, deren Elemente alle nichtnegativ und 
nicht alle Null sind, wenn I= .1* ıst; dann sind aber alle Elemente der Lösung positiv '”). 


Bezeichne jetzt 1” die Maximalwurzel für die Matrix "’St und z* den z-Wert, der sich 
aus W,(2)= .1* eindeutig ergibt. Dann ist, wie bereits gezeigt, z* reell und > 1, also als z, 
in (1. (45) brauchbar. € und D (oder im Fall von m Gleichungen, m Konstanten CC, C®,..., 
CP?) bestimmen sich bis auf einen gemeinsamen Faktor aus dem homogenen Gleichungs- 
system (42) mit 2 z*, das die Matrix 


"R—-W,e@!)E=MRKR—-AE .. . . 2 20202 2 0.46) 


besitzt, und können nach dem vorstehenden Satze sämtlich positiv gewählt werden. Eine 
andere Wahl von z, ist nicht möglich; selbst wenn Gl. (41) noch eine andere Wurzel hätte, 
die ebenfalls reell und > 1 ist, so wäre diese für den Ansatz (45) unbrauchbar, da dann € 
und D nach dem oben angeführten Satze sicher nicht einerlei Vorzeichen haben. z, ist also 
immer eindeutig bestimmt, wie groß auch m, die Anzahl der Differentialgleichungen des 
Systems sei. Die praktische Berechnung von JA aus Gl. (453) ist auch insofern wesentlich 
erleichtert, als man nur die Wurzel vom größten Absolutbetrag zu berechnen hat. 


Der in € und D noch unbestimmte gemeinsame Faktor ist erstens so zu wählen, daß 
C und D beide positiv sind und so groß ausfallen, daß 


ai. War u BE ea Ki cn — 5 
(47) 
Var. Wer Bein Frei Fra ns 


(vel. (27)) gilt, was stets möglich ist. 


», (4. Frobenius: Sitz.-Ber. d. Preuß. Akad. d. Wiss., Math.-phys. Kl. 1908, S. 471 bis 476: 1909, S. 514 bis 518; 
1912, S. 456 bis 477. 


10) Siehe die dritte der unter ®) angeführten Arbeiten”S. 4654 bis 465. GC ist die Einheitsmatrix. 
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Kle- Ist 4=0, dann sind Polynome in » eine partikuläre Lösung des inhomogenen Systems (40). 
EnSs Liegt J nahe bei Null, dann werden die zweiten Summanden auf den reehten Seiten von (45) 
Be- sehr groß und überwiegen für die in Frage kommenden Werte von n stark gegenüber den 
ersten. Es empfiehlt sich daher, in einem solchen Fall eine Lipschitzkonstante so zu 

j vergrößern, daß A=0 wird. 


5. Beispiel. Als Beispiel betrachten wir das System 


t=uae"’—e*’+e’’cos’y=f(x,Yy) | 


| E re 
die _ SNY ı „x. am 
y= cos y re sı1n cos Yy Ar, My) 
mt 
- l 4 ” 
15). =>——. 7.=-,=09. a=0070600 
16 y,* ’ı / 

Are 
ses , unter den Anfangsbedingungen 
Ist, 

ist 0621851, 20 | 
of t ” - für f 0. f i 5 N ; . \ . HAN. 
er yo) 9 0,174496 | 

D 6 | 

y n 6 

be- Die Integralkurven .z(M), y(t) sind von C. Störmer für f=0.., 8 mit ho, berechnet 
Z- . —— 

all worden. Abb. 2 und 3 geben den Verlauf wieder''). 





he 
re 
"n, N ji 
& y [ 
3 h 
ıls ( / 05 
IX 5 f “ 
nd N) 7 2 24 02 7 04 21 06 07x 
65 99 # 
i Pe FERE 5 
et FE hr 
*h 0 — , 
= -0,3 . Ct 55 
'S- Abb. 2 Abb. >. 
Für die Fehlerabschätzung ist es im allgemeinen erforderlich, daß die eigentliche Rech- 
) nung bereits völlig durchgeführt sei. Obwohl Störmer mit dem Näherungsgrad r = 4 rechnete, 
m wird der Fehlerabschätzung r =3 zugrunde gelegt, da die vierten Differenzen schon sehr stark 
z mit Abrundungsfehlern behaftet sind und die Zahlenwerte kaum beeinflussen. Wir beschränken 
f' uns auf das Intervall 0=t=3 der unabhängigen Veränderlichen, also auf 24 Sehritte. Aus 
m den Zahlen der durchgeführten Rechnung kann man vermuten, daß; 
Ss .n_ x 3 
| VD <rN<OGH5, 0: yıh) "4% We IV I AT) re 
h 
eine rohe Eingrenzung von x und y ist (vgl. Abb. 2). Die folgende feinere Fehlerabsehätzung 
zeigt, daß diese Schranken nicht zu eng sind. Für x... %,. 4x, 4, nehmen wir obere Schranken 
3 5 & x N X R/] 
von 
| 5 » P [3 1 | [3 | » 
x (A h) Y) max» fu («u ’ Yy) max» Ix (. ’ Y) max) Ju (4 ) Y MAX 
) für den Bereich (50). Es ist hier übrigens f,, = 9x, also 2, 4x. Als Zahlenwerte ergeben sıch 
2. = 0AWI, 2,0148, Ar =0,1418, A, 0,9509. 
> 11) Die Kurve (x) ist auch abgebildet bei C. Störmer: Ztsehr. f. Astrophys., Bd. I, 1930, S. 237 bis 274, und 


a von rechts. Die 


zwar ist es in Abb. 9 (S. 248) die dritte Kurve von rechts und in Abb. 13 (S. 255) die zweite Kurve 
Bezeichnungen der Konstanten a,y,yı entsprechen denen bei Störmer., 








a. ge 
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d’f d’g r az 
Um m und ip abzuschätzen, gehen wir, da f und 9 durch verhältnismäßig 
IMAX MAX 


komplizierte Ausdrücke gegeben sind, hier nieht so vor, wie am Schluß von 3 angedeutet 
wurde, Wir beachten vielmehr, daß 


fl Wi ME "y_ A* 7 j? 7 
der am h ’ dt At h* 


ist (£=h’f, 7 = h’g) und wählen eine obere Schranke aller 4’& und 4*), die im Differenzen- 
schema der Rechnung stehen. Dort findet man, daß 


W- u. 6 | A Pr ® 6 
u max > 6-10 , j N Iimax OÖ 6-10 ; 
Daher wollen wir annehmen 


d’f __ d’« 
/ <8.6-10 = 1513, J 


dt max d " 


<1513. 
max , 
Da wir nur den reinen „Formelfehler* untersuchen und deshalb von allen Abrundungsfehlern 
absehen wollen, setzen wir die Anfangswerte &,, X, X, #5 Yos Yıs Ya, Y, als fehlerlos vor- 
aus, setzen also 


Sl = Lu = 1, =] 
Aus Gl. (39) folgt dann X, =e,, Y 


-h=-h-=0. 
‚= d,. Und zwar ergibt sich aus Gl. (37) 


= d,=0INT 10°, 


Die schrittweise Berechnung der X,, Y, nach den Gl. (40) ergibt schließlich folgende Zahlen- 
werte: 








v X, + 10° Y, + 10° v X, + 106 Y,.10° 
D 2.384 2.391 15 75.83 86,90 
6 4.799 4.835 16 91.79 107.6 
7 8.073 8.180 17 110.2 132,2 
S 12.26 25 | 18 131.5 161.5 
9 17.43 18.03 19 156.0 196,4 

10 23.68 24,81 0 184,2 237,8 

11 31.12 33.08 21 216,6 2871 

12 39,87 13.10 22 254.0 345.7 

13 0.11 9,15 2 297,2 415.3 

14 62.02 69,61 24 346,8 497.9 


Um die Fehlerabschätzung auch nach der Methode der majoranten Partikularlösung 
durchzuführen, berechnen wir aus Gl. (43) und (44) die Wurzeln A und z und finden 


A, 0,015 484 :, — 1,182 


2 1.120. 


A, = 0,007 031 


Wir wissen bereits, daß nur die Maximalwurzel z, brauchbar ist und setzen daher z, — 1,182. 
In der Tat liefert Gl. (42) für z, das Verhältnis U: D--+0,2831, für 2, aber ein negatives 
Verhältnis. Die Determinante der Lipschitzkonstanten ist 4 0,4459. Aus den Gl. (45) wird 

K. = U2," 19,97 - 10° ® mit 2,=1,182 


»1). 
Y'—=Dz"--19,79:10°% | und 0: D= 0,2831 | m 


Die Bedingungen (47) sind erfüllt, wenn man 
©4597 :10°%, D-16244 -10°*® 


setzt. Für » = 24 erhält man X,, > 2500 -10°°%, also einen mehr als 7 mal so großen Wert 
wie bei der rekursiven Fehlerabschätzung, so daß man hier auf die Benutzung der Formel (45) 
besser verzichten wird. Bei kleineren Werten von z, können aber beide Verfahren durchaus 
konkurrieren. 424 
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Über die kleinen Schwingungen eines Pendels 
mit oszillierendem Aufhängepunkt. 


Von Artur Erdelyi in Brünn. 

Einleitung. Der Stabilisierung instabiler Schwingungsvorgänge durch Überlagerung 
kleiner, aber sehr rascher Schwingungen kommt in der Technik eine stets wachsende Be- 
deutung zu. Es ist bereits darauf hingewiesen worden, daß das einfachste mechanische Modell 
eines soleherweise stabilisierten Systems als ein Pendel gedacht werden kann, dessen Auf- 
hängepunkt in kleine, aber rasche Schwingungen versetzt wird. So dürfte die Behandlung 
dieses Problems mehr Interesse verdienen, als es auf den ersten Bliek der Fall zu sein scheint. 

Die umfangreichste Behandlung dieses Problems rührt von Paul Hirsch!) her, der 
sich in der unten zitierten Arbeit mit den Gleichgewiehtslagen solcher Pendel beschäftigt 
hat. Leider geht aus der ganzen Art seines Näherungsverfahrens der Gültigkeitsbereich seiner 
Resultate nicht klar hervor. 

Hirsch behandelt das Problem statisch, indem er den „Effekt“ der Oszillation, d. h. 
len Mittelwert der Beschleunigungen des Pendels während einer Schwingungsdauer des Auf- 
hängepunktes berechnet. Die Gleichgewichtslagen sind dadurch bestimmt, daß für sie die 
mittlere Beschleunigung die Erdbeschleunigung gerade aufhebt. Ohne Kenntnis der wirklichen 
Bewegung des Pendels aber kann man diese Berechnung nur unter vielen und zum Teil un- 
kontrollierbaren Vernachlässigungen durchführen. Die Existenz kleiner Schwingungen ist die 
erundlegende Voraussetzung hierfür; es kann aber nicht gezeigt werden, daß solche Schwin- 
ungen wirklich möglich sind. Auch die Voraussetzung, daß der Schwerpunkt des Pendels 
bei der Oszillation seine Mittellage unter einer solchen Richtung passiere, daß in diesem 
Augenblick die vom Antriebsmechanismus ausgeübte Kraft verschwinde, ist nicht ohne weiteres 
einleuchtend. Ebensowenig erscheint es uns zwingend, daß die Bewegung des Schwerpunktes 
in erster Näherung geradlinig erfolgt. Dazu kommt noch, daß die Gleiehgewichtslage bei 
dieser statischen Behandlung gar nicht definiert werden kann, wie Ph. Bock gezeigt hat’). 
Freilich geht aus der von Ph. Bock gegebenen Darstellung des Problems und der gegen- 
wärtigen Arbeit hervor, daß die von Hirsch erzielten Resultate für sehr rasche und sehr 
kleine Schwingungen des Aufhängepunktes gültig sind. Trotz dieser Unsicherheit ist es aller- 
dings ein Vorteil der von Hirsch angewendeten Methode, daß man mit ihrer Hilfe ziemlich 
allgemeine Schwingungsformen des Aufhängepunktes untersuchen kann. 

Professor Weyrich hat im mathematischen Seminar der Deutschen Technischen Hoch- 
schule in Brünn bei der Behandlung der Mathieuschen Funktionen und verwandter 'Tran- 
szendenten darauf aufmerksam gemacht, daß im Falle einer linearen, harmonischen Oszillation 
des Aufhängepunktes in der Schwererichtung die Differentialgleichung der Pendelbewegung 
für kleine Schwingungen eine inhomogene Mathieusche Differentialgleichung ist. Im An- 
schluß daran hat Ph. Bock (Dissertation, Prag 1932) nachgewiesen, daß die Hirschschen 
Resultate das erste Glied einer asymptotischen Reihe sind, so daß man das Gebiet ihrer 
Gültigkeit und auch die Grenzen ihrer Leistungsfähigkeit klarer beurteilen konnte. 

In der vorliegenden Arbeit wird, ausgehend von der Differentialgleichung der Bewegung, 
die Theorie in zwei Richtungen weitergeführt. 

Es gelingt unter Annahme einer beliebigen (periodischen) Bewegung des Aufhängepunktes 
in der Schwererichtung die Existenz kleiner Schwingungen um eine Gleichgewichtslage sicher: 
zustellen. Für eine spezielle Bewegung des Aufhängepunktes (Annäherung der Sinuskurve 
im Weg-Zeit-Diagramm durch sich in den Nullstellen der Sinuskurve berührende Parabel- 
bogen)*), die die harmonische Bewegung hinreichend gut annähert, kann man die Differential- 
gleichung mittels elementarer Funktionen integrieren und die Gleichgewichtslagen numerisch 
bestimmen. Es zeigt sich, daß es eine kritische Frequenz der Oszillationen des Aufhänge- 
punktes gibt, unterhalb welcher sich die Stabilisierungserscheinungen überhaupt nicht zeigen. 
Für alle Frequenzen oberhalb der kritischen gibt es zwei in bezug auf die Schwereriehtung 
symmetrisch gelegene Gleichgewichtslagen. Die Gleichgewichtslagen werden definiert als der 
Mittelwert der Pendelstellungen. Ist © der Winkel, den die Pendelstange’) mit der nega- 
tiven Schwererichtung einschließt, so ist die Gleichgewichtslage gegeben durch: 

t 


a = lım \ dt. 
Is 
0 


I), Paul Hirsch: Das Pendel mit oszillierendem Aufhängepunkt. Zeitschrift für angewandte Mathematik und 
Mechanik, Band 10 (1930), S. 41—52. 
2) Dissertation. Prag 1932. 
3) Die Anregung hierzu verdanke ich Herrn Professor Weyrieh. 
3) Als Pendelstauge bezeichnen wir die_Verbindungsinie des Pendelschwerpunktes mit dem Aufhängepunkt. 
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Durch «diese Festsetzung wird die Gleiechgewichtslage erst definiert. Aus der aufzustellenden 
Differentialgleichung folgt dann unter Heranziehung der Existenzsätze die Möglichkeit kleiner 
Schwingungen um solche Gleichgewichtslagen. Es ist beabsichtigt, in einer zweiten Mitteilung 
die Gleichgewichtslagen für allgemeinere Schwingungen des Aufhängepunktes zu ermitteln 
und auch die Schwingungen zu behandeln, die nicht als kleine Schwingungen um eine Gleich- 
rewiehtslage beschrieben werden können. 


Als ganz neuer Gesichtspunkt wurde in der vorliegenden Arbeit die Dämpfung des Systems 
in Rechnung gezogen, die bei den vorgenannten Autoren vernachlässigt wurde. (Beim mecha- 
nischen Modell des Pendels macht sich diese Dämpfung als Luftreibung und Reibung des 
Pendelgelenks bemerkbar.) Es zeigt sich zwar, daß die Dämpfung, solange sıe klein ist, die 
numerischen Resultate nicht wesentlich beeinflußt, wie das nicht anders zu erwarten ist; ab- 
vesehen aber davon, daß sie den Schwingungsvorgang wesentlich beeinflußt, wenn sie nur 
einigermaßen groß ist, ist ihre Einführung deshalb von Wichtigkeit, weil sie die ganze Theorie 
von einer gewissen Willkürlichkeit und Unbestimmtheit befreit, die ihr sonst anhaften würde, 
da durch sie einwandfrei entschieden werden kann, welche partikulären Integrale der Diffe- 
rentialgleichung für die Herstellang der Lösung heranzuziehen sind, was bei Vernachlässigung 
der Reibung von vornherein nicht ohne weiteres möglich ist. 


I. Allgemeiner Teil. 


Problemstellung. Es sollen die kleinen Schwingungen eines (physischen) Pendels um eine 
Mittellage behandelt werden, wenn der Aufhängepunkt desselben eine vertikale Schwingung 


mit der Amplitude 4 und der Schwingungszahl » [Schwingungszeit 


‚4 h 
ER: e ausführt. Die Entfernung des Pendelschwerpunktes vom Auf- 


hängepunkt sei s, sein Trägheitsradius in bezug auf den Aufhänge- 
punkt k und die reduzierte Pendellänge I (siehe Abb. 1). Hat man 
ein mathematisches Pendel, so bedeutet !=s=%k die wirkliche Länge 
des Pendels’). 


Die Differentialgleichung der Bewegung. Wir stellen die La- 
erangesche Gleichung für die Bewegung des Pendels auf. », sei 
1» ass. die Elongation des Aufhängepunktes, © der von der Pendelstange 


und der negativen Schwererichtung eingeschlossene Winkel (siehe 





e bb . . 1} ® rmm\ . . . 
.' Abb. N). M die Masse des Pendels,. kM sein Träzrheitsmoment in 
‚ ö ze! ah b Id) . do 2 
bezuz auf die Schwingungsachse. Wie üblich, ıst Ar „J, ie :9 usw. gesetzt. 


Das Pendel führt eine Translation in der negativen Schwererichtung mit der Geschwin- 
digkeit ;; und eine Rotation um den Aufhängepunkt mit der Winkelgeschwindigkeit © aus. 
Seine kinetische Energie ist ®): 


LM? +2 MO ssin 9+Mk®-OY— 5 Mi? +27s9sin9+R AR. 


— 
— 


Die potentielle Energie des Pendels ıst: 
U=Mgıy+scos% . 
Wenn man L und U in die Lagrangesche Differentialgleichung 


o(Lh—-U) d (“ 4 0 
ie] d\8 0 


einsetzt. bekommt man: 


ko) s(g y)-sın O zul), 


5) Um den Vergleich mit den Hirschschen Resultaten zu ermöglichen, sei hier die Zusammenstellung der Be 


zwwiehnungen gegeben: 
In der vorliegenden Arbeit: u. 2 l, 4, n. 


In der Arbeit von Hirsch: BEE, ar u 
6) Das erste Glied in der Klammer rührt von der Translation, das dritte von der Rotation her, das zweite Glied 


ist bedingt dureh die Coriolisbeschleumgung. 
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k’ 
Bedenkt man, daß die reduzierte Pendellänge 1 = - ist, so hat man: 


/ 


ro, 
EEE A NE 


.. { 
oJ — ' 
Zieht man noch die Reibung in Betracht’), so wird die Gleichung lauten: 


Z 


+22nG BE ee ee 


wo z eine (dimensionslose) positive Konstante ist. 


Wir nehmen zunächst heuristisch an, daß kleine Schwinzeungen um eine von der Zeit 
unabhängige „Gleichgewichtslage* (im Sinne der in der Einleitung gegebenen Definition) 
möglich sind und führen diesen Ansatz analytisch durch, um dann zu zeieen, daß die zu- 
nächst hypothetisch angenommenen Schwingungen wirklich existieren. Wir zerlegen also 
den Winkel 

=a+ty, 


Wo 
! 
I \[ 
«a lim ‚\oat, 
I>ı% 
u 
also 
! 
I ( ” 
im \oat ee 
u 
und 
PR. 
Y .) 


dann ist cosy=1, sing =gy, sin(a+g)=sina+ cosa-gy. Gl. (2) geht dann über in 


| ,. 9 9-Nn . 
q +2ıng - U cosa:g — -sIina . . TE . } u 7 (4), 


eine inhomogene Gleichung, die aber dafür linear ist. Von den Koeffizienten dieser Differential- 
eleichung sind alle bis auf 5) konstant, dieses aber ist periodisch mit der Zeit veränderliech. 


Wir lösen zunächst die homogene Gleichung 


m .). 7 J Be . @ \ ’ v0 ” 
ln N Cosa, a u ee ee 


und führen sie durch die Substitution 
TER AD ee er ER N NT 
in die selbstadjungierte Form 


5 9-1 al z 
ih) wurde ran 


über, eine Differentialgleichung vom Hillschen Typus’°). 


Ist z(f) die allgemeine Lösung von (7), so ist y(f) = err"t.y(h die entsprechende 
Lösung von (5). Da (4) das partikuläre Integral  tga besitzt”), so ist ihre allgemeine Lösung 


ort . 2 ee aa ana er 


7), Wir setzen die Reibung der Winkelgeschwindigkeit proportional. 

8) Bezüglich der Theorie der Hillschen Differentialgleichungen sei verwiesen anf die Monographie von M. J. 
O0. Strutt: Lamesehe, Mathieusche und verwandte Funktionen in Physik und Technik. I. Band, Heft 3 der Ergeb 
nisse der Mathematik und ihrer Grenzgebiete. Berlin 1932, Springer. Da sind neben einer kurzen Zusammenstellung 
der Ergebnisse auch reiehe Literaturangaben über die Theorie und Anwendungen dieser Differentialgleichungen zu 
linden. 

9) Dieses partikuläre Integral hat keine unmittelbare physikalische Bedeutung. Man kann nicht etwa diese 
Zerlegung auffassen als eine in erzwungene und freie Schwingungen, wie diese Begriffe hier überhaupt keine Be- 
deutung haben. 
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Lösung der Differentialgleichung. Aus der allgemeinen Theorie der Hillschen Differential- 
eleichung 
d’u 


L$D(» ; — I) 
je’ Pi2rnt)-uw i 


wo ® eine gerade periodische Funktion ist (Pix +27) P(x)), benützen wir folgende Sätze 
(siehe z. B. bei Strutt a.a.0. Ss. 131F.): 


I. Das Theorem von Floquet'): Die Hillsche Differentialgleichung besitzt zwei (im 
Falle «==0 linear unabhängige) Lösungen von der Form: 


u,()= e""t Fi2nnt) u,()=e-HriF(—2nnt), 


| en ' 
wobei F(2rznt) eine modulo 7 - periodische Funktion von # ist. 
2. Zur Berechnung von « hat man eine Gleichung von der Form 
Bee . ee ra er 


wo f eine reelle Funktion ihres Argumentes ist. (Im allgemeinen hängt die rechte Seite von 
allen verfügbaren Parametern der Differentialgleichung ab, die in ® enthalten sind; in unserem 
Falle ist nur cos«a verfügbar, alle anderen sind durch die Daten des schwingenden Systems 
zereben.) Zunächst ist «= Ar&oj|[fteos a)] unendlich vieldeutig. Da indessen die Wurzeln 
dieser Gleiehung sich nur um ganzzahlige Vielfache von 2 ri unterscheiden, so daß also der 
diesen Summanden entsprechende Teil des exponentiellen Faktors in den Lösungen «,(f) und 
u,(#) periodisch modulo 7 in £ ist und daher zu der periodischen Funktion F(27nt) geschlagen 
werden kann, genügt es, aus den unendlich vielen Wurzeln der Gl. (9) einen Wert als Hauptwert 
herauszugreifen. Als solchen wählen wir denjenigen Wert, für den 


Reu=0 und 0= mu<2a"). 


3. Die Sätze von 0, Haupt'’) besagen in unserem Falle, daß für jedes vorgegebene 
reelle « die Gl. (9) abzählbar unendlich viele reelle Wurzeln 


cos a, =0,1,2,...) 


hat, und daß die zu cosa,,„_, und cosa,, gehörenden Lösungen der Differentialgleichung 
eenau » Nullstellen in dem Intervall (0, T) haben. Die Gleichung 


Go] u —=f(ecos a) 
hat ebenfalls abzählbar unendlich viele Wurzeln 
COS Ay (=1,82.:..); 


und es haben wieder die zu cosa,,_, und cos a,, gehörenden Lösungen der Differential- 
eleichung » Nullstellen in dem Intervall (0, T). 


Unsere Aufgabe ist es, solehe Lösungen z(f) der Differentialgleichung (7) zu finden, 
daß (S) eine für alle Zeiten kleine Funktion mit dem Mittelwert O ergibt. Wir müssen dabeı 
zwei Fälle unterscheiden, je nachdem tga=0 oder tga=-0 ist. Wir führen hier den Existenz- 
beweis zunächst für den allgemeinen Fall tga=+0 und werden dann die Diskussion des Falles 
tz a—=0 nachtragen. Für tg a=+-0 zeigen wir zunächst, daß obige Aufgabe identisch ist mit 
der Bestimmung solcher von Nullstellen freier Lösungen der Differentialgleichung (7), für 
die u=x Ist. 


Zunächst muß die Lösung z(f) frei von Nullstellen sein. Hätte nämlich y(f) eine Null- 
stelle #,, so wäre an dieser Stelle g (#,): tga, also gewiß nicht sehr klein. (Den Fall 
tz a0 haben wir ausgeschlossen.) 


Jetzt zeigen wir, daß wir nur Lösungen brauchen können, für die u=x ist. Aus dem 
Floquetschen 'Theorem folgt, daß y sich linear zusammensetzt aus 


z(d=etrtFeant) und „d=ertrtiF(—2ant) . . .. 0.0. (AO. 


Wir müssen durch die Wahl von cosa erreichen können, daß (8) die oben angegebenen 
Kigenschaften besitzt. 


10), Floquet: Sur les equations differentielles lineaires 5 eoefficients periodiques, Paris, These, 1879, Ann. 
Keole norm. 12 (1883), S. 47-8, 

11) Inter Ne und Am seien der reelle und imaginäre Teil einer komplexen Größe verstanden. 

12) 0, Haupt: Über lineare homogene Differentialgleichungen zweiter Ordnung mit periodischen Koeffizienten, 
Math. Annalen Bd. 79 (1919), S. 278 bis 28. 
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tıal- Aus (6) hat man dann 
w()=e-*tWntFlö2nnt) und y()=el-*-WntF(—2ant) . . . (MM. 
Wählt man cosa so, daß Ne u>x, so ist y,(f) instabil, d. h. ihr absoluter Betrag wächst 
ätze mit zunehmendem ?t über alle Grenzen, und y,(f) gedämpft. (In diesem Falle ist nämlich 


Nel- #2 +tm)>O und Rel— x — u)<0.) Dieser Fall, wo nach ($) auch 9(f) über alle Grenzen 
wachsen oder gegen — tga gehen muß, kann gewiß keine kleinen Schwingungen liefern. 

(im Wählt man aber cosa so, daß NReu<z, so ist Rel- x +u)<O und erst recht 
el 2 u)<O. (11) zeigt dann, daß beide Lösungen der homogenen Gleichung abklingen 
und daß die Lösung von (4) sich nach hinreichend langer Zeit auf —tga reduziert. Da wir 
den Fall tg«-- 0 einer späteren Diskussion vorbehalten haben, so können wir wieder sagen, 
daß auch Re u < x keine kleinen Schwingungen mit dem Mittelwert O liefert, wie sie (3) verlangt. 

Erst wenn cos«a so gewählt wird, daß Re u —= x, hat man eine stabile, nieht abklingende 
Lösung, von der wir sogleich zeigen werden, daß sie sogar periodisch ist, und eine zweite, 


(9) die allerdings abklingt, mit der Dämpfung 2x. 

von Da Eoj u —=f(cosa) reell ist, wie wir schon früher betont haben, muß mu 07 sein, 
em « selbst muß dann die Form haben 

MS n=x+tonrt Rh © 5% Re, u T ; ; ; 3 
eln 

der Nach der Definition des Hauptteiles von « können wir uns auf die Werte 00 und o = | 
ınd beschränken. Wir werden sogleich zeigen, daß auch o = 1 auszuschließen ist. 

Au Aus den oben zitierten Sätzen von Haupt folgt nämlich, daß zur kleinsten Wurzel der 
ert (Gl. (9) eine nullstellenfreie Lösung z, (f) gehört. Für diese Lösung ist * 


v, d=F(2nnt) 
also, wie wir behaupteten, eine periodische Funktion mit der Periode T. Zu allen anderen 
un Wurzeln der Gl. (9) gehören Lösungen der Differentialgleichung, die Nullstellen haben, also 
für unsere Zwecke nicht brauchbar sind. 


Die andere Möglichkeit 
u=xtin 


ng 
scheidet aus, weil nach den oben angeführten Sätzen von Haupt zu den Wurzeln der Gleichung 
Col (2 +im)= - Lojx=f(cos a) 
Lösungen der Differentialgleichung gehören, die alle Nullstellen haben. 
Wir haben also gesehen, daß die Gleichung 
al- Co] x = f(eos a) 
immer eine kleinste reelle (nicht negative) Wurzel cosa besitzt, und daß für diesen Wert von 
n, cos a immer eine nullstellenfreie periodische Lösung von (5) existiert, die wir etwa mit y, (f) 
el bezeichnen wollen. Von ihr zeigen wir, daß (8) tatsächlich eine Funktion mit dem Mittel- 
12 wert 0 ergibt, also eine kleine Schwingung um eine stabile Gleichgewichtslage darstellt. 
es 
it Day, E+T)=y,(t), ıst 
ir f n 
u. 22 If 
lım - -| y‚dt= T | y,dt. 
t> I t « « 
l- 0 0 
ll Zr u ' ' | 
Da die Funktion y,(f) das Vorzeichen nicht wechselt, so kann man sie so normieren, daß 
r 
m | 
T [w. dt=1. 
0 
). 
Dann ist plt)=tga-[y,(l)—1] 
n 
die unserer Bedingung genügende Lösung von (4). 
n. Ist dabei die kleinste Wurzel cosa<1, so haben wir Schwingungen um eine Gleich- 
gewichtslage «=F0 vor uns. Diejenige Frequenz des Aufhängepunktes, für die gerade 
” cosa—1 wird, ist die kritische Frequenz; denn für Frequenzen, für die cosa>]1, ist die 


Stabilisierung gar nicht möglich. 
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Kiniges sei noch über die Rolle gesagt, die die Dämpfung bei diesen Betrachtungen 
spielt. Führt man den Grenzübergang z>0 durch, so gewinnt man für das ungedämpfte 
System (2<=0) das Resultat, daß für dieses die kleinste Wurzel der Gleichung 


l f(cos a) 


genommen werden muß, zu der eine periodische Lösung gehört. Hätte man von vornherein 
z 0 gesetzt, so könnte man aus den vorhin angegebenen sehr einfachen Überlegungen nur 
schließen, daß Re a0. u müßte dann rein imaginär sein, könnte aber dafür alle imaginären 
Werte annehmen. Man hätte dann nicht nur periodische Lösungen zu betrachten, sondern 
alleemeiner stabile, d.h. solehe Lösungen, die stets innerhalb endlicher Schranken verbleiben '?). 
Man hat sieh bisher von vornherein auf periodische Lösungen beschränkt, ohne eine hin- 
reichende Begründung für das Ausschließen der anderen stabilen Lösungen (die, nebenbei 
gesagt, periodische Funktionen zweiter Art sind) angeben zu können. Hierzu ist bei ver- 
schwindender Dämpfung auch gar kein Grund vorhanden. Man würde dann nicht Gleich- 
vewichtslagen, sondern ganze Gleichgewichtsgebiete bekommen (etwa so, wie man bei der 
Behandlung der Schüttelerscheinungen elektrischer Lokomotiven, die auf eine homogene 
Hillsche Differentialgleichung führt, nicht diskrete Schüttelgeschwindigkeiten, sondern ganze 
Bereiche von solchen festgestellt hat)'’). Die richtige Stabilisierung des Vorganges wird 
also gerade dureh die Dämpfung bewirkt, was natürlich nieht hindern kann, in Fällen, wo 
sie sieher sehr klein ist, in der fertigen Lösung z gegen Null gehen zu lassen. 


Diskussion des Falles tg «=0. Wir haben früher den Falltg «=0 ausschließen müssen 
und holen hier seine Diskussion nach. Wenn wir zunächst cosa=1, a= 0" annehmen, geht 
(4) über in 


| 4 7 
q +2ung I ] 2 v0 
Mit oll)=e-«rtytt) 
( 1 
bekommt man y(t) re ] "+e’n), 0. 


Die Lösungen dieser Differentialgleichung haben nach dem Theorem von Floquet die Form 


.,d=ert.Feand und y,(t)=e-Frt.P(—2nnd). 
Dann ist u u 

9, d=en-»+Pmnt.F(2ant) und p,(t)=el-».-Mni.F(—Znnt), 
wobei Boj m = fleos 0) = fl). 


Für solche Werte von n, für die die kleinste Wurzel der Gl. 9) cosa>1 ist, wo also die 
Stabilisierung nicht möglich ist, ıst 


Golf u >kojz2, u>z. 


y,(f) wächst unbegrenzt, 9,(f) klingt ab. In diesem Falle gibt es wohl kleine Schwingungen 
um die Gleiehgewiehtslage «0, die aber labil sind. Für solche Werte von n, für die 
cosa<f1 wird, für welche die Stabilisierung also möglich ist, haben wir 


Goju<kojr, u<a. 
Beide Lösungen klingen ab. In diesem Falle sind also gedämpfte Schwingungen möglich. 
Nur ın dem Falle, daß 

bot v=bolx, 
wo also Goir=f(l), 


stellt die eine Lösung einen stationären Schwingungszustand dar. Dieser Fall aber ordnet 
sich der allgemeinen Theorie unter. 


Nimmt man jetzt «= 1S0° an, also cosa=-— 1, so wird man auf dieselbe Differential- 
eleichung geführt, die auch der Untersuchung von Schüttelerscheinungen elektrischer Loko- 
motiven zugrunde liegt. Bezüglich der Literatur dieses Problems sei auf Strutt l. e. hingewiesen. 
Die für uns in erster Linie in Betracht kommenden Arbeiten sind die von Meißner und 

13) Die analytische Durehreehnung dieses Falles scheint uns im allgemeinen auf große Schwierigkeiten zu 
stoßen. Sie würde vermutlich ergeben, daß die niehtperiodischen stabilen Lösungen keine kleinen Schwingungen mit 


dem Mittelwert 0 ergeben. : 
14) Freilich haben hier diese Gleichgewiehtsbereiche keine physikalische Bedeutung. 


\ „rr 
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neen Müller. Es genügt hier darüber folgendes zu sagen: Wenn die Frequenz des Aufhänge- 
OD .. r \ ® ® .. . . . 
pfte »unktes wächst, wechseln Gebiete ab, wo beide Lösungen stabil sind (bzw. abklingen, wenn 
eine Dämpfung vorhanden ist), mit solchen, wo die eine Lösung abklingt, die andere aber 
ınendlich anwächst (Schüttelgebiete). Getrennt werden je zwei dieser Gebiete durch 
requenzen, bei denen die eine Lösung rein periodisch ist, die andere abklingt (wenn das 
‚ a von + ı0 
rein System gedämpft ist). 
nur R 3 : 
Brane II. Spezieller Teil. 
lern Schwingungen des Aufhängepunktes. Bisher wurde von den Schwingungen des Aufhänge- 
13 . . Ä : ? i ; 5 j 
N). punktes nur vorausgesetzt, daß sie periodisch sind. Alle Aussagen, die wir bisher gemacht 
hin- haben, gelten daher für beliebige (periodische) Schwingungen. Die explizite Durchrechnung 
ıbei ist aber im allgemeinen Falle nicht durchführbar. Wir müssen eine spezielle Form der 
ver- Schwingungen des Aufhängepunktes annehmen, um mehr Einzelheiten von dem Verhalten des 
ich- Pendels zu erfahren. und wählen dazu eine Form, für welche die Hillsche Differential- 
der oleichung explizite lösbar ist, und die dabei die harmonische Schwingung hinreichend gut 
ene annähert. 
HIZe Den Arbeiten von Hirsch und Bock liegt die sehr naheliegende Annahme zugrunde, 
vird daß der Aufhängepunkt eine harmonische Bewegung vollführt, daß also 
wo 
n=4A,-sındnnt. 
In der vorliegenden Arbeit soll angenommen werden, daß die Zeit-Weg-Kurve des Aufhänge- 
BER punktes aus Parabelbogen zusammengesetzt ist, welche die Sinuskurve so annähern, daß dıe 
eht kinetische Energie der Ersatzbewegung ebenso groß ist wie 
die der ursprünglichen. Die Parabelbogen sind so aneinander- 
T 
oestückelt, daß sie an den Nahtstellen t=0o 5 (0-0, +1, +2....) 
samt Tangente stetig ineinander übergehen (siehe Abb. 2). Die 
Beschleunigung ist dann abteilungsweise konstant. 
Wenn, wie im folgenden überall I das Intervall 0-7 
en iu m 
„ T, II aber das Intervall „ T=t= T bedeutet, so ist Abb. 2. 
rm 5 - Ä 
N 2 
in I» 16.4 „| rl, 
(13) 
in II 16 mei z r\(1 Fr) 
h Dann ıst in I 5, 32 a n“,\ 
lie (14) 
in II 32 An? J 
Tı2 
| Die auf die Masseneinheit bezogene kinetische Energie | 7? dt ist 
.. A 0 
lie de 
im Falle der harmonischen Bewegung \ (2 an A, cos2ant® dt n’A’-n, 
0 
TI2 
. y y R - ; | ’qv r 32 .r 
im Falle der Ersatzbewegung 16.4 |; l— 21 di=— 4’ 
1. . . _ / +.) 
0 
Durch Gleichsetzen bekommt man 
> > R 
4 "— 39 “ Ss D . . . . . . . . . . . . (1»). 
et Die Amplituden der beiden Schwingungen sind nahezu gleich. 
l- Berechnung der Lösung. Setzt man (13) bzw. (14) in (7) ein, so bekommt man für y die 
4 Differentialgleichung: 
n. in 1 ze n’(2aocosa+ 2’—B32acosa)yV 
ıd bi , N br a ee 
in II yo n’(2aocosa+zx’+32acosa)y=V 
zu 
it A Eu 7 
| mit Am, 9m Be. ea Ne A 
16 


u 
Pl 
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a ist die „relative Amplitude“ der Schwingungen des Aufhängepunktes (bezogen auf die Ein- 
heit der Pendellänge), o eine dimensionslose Konstante, die dem Quadrate der Schwingungs- 
zahl des Aufhängepunktes umgekehrt proportional ist, » eben diese Schwingungszahl. 

(16) ist eine spezielle Hillsche Differentialgleichung, deren periodische Funktion ab- 
teilungsweise konstant ist. Solche Differentialgleichungen hat schon Meißner’) betrachtet 
(siehe auch bei Strutt, S. 39/40). Ihre Lösung kann explizite hergestellt werden, ebenso 
die zu ihr gehörige Gl. (9). 


Setzt man 


32 a0C0OSsa—+ 2° - H2acos a 4 h? \ 
| (18). 
und 32 a0cosa- x” +32 acos a 11° 1°) | 
so nımmt (16) die Form an: 
t+(2hn)”’y=0 in ı\ 
" ' (19). 
+(@kn)”’y=0 in I1J 
Da Ah sowohl als auch % konstant sind, ıst 
in I y(H=e,cos2?hnt—+e,sin2hnt, 
in II xd=c,cos2knt+c,sin2knt. 
"ür /= , T muß die Elongation sowohl als auch die Geschwindigkeit stetig sein 
4 Eı eu: en R 
(> T-0)=7(5 7+0), 42 r-0)=7(5 7+0)"), | m 
- ” - vr & Bus 9 _v), 


ferner nach (10) und (12): zT) = e*y(V), (T)=e* xI0) 


Die vier Rand- bzw. Stetiekeitsbedingungen ergeben ein homogenes, lineares Gleichungssy stem 


für die Koeffizienten e,.€,., €. €, 


c,cosh + e,sinh „cos k c,sin k 0 

ec, h-sinh-te,h-cosh + e,k-sink ec, k-cos k 0 ” 

6:8* c,„cos2k e,sin2%k 0 Dh 
e,h.e* Hc,k-sin2k —c,k-cos2k=0 





das eine nichttriviale Lösune nur dann besitzt. wenn seine Determinante 


cosh, sın bh. cosk, sin k 
hsınh, hcosh, ksınk, k cos k 

e” . 0 { cos? k, sın? k 

0 ne . ksin2k, kcos2k 


verschwindet. Daraus geht die Gleichung hervor, die der Gl. (9) der allgemeinen Theorie 
entspricht (Ak und Ak enthalten cos «a): 


Koi x --cosh-cosk 


I /h h: 
m .1 .1 . J») 
| k h s 1n hı sın I . . . . . . . (22). 


») 


Rechnet man aus (21) unter der Voraussetzung des Bestehens von (22) die Konstanten e, bis 
ce, aus, so bekommt man bis auf einen vorderhand unbestimmten Faktor b: 





‚A ’ 
ın I yd=b (sin h cos + h -sın h-sın 1) - cos 2hmndt 
| I /Ah J 
- eins = | |. sin h - sin k sin®hni\ 
| 2 \h k: J 
(25) 
FR A, 
nn II yN=b (e’ sin h + h sin :)- cos „(nt 1) 


h 


1; (ecosh  cosk)-sink(2nt 1)! 


15) E. Meißner: Über Sehüttelsehwingungen in Systemen mit periodisch veränderlicher Elastizität. Schweizer 
sauzeitunge, Bd. 72 (1918). S. 95 bis 99. 





16) Dieses A ist nieht mit dem früher ebenso bezeichneten Trägheitsradins zu verwechseln. 


I 
7), XI, T—N (ir)! 


(! T-+0| (! T) 
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), ist aus der Forderung 


. 
l . 
„|e et 
M 

u bestimmen, die für b immer einen endlichen von Null verschiedenen Wert liefert, wie wir 


‚siehe S. 239) ganz allgemein gezeigt haben. 


Der Fall verschwindender Dämpfung. Wie wir im allgemeinen Teil gezeigt haben, bekommt 
man die Gleichgewichtsbedingungen des ungedämpften Systems, wenn man in der Grund- 
sleiehung des gedämpften Systems z gegen 0 gehen läßt. Dabei geht (22) über in 


I /h k 
| cos h - cos k | + )- sin h ME (25 
2 \&k hı :- ID) 
und die Gleichungen (15) in 
h’— Sal o)ecosa und K=S8sii tja . . 5.047 IB). 


Man kann sich leicht überzeugen, daß (25) eine gerade Funktion von h und k ist, also 
eine eindeutige reelle Funktion von A’ und k°. In Abb. 3 ist die durch (25) dargestellte 
Kurve eingetragen in ein Koordinatensystem, dessen Abszisse h” und dessen Ordinate A? ist. 
Nach (26) ist 

k? I-+0o 
h? | o 


/)”7 
(ws ), 


Die in Frage kommenden Werte h, k bekommt man, wenn man die Kurve (25) zum 
Schnitt bringt mit der durch (27) dargestellten Geraden, die durch den Ursprung geht und 
mit der positiven h’-Achse einen zwischen 15° und 135° liegenden Winkel einschließt. 
Da nämlich O<n<x, ist © >6o>0, woraus unter Zuhilfenahme von (27) obige Behaup- 
tung folet. 


























7 " n 90° 
/ (XL: 
04 N N + + + pi } 1 80° 
/ 
104 a g° 
03 a + | . + 77 + N | 60° 
Il I i + + + 4 sn? 
m v7 50 
02 N + Y t + r T 1 40° 
I 1,71 - 30° 
ga —/4 11 20° 
ARER ER ER ER TR ER 
BEN || 
0 02 04 5 06 08 o wm 05 10 L r ” Z 
zZ I 


Abb. 3. Abb. 4. Abh. ». 


h?” und %? können dabei nicht beliebige Werte annehmen, wenn eine Gleichzewichtslage 
wirklich existieren soll. Nach (26) ist h”  k?=-16acosa: da aber |cos «a I,muß !h? 4° I6«. 
Die für uns in Betracht kommenden Punkte (h?, k?) liegen innerhalb des durch die Geraden 


B”’— 6a und A? - K’ er een 


begrenzten Parallelstreifens, dessen Breite noch von der relativen Amplitude abhängt (zu ihr 
proportional ist). Dieser Parallelstreifen liegt symmetrisch in bezug auf die Gerade h’ - 1° °0, 
die auch die Symmetrieachse der Kurve (25) ist. 


Für 5 =0(n=x) berührt die Gerade (27) unsere Kurve im Punkte W’=k’"-0. In 
diesem Falle muß auch cosa 0 sein, also «= 9%". Bei unendlich hoher Frequenz des 
Aufhängepunktes schwingt das Pendel um die horizontale Lage! (Dies wurde bereits von 
Hirsch mit anderen Methoden gezeigt.) 


Wächst jetzt o (nimmt die Frequenz des Aufhängepunktes ab), so wird sich der nicht- 
triviale Schnittpunkt der Geraden mit der Kurve (25) immer mehr vom Nullpunkt entfernen. 


Io” 
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Solange sich dieser Schnittpunkt innerhalb des durch (28) bestimmten Parallelstreifens bewegt, 
gibt es immer eine Gleichgewichtslage, die man aus (26) zu 


h?’ — k° 


COS A = 
15 a 


berechnet und die bei einer Frequenz auftritt, die sich aus 


h?— k? 
h? — k? 
ergibt. 

Diejenige Frequenz, für die cosa=1 wird, für die also der Schnittpunkt der durch (27) 
dargestellten Geraden mit der durch (25) dargestellten Kurve genau auf die Begrenzung des 
Parallelstreifens fällt, nennen wir die kritische Frequenz. Für diese Frequenz wird «=, 
das Pendel stellt sich vertikal. 

Überschreitet o diesen Wert o,, d.h. wird die Frequenz kleiner als die kritische, so 
haben (25) und (27) keinen nichttrivialen Schnittpunkt innerhalb des durch (28) bestimniten 
Parallelstreifens, das Pendel hat überhaupt keine Gleichgewichtslage (bis auf die triviale 
Ad IS"). 

Die kritische Frequenz hängt noch von a ab und kann bestimmt werden, indem man 
in h und k für cosa Eins setzt. 

Dann ist 

h=ySa(l—o,), k=iySa(l-+o,). 


Setzt man diese Werte in (25) ein, so bekommt man 


E On, ’ BR 
I =cosySa(l  0,)-Kolj BEREIT, le Sa(l- 0,)-SinySa(l+o, . (29). 
Aus (17) folgt 
wr; | 
ee 
7 (0, 


/u jedem a kann man aus (29) 6, berechnen und aus (30) die kritische Schwingungszahl n 
wenn noch die Pendellänge I! gegeben ist. 

Die für die kritische Frequenz maßgebende Größe o, wurde für einzelne Werte der 
relativen Amplitude berechnet und das Ergebnis in Abb. 4 eingetragen, die o,„ als Funktion 
von a zeigt (voll ausgezogene Kurve). Es zeigt sich, daß o, mit a wächst, daß also (siehe (30)) 
die kritische Frequenz bei wachsender relativer Amplitude abnimmt. 

Bei einer bestimmten vorgegebenen relativen Amplitude kann jetzt die Gleichgewichts- 
lage als Funktion der Schwingungszahl des Aufhängepunktes wie folgt bestimmt werden: 
aus Abb. 4 oder aus den Gl. (29) und (30) entnimmt man die kritische Frequenz, für die 
«VO ist. Gleichgewichtslagen sind vorhanden für alle » >n, und sind in der oben an- 
verebenen Weise zu ermitteln. Man kann sie dann als Funktion von » auftragen. In Abb. 5 
sind für a=0,2 und @a=0,5 die Gleichgewichtslagen als Funktion der dimensionslosen Größe 
» ] !/g aufgetragen. Wie alle Gleiehgewichtskurven haben auch die beiden aufgezeichneten 
die Gerade aM" zur Asymptote. 


Näherungsweise Berechnung der kritischen Frequenz und der Gleichgewichtslage bei raschen 
und kleinen Schwingungen des Aufhängepunktes. Es wurde bereits darauf hingewiesen, daß 
bei sehr hohen Frequenzen des Aufhängepunktes das Pendel nahezu horizontal schwingt, 
daß also eos a sehr klein geren Eins ist. Nach (26) müssen dann auch A? und A? sehr klein 
sein. Andererseits sind bei einer geringen relativen Amplitude a, h? und %k? stets sehr klein, 
was für Werte auch cosa annehmen mag, da jedenfalls o<1, wie aus Abb. 3 oder 
auch dem Diagramm der kritischen Frequenzen ersichtlich ist. In beiden Fällen, die alle für 
uns interessanten umfassen, sind A? und Ak? klein, und man kann eine erste Näherung erhalten, 
wenn man die rechte Seite von (25) nach steigenden Potenzen von h und k entwickelt und 
die Glieder 3. Ordnung in 4? und A? vernachlässigt. So bekommt man 


| LEN 
| 1 —- ("+ +5, Bh'+ 05° 7 +3k’) +0 (h*, k*). 
Wenn man (26) in diese Gleichung einführt, bekommt man als erste Näherung für die Gleich- 
vewichtslage 


(31). 


EOS (A 
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egt, Wir wollen zunächst zeigen, daß (31) genau die Hirschsche Formel liefert, die also 


als erste Näherung der Lösung der von uns aufgestellten Gleichung erscheint. Nach (17) ıst 


A 2» 7) | Il 
ı=Tı (7 BA Tim 


Berücksichtigt man noch, daß nach (15) 


. 3° » 
Po. 
27) so wird 
. gi 
les cos Aa 3_3 =— (la). 
-() a 2, m“ 


Wenn man in der Formel (14b) auf S. 485 in der Arbeit von Hirsch 9-0 setzt, wie 


so n ; “17: .- , 
es unserem Falle des in der Schwererichtung oszillierenden Pendels entspricht, und für A 
en ur GG 4” ’ 
vemäß (13) auf S. 4 
ale ur 
Br. SP 
K= F: 
an J 
einführt, so geht die Formel (14b) der oben zitierten Arbeit in unsere Formel (31a) über. 
Es zeigt sich also, daß die Voraussetzungen der Arbeit von Hirsch sind: entweder 
kleine Amplitude oder aber hohe Schwingungszahl des Aufhängepunktes, Voraussetzungen, 
von denen man sich leicht befreien kann; dann aber auch die Voraussetzung, daß das Pendel 
um eine Gleichgewichtslage kleine Schwingungen ausführe, eine Voraussetzung, die auch 
9) unserer Arbeit zugrunde liegt und deren Aufgeben das Problem wesentlich komplizieren würde. 
Zum Vergleich der Hirschschen Resultate mit unseren genaueren Ergebnissen sind 
die kritischen Frequenzen, wie sie sich aus (31) ergeben, in Abb. 4 eingetragen (strichlierte 
)). (Gerade). Für cosa=1 bekommt man 
2a 
&, =9 
+.) 
er die Gleichung einer Geraden, die die richtige Kurve der kritischen Frequenzen in a0, 
IN s—=0 berührt und sie für kleine a und o gut annähert. Aber schon bei a = 0,2 ist der Fehler, den 
)) man begeht, wenn man die voll ausgezogene Kurve durch die striehlierte Gerade annähert, 
merklich groß. 
'S- - 
B Berücksichtigung der Dämpfung. Ist die Reibung nieht Null beziehungsweise so klein, 
” daß man sie vernachlässigen kann, so hat man an Stelle der Gl. (25) 
n- | 
Is‘ Ze l e7 I: n r 
Go] z =cosh -cosk 7 + n |sinh - sink WE er rer  . (22), 
3€ ; pn! v 
»n ‘ u ® ) 
und an Stelle von (27) gewinnt man aus (18) 
.& BE x” | +0 (32) 
) > . . . . . . . . . . . . . . ı—]. 
n Ih? —+ 7° 1 —o 
E Der Schnittpunkt der Kurve (22) mit der Geraden (32) ergibt h’ und k’; dann ıst 
N h? I? 
COS A R 
: 16a 
T | . ; R b a 
£ Aus (22) geht hervor, daß unsere Kurve den Punkt W" ki, d.h. Wh" * ent- 
1, hält, und daß sie symmetrisch in bezug auf die Gerade h?  Kk? ist. 
Offenbar ist (32) die Gleichung einer Geraden, in der h’-, A*-Ebene, die durch den 
Scheitel der Kurve (22), nämlich durch den Punkt h?”=Kk? = eeht und mit der posı- 
tiven A?-Achse einen zwischen 5’ und -— 135° liegenden Winkel einschließt. 


- Ist x zahlenmäßig gegeben, so kann man die zugehörige Kurve (22) zeichnen und die 
Ermittlung der Gleichgewichtslagen, der kritischen Frequenzen geht ebenso vor sich wie 
früher. Für #2 ist die Kurve in Abb. 3 wirklich eingezeichnet und in Abb. 5 die Gleich- 
gewichtslagen für z=2, a 0,5. 











Ber. 
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Näherungsrechnungen für kleine Reibung. Ist x klein, so befindet sich der Scheitel der 
Kurve (22) in der Nähe des Ursprunges. Für große Frequenzen oder kleine Amplituden 
befindet sich der in Betracht kommende Schnittpunkt in der Nähe des Scheitels, also wieder 
in der Nähe des Ursprunges, so daß man A’ und k° klein voraussetzen darf. Entwickelt 
man wieder (22) nach Potenzen von h und k und vernachlässigt die Glieder 3. Ordnung in 
h" und A”, so bekommt man 


Koi x | (h”’—+&k”) nah - 10Rh?’k?+3K’) +0(h®, k*). 
Setzt man in diese Gl. (18) ein und entwickelt die linke Seite nach Potenzen von x, so be- 
kommt man bis auf Glieder höherer Ordnung 
(b-+-x°)o m 
COS A: en 
la 

eine Formel, die für < © in (31) übergeht und zeigt, daß bei kleiner Dämpfung der Einfluß 
derselben numerisch gering ist und daß dieser Einfluß in einer Verkleinerung von a besteht. 
Die kritische Frequenz wird durch die Reibung erhöht. 


Näherungsrechnungen für große Reibung. Der Schnittpunkt der Kurve (22) mit der 
Geraden (32) fällt bei großen Frequenzen oder bei kleinen Amplituden des Aufhängepunktes 


in die Nähe des Scheitels Ah? A’ Ra, so daß in diesem Falle 4A’ —+x°” und 4k’—+x° 
klein sind. Ist die Reibung groß, so kann man annehmen, daß 

Ih’? 32 ll o)acos a z 
une 

Ik’+x° 32 (1 +0)acosa a". 


Entwiekelt man die rechte Seite von (22) in die Taylorsche Reihe ın der Umgebung der 


| 2 | 
Stelle h?’= k? = und vernachlässigt wieder die Glieder 3. Ordnung in [h’-+ | #); 
| a 
Ik: | x’), so kommt man auf die Gleichung 
Be ,,»e\Sin«, | 4 “’\"&vjx 
+ —)+He +7) + PR +Zz)tHR+G ) 23 


Setzt man hierin 


h?- x" —=S(0 l))acosa 
und 
c#; 
L? a: S(o-+Il)acos a 


und vernachlässigt man höhere Potenzen von o (entsprechend der Tatsache, daß für hohe 
reuenzen o sehr klein ıst), so kann man als erste Näherung bekommen: 


HM 


0» 
vos a ==, 31 
Sa(?-+x) .. 
oder, wenn x 5: 
0x” 
cos a | Ih DD 2 ed 18 „2; GES Be. BE SEe, PRESENT EB ee Er 
Sa e ) 


Die Formel zeigt qualitativ dieselbe Abhängigkeit von o und a wie (31). Für große z, 
für die sie allein gültig ist, zeigt sich wieder, daß durch die Dämpfung cos a vergrößert, 
a also verkleinert wird. Der Einfluß ist hier ganz beträchtlich. 

Die kritische Frequenz ergibt sich aus (34a), indem man cos a1 setzt 
Sa 


4° 





New 
Mech 
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ıden 
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he- 
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(jerenüber dem Werte der kritischen Frequenz bei ungedämpftem System zeigt es sich, daß 
& beträchtlich abnimmt, also die kritische Frequenz, die geringste Frequenz, bei der diese 
Stabilisierung noch möglich ist, stark zunimmt. Ja, wenn die Reibung über alle Grenzen 
wichse, würde auch die kritische Frequenz, unterhalb deren sich die hier behandelten Er- 
scheinungen gar nicht zeigen, unbegrenzt zunehmen, die Stabilisierung wäre gar nieht möglich. 
‘Selbstverständlieh ist hier die Stabilisierung von Schwingungen gemeint. Daß soleh ein 
Pendel in jeder Stellung in statischem Gleichgewicht verharren kann, ist einleuchtend.) 


Diskussion des Falles tg «= 0. Der Fall tga=0 ist ein Ausnahmefall der allgemeinen 
Theorie und muß auch bei diesem Beispiel speziell diskutiert werden. 
Wenn wir zunächst eosa-==1, a=-0" nehmen, so gehen die Gl. (18) und (19) über in 


Ih"?=32ao+ x: — 3a 
IE=-B»ao+xr’+ 3a, 
bzw. in | 7+(2hn)’y- 0 
in II y+(2kn)’y=Vd. 


Es ist in diesem Falle 
h? k? 16a 


| 
h? + k? l6ao " 
Da do <on x, so Ist zu >w’-+K >, 
h”+k?: „ x” ist die Gleiehung der Tangente unserer Kurve (22) in ihrem Scheitel. 


(In Abb. 3 strichliert.) Die Gerade h? - k?= 16a steht senkrecht auf dieser strichlierten Ge- 
raden und schneidet diese sowohl als auch unsere Kurve (22), erstere für 0-0, letztere 
für 0:=64,. 

Wächst nun so von Null, so wird sich der auf A?’ k’= 16a befindliche Punkt immer 
mehr von der strichlierten Tangente entfernen. Solange 6 < o,, befindet er sich im Gebiet 
zwischen Tangente und Kurve, wo sicher Re u< x, beide Lösungen klingen ab, das Pendel 
kann keine stabilen Schwingungen um die Nullage ausführen. Das sind nämlich gerade die 
Schwingungszahlen oberhalb der kritischen, und das Pendel hat eine von Null verschiedene 
Gleichgewichtslage, die wir früher berechnet haben. Wenn der wandernde Punkt gerade ın 
den Schnittpunkt der Geraden h? — k’=- 16a mit der Kurve (22) fällt, ist o = o,, wir haben 
die kritische Frequenz vor uns, das Pendel führt um die Nullage stabile Schwingungen aus. 
Überschreitet o diesen Wert, liegt der Punkt jenseits der Kurve, so ist Neu >z, die eine 
l,ösung klingt ab, die andere wächst über alle Grenzen, wir haben labile Schwingungen vor 
uns, für alle Frequenzen, die kleiner sind als die kritische. 

Über «= 180° siehe das im allgemeinen Teil Gesagte. 


Zusammenfassung. Wir haben gezeigt, daß ein Pendel, dessen Aufhängepunkt in vertikale 
Schwingungen versetzt wird, um eine schräg nach oben gerichtete Mittellage Schwingungen 
ausführen kann, wenn die Frequenz des Aufhängepunktes einen Schwellenwert überschreitet, 
den wir als die kritische Frequenz bezeichneten und von der wir zeigten, daß sie von der 
Amplitude der Schwingungen des Aufhängepunktes und der Dämpfung des Systems abhängt. 
Es wurde eine Methode entwickelt, die es gestattet, bei vorgegebener Amplitude der Oszilla- 
tionen des Aufhängepunktes und Dämpfung des ganzen Systems die kritische Frequenz sowie 
die Gleichgewichtslage für jede Frequenz oberhalb derselben zu berechnen. Näherungsformeln 
für sehr kleine und sehr rasche Schwingungen des Aufhängepunktes ergaben bei verschwin- 
dender Dämpfung die von Hirsch aufgestellte Formel und zeigten, daß sie für gedämpfte 
Systeme eine zu kleine kritische Frequenz und einen zu großen Gleichgewichtswinkel ergibt. 
Die für gedämpfte Systeme geltenden Näherungsformeln zeigten noch, daß mit unbegrenzt 
wachsender Dämpfung die kritische Frequenz über alle Grenzen wächst, eine Stabilisierung 
also nieht mehr möglich ist. 





Zum Schlusse sei mir gestattet, Herrn Professor Weyrich meinen Dank auszusprechen, 
der nicht nur die Anregung zu dieser Arbeit gab, sondern in eingehenden Besprechungen mir 
auch Gelegenheit bot, über wichtige Punkte Klarheit zu verschaffen. 


Brünn, Mathematisches Seminar der Deutschen Technischen Hochschule, 27. März 1935. 
+.) )* 
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Konstruktion eines Fräsers, der eine gegebene 
Schraubenfläche erzeugt. 
Von Othmar Baier in München. 


A. Erklärung der Konstruktion. 


1. Vorbemerkungen. Eine Schraubenfläche [A] sei vorgegeben. Ihre Achse sei a und 


ihre Ganghöhe h. Die Schraubentläche [A] kann stets durch Verschraubung einer auf ihr 


liegenden Kurve (C) erzeugt werden. Gesucht ist der Meridian desjenigen Fräsers [F] (d. i. 
derjenigen Drehfläche), der bei Drehung um die noch beliebig vorgebbare Achse b in ein 
Werkstück die gegebene Schraubenfläche [A] fräst, wenn sich das Werkstück bei einer 
Volldrehung um die Achse a um die Ganghöhe Ah in Richtung der Achse a verschiebt'). 

Der Fräser [F] wird die Schraubenfläche [A] in den Punkten einer Kurve, der Ein- 
eriffslinie, berühren. Falls es gelungen ist, die Eingriffslinie zu ermitteln, kann man aus 
ihr leicht die Form des Fräsermeridians konstruieren, da der Fräser [7] durch Drehung der 
Kineriffslinie um die Achse 5b entstanden gedacht werden kann. 

Es soll eine Methode angegeben werden, mit der die Eingriffslinie rasch punktweise 
ermittelt werden kann. 

Die Normalen der Schraubenfläche [A] in den Punkten der Eingriffslinie müssen als 
Normalen von [F| die Fräserachse b treffen. Es handelt sich also darum, diejenigen unter 
den Normalen der Schraubenfläche [4] ausfindig zu machen, welche die Fräserachse b 
schneiden. Die Gesamtheit ihrer Fußpunkte auf der Schraubenfläche [A] bildet dann die 
Kıngrıffslinie. 

2. Hilfsbetrachtung (vgl. Abb. 1). Wir konstruieren zunächst in einem beliebigen Punkt P 
der Kurve (€) die Normale » «der Schraubenfläche [A4]’). Wird die Normale » längs der 














Abh. 2. 


Bahnschraubenlinie des Punktes P verschraubt, so kommt sie auch in eine Lage n, in der 
sie «lie Achse 5b schneidet?). Dieser Schnittpunkt sei mit 9 bezeichnet. Der Fußpunkt P 
der Normalen » ist dann nach dem Obigen ein Punkt der Eingriffslinie. Es wäre demnach 
lie dureh Verschraubung der Normalen » entstehende Schraubenftläche 5 mit der Fräserachse 
zum Schnitt zu bringen. Um weitere Punkte P,, P,,.... der Eingriffslinie zu erhalten, 
müßte man jeweils die durch Verschraubung der zugehörigen Normalen n,, n,. 
entstehenden Sehraubenflächen mit der Achse 5b schneiden. 

', Lösungen für das analoge Problem bei Zylinderflächen, Drehflächen und, falls die Fräserachse b der Schrauben 
achse a parallel ist, bei Sehraubenflächen hat S. Finsterwalder in seinen Vorlesungen über Darstellende Geometrie 
an der Technischen Hochsehnle München angegeben. Dadurch wurde ich zur Behandlung der gestellten Frage angeregt. 

-) Die Normale »n wird als Lot im Punkt P auf die Tangentialebene der Schraubenfläche [.4| gefunden. Die 


"angentialebene in P ist bestimmt dureh die Tangente an die Kurve (C) in P und die Tangente an die durch P gehende 
Bahnschraubenlinie. 


>) Man wählt, falls verschiedene solche Lagen existieren, zweckmäßig diejenige, die der Ausgangslage am 
nächsten ist. Falls keine solehe Lage existiert das kann in Ausnahmefällen vorkommen ‚kann bei der gegebenen 
lage der Achse b kein Fräser die ganze Sehraubenfläche erzeugen. 
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3. Konstruktionsprinzip (vgl. Abb. 2). Die obige Konstruktion der Eingriffslinie läßt sich 
‚och sehr weitgehend vereinfachen: Anstatt die Schnittpunkte der Achse b mit den durch 
Verschraubung der einzelnen Normalen n,, n,, ... . entstehenden Schraubenflächen aufzu- 
:uehen, kann man ebensogut die Schnitte der Normalen »,, n,,.... mit der dureh Ver- 
schraubung der Achse b entstehenden Regelschraubenfläche [B] (Schraubungsachse a, Gang- 
höhe A) aufsuchen. Die Normale n schneide diese Schraubenfläche [B] in einem Punkt @. 
Durch diesen Punkt 0 geht dann eine Gerade b* der durch Verschraubung von b entstandenen 
Schraubenfläche [B]. Man schraubt nun die Gerade b* so weit zurück, bis sie wieder mit b 
usammenfällt. Übt man dieselbe Schraubung auf den Punkt P und die Normale » aus, so 
selangen sie in die gesuchte Lage P, n. Diese Abänderung der Konstruktion hat die ver- 


»infachende Folge, daß man die Normalen »,, n,,..... mit der einzigen Schraubenfläche [BP] 
‚ıı schneiden hat, um die gesuchten Punkte P,, P,, .... der Eingriffslinie zu finden. 


4. Praktische Ausführung der Konstruktion. Die praktische Ausführung der Konstruktion 
vestaltet sich zweckmäßig folgendermaßen: Es wird passend die Grundrißtafel senkrecht zur 
Schraubenachse a, die Aufrißtafel parallel zur Fräserachse b angenommen. Dann graduiert 
man die Normale », d. h. man teilt ausgehend von einer beliebigen zur Grundrißtafel parallelen 
bene & die Ganghöhe h in eine gewisse Anzahl » von gleichen Teilen und sucht diejenigen 

us | h h 
Punkte der Geraden » auf, deren jeweilige Höhen über &e 0, + 8 2 ea sind. Dazu 
ist eszweekmäßig, die Aufrisse einer Reihe von zur Grundrißtafel parallelen Ebenen zu zeichnen, 


h h i 
die von e die Abstände + —, +2 ren besitzen, dann deren Schnittpunkte mit »'"*) zu 
v ’ u 


ermitteln und hieraus die Grundrisse der Schnittpunkte auf »’ zu finden. Die so gefundenen 
Punkte seien durch ihre Höhenmarken bezeichnet. Man konstruiert nun den Schnitt der 
dureh Verschraubung der Achse b entstehenden Schraubenfläche [PR] mit der Ebene & ((uer- 
schnitt einer Regelschraubenfläche). Durch den Grundriß a’ der Schraubenachse a und den 
Grundriß des Sehnittpunktes von b mit e zieht man eine von a’ ausgehende Halbgerade y 
und teilt den Vollwinkel mit dem Scheitel a’ von g ausgehend in » gleiche Teile und schreibt 


I Ir 


zu den einzelnen Halbstrahlen den mit g eingeschlossenen Winkel: 0, + 3 


v v 
Dabei ist, je nachdem die vorliegende Schraubenfläche rechts- oder linksgängig ist, der Winkel 
von g aus im Zeigersinn oder Gegenzeigersinn positiv zu zählen. Die Schnittkurve samt 
Winkelteilung und Bezifferung überträgt man auf Pauspapier. Die horizontalen Querschnitte 
der durch Schraubung von b entstehenden Schraubenfläche [B] sind kongruent, und man erhält 


h 

den Querschnitt mit der horizontalen Ebene, die von e den Abstand 1 hat, indem man die 
) - 

Pause so lange um a’ dreht, bis der Halbstrahl mit der Marke l auf den Strahl mit der 


Marke 0 zu liegen kommt. Um nun den Schnitt der Geraden n mit der Schraubenfläche [3] 
zu erhalten, dreht man die Pause so lange, bis der Schnitt von »’ mit der gedrehten Schnitt- 
kurve die gleiche Marke (Faktor von A, bzw. 27) zeigt, wie der über 9 (Marke 0) liegende 
Halbstrahl. (Dies kann durch Interpolation rasch und mit genügender Genauigkeit erreicht 
werden.) Dieser Schnittpunkt ist nach dem Obigen der Punkt @. Man dreht jetzt den 
Punkt 9’ zurück, bis er auf b zu liegen kommt, und erhält damit den Punkt 9°). Des- 
eleichen dreht man den Punkt P’ um den Winkel 0 «a’ 9 zurück und erhält den Grund- 
"iß P’ des gesuchten Punktes P der Eingriffslinie. Nach bekannten Methoden kann nun der 
Aufriß der Eingriffslinie und aus den beiden Rissen der Eingriffslinie der Meridian des 
Fräsers punktweise konstruiert werden. Da 9 P auch Normale des Fräsers ist, kann man 
in jedem konstruierten Punkt des Fräsermeridians dazu noch dessen Tangente angeben). 


Ist speziell die Fräserachse b parallel zur Schraubenachse a, so wird die Schrauben- 
äche [B] zu einem Drehzylinder, dessen Schnittpunkte mit den Normalen »n,, n,, .... sofort 
konstruiert werden können. Die weitere Konstruktion verläuft dann wie oben. 


%, (Grundrisse sind durch einen Strieh, Aufrisse durch zwei Striche bezeichnet. 

5), Durch den achsennächsten Punkt des Querscehnitts werden die Punkte des Quersehnitts in zwei Klassen ge- 
teilt. Die Erzeugenden der Schraubenfläche |B]| dureh die Punkte der ersten Klasse treffen die Kehlsehraubenlinie 
von [RB] unterhalb, die dureh die Punkte der zweiten Klasse oberhalb der Quersehnittebene. Je nachdem @ zur ersten 


oder zur zweiten Klasse gehört, muß auch @ auf b unterhalb oder oberhalb des Treffpunkts von b mit der Kehl- 


schraubenlinie liegen. Dadureh ist Q eindeutig bestimmt. 

%) Man wird dureh Wahl einer geeigneten Lage der Achse b zu der zu fräsenden Schraubenfläche [A] eine 
möglichst günstige Form der Eingriffslinie zu erreichen suchen. Es kann unter Umständen eintreten, daß der Fräser- 
neridian sowohl Selbstüberschneidungen als auch Spitzen aufweist. 
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B. Durchführung der Konstruktion des Fräsers an einem Beispiel. 
Zur Erläuterung der Konstruktion wurde als Beispiel für die Schraubenfläche [A] ein 
Spiralbohrer gewählt. Die beiden Aufnahmen (Abb. 3a, b) zeigen den Spiralbohrer und ein 


f 


> 





nu! 3 
Abb. 3a. Abb. 3b. 


Holzmodell des Fräsers. In Abb. 3a ist die Eingriffslinie in die Bohrernut weiß eingezeichnet. 
Abb. 3b zeigt den Fräser im Eingriff. 

In Abb. 4 ist die Konstruktion eines Punktes der Eingriffslinie durchgeführt. Im Grund: 
ri ist der Querschnitt des Bohrers (Maßstab 1:2,5) mit der horizontalen Ebene, die den 
kürzesten Abstand der Schraubenachse a und der Fräserachse b enthält, gezeichnet. Der 
Aufriß dieser Ebene ist mit O0 bezeichnet. Sie spielt die Rolle der Ebene & (vgl. oben!). Zur 
Konstruktion verwendet man als erzeugende Kurve (C) für einen Teil der Schraubenfläche 
zweckmäßig die Schneidkante des Bohrers. In der Abb. 4 wurde die Schneidkante in den 
Punkt P des Querschnitts geschraubt und ihre Ver- 
\ / längerung in dieser Lage mit s bezeichnet. Zur Kon- 

struktion der Normalen » der Schraubenfläche [A] im 
; Punkt P verschiebt man sowohl die Tangente an die 
Bahnschraubenlinie in P als auch die Tangente an die 
Kurve (C) in P — diese fällt hier mit s zusammen 
parallel, bis beide Gerade durch die Spitze S des Rich- 
tungeskegels gehen (der Abstand des Punktes S von der 
Ebene O0 ist l/ 2m). Diese beiden Geraden durch S be- 
stimmen dann eine Ebene, die parallel zur Tangential- 
+, ebene des Punktes P’ ist. Man sucht die Spurparallelen 
SD dieser Ebene auf und konstruiert in bekannter Weise 

die Normale ». Ein Teil des Querschnitts der Schrauben: 

läche [B] mit der Ebene 0 ist in der Abbildung mit den 
3. zugehörigen Marken eingetragen. Der Faktor 27/16 ıst 
bei den Querschnittmarken weggelassen. Die Aufrisse 


. r h h _ 
der Ebenen mit den Abständen + ,.. +2,45: +94; 
16 16 16 


IS) 
— 








r 


von der Ebene O0 sind eingezeichnet. Der Faktor h/ 16 
ist bei den Höhenmarken ebenfalls weggelassen. Im 
Grundriß ist die Lage des gedrehten Querschnitts der 
Schraubenfläche [B], der den zu P gehörigen Punkt W 
liefert mit den zugehörigen Marken gestrichelt einge» 
{ zeichnet. Durch Drehung erhält man in der oben an- 
b “u | X verebenen Weise aus den Punkten P’ und 9 die Punkte 
- Nr P’ und © und aus P’, P” und P’ den Punkt P”. Die 
wahre Gestalt des Fräsermeridians ist der Umriß des 
Ks 7 37 Fräsers im Aufriß. Den zu P gehörigen Umrißpunkt 7, 
2 2 bestimmt man, indem man die Ebene (LP b) um b dreht. 

bis sie parallel zur Aufrißtafel ist. Bei dieser Drehung 

kommt der Punkt P in die Lage P,. Die Senkrechte auf P, 0’ in P, ist die Tangente des 
Meridians in P,. Um weitere Punkte der Eingriffslinie zu finden, wurden außer der verschraubten 
Schneidkante noch Punkte des Quersehnitts benützt. Die Risse der sich so ergebenden Ein- 
eriffslinie sind punktiert eingezeichnet, der Fräsermeridian ist durchgezogen ’). 439 





wo 


« 


) Da die Aufnahmen aus verhältnismäßig geringer Entfernung erfolgten, scheint bei Abb. 3b infolge der 
perspektivischen Verzerrung der Fräser zu weit über die Bohrernut herauszuragen. Dies ist jedoch nicht der F all, 
wovon man sieh dureh Abb. 4 überzeugen kann. 
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KLEINE MITTEILUNGEN 


Zum Landestoß von Seeflugzeugen. |ı 
iner inhaltreiehen Arbeit über Stoß- und Gleit- 
vorgänge an der Oberfläche von Flüssigkeiten (diese 
Jeitsehrift Bd. 12, S. 193, 1932 bzw. auch Verhandl. 
d. 3. intern. Kongresses f. techn. Mechanik Bd. 1. 
S, 215. 1930 und Zs. f. Flugtechnik u. Motorluftsch. 
Bd. 22, S. 1. 1931) behandelt H. Warner an her- 
‚orragender Stelle den Landestoß von >Seetlug- 
euren (diese Zeitschrift 1. e.. S. 205). Zu der In- 
teeralgleichung, die nach H. Wagner die ent- 
cheidende Rolle für dies Problem spielt, möchte 
ch die folgende Bemerkung machen. 

In den Bezeichnungen der angeführten Arbeit 
st 7, (x) die Gleichung des Bodens des Schwimmers, 
ler als zylindrischer Körper behandelt wird. Der 
\nfangspunkt des Koordinatensystems ist dabei 
in den tiefsten Punkt des symmetrischen, gekielten 
schwimmerbodens gelegt: .r ist die horizontale Ent- 
(ernung (parallel der ursprünglichen. Wasserober- 
tläche vor dem Eintauchen des Schwimmers), 7, ist 
die vertikale Entfernung von diesem Punkt. V sei 
die vertikale Eintauchgeschwindigkeit des Schwim- 
mers,. die im allgemeinen von der Zeit f abhängig 
ist. ec ist die halbe Breite der in einem bestimmten 
Augenblick gerade benetzten Fläche des Schwimmers. 


! 


k, [ r [3 [3 . * 
r also die halbe Verbreiterungsgeschwindigkeit 


z V 
der Aufschlagfläche. Dann ist mit u Et 
7b 
x dt 
" ul(e)de 
n,\£) = 
ce? 
Via 
« DB” 


Wenn man (ec) aus dieser Volterraschen Integral- 
rleiehung erster Art ermittelt hat. so kann man. 
wie H. Waener gezeigt hat, alles Weitere, insbe- 
sondere die Kraft auf den Schwimmerboden. leicht 
ausrechnen, 

H.Waener eibt nun die Lösung dieser grund- 
lerenden Integralgleichung in dem Spezialfall an, 
daß nn» (r) in Form eines Polynoms in x gegeben 
ist, wobei dann #(c) ebenfalls ein Polynom wird, 
dessen Koeffizienten durch Einführen eines ent- 
sprechenden Ansatzes in die Integralgleichung be- 
stimmt werden (. e. Gl. (36a) (36b)). Wenn die 
Kenntnis einer solchen Lösung praktisch auch oft 
renügen wird. so wird trotzdem eine allgemeine 
l,ösunge nicht ohne Interesse sein. 

Die vorliegende Integralgleichung ist nahe ver- 
wandt mit der Abelschen Integralgleichung, so 
daß sich ihre allgemeine Lösung tatsächlich leicht 
anzeben läßt (ef. Frank- v. Mises: Die Differen- 
tial- u. Integralgleichungen d. Mechanik u. Physik. 
Bd. 1. S. 484, Braunschweig 1930). Es ergibt sich 

> 
"1 | )1 
ulc) = bio o, u | E 6 \X \ G L\ 
“N Jäzl z |\yr—al 


0 


-) 


“, ist dabei der Winkel des Schwimmerbodens 
veeren die r-Achse im Punkte r==0. Im Falle 
eines geradlinig gekielten Bodens ist übrigens 


77 : 
:eonst =tra,, so daß dann in der obigen 


lösung das Integral wegfällt. 

Von Praktikern ist mir erklärt worden, daß recht 
komplizierte Bodenformen n,, (x) aus konstruktiven 
(“ründen in Betracht kommen könnten, so dab die 
srößere rechnerische Bewegungsmöglichkeit, die 
dureh die Lösung für beliebige 1 (z) gegeben ist, 
für die Ausnutzung der Wagnerschen Theorie 
willkommen sein dürfte. 
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Über die durch Krümmung _ steifge- 
machten Meßbänder. \iuclht man ein derartiges 
Meßband umzubiegen, so wie es in der beistehenden 
Abb. angedeutet ist. so kann man am Meßband 
deutlich fünf verschieden zeformte Teile wahr 
nehmen. In den Teilen AB und EF bleibt die 
Krümmung des Querschnittes erhalten, und diese 
Teile nehmen offenbar fast gar keine Deformations 
energie auf. Die Stücke BC und DE stellen 
„Übergangsstücke“ dar. und im Stück CD ist die 
Krümmung im Querschnitt vollkommen verschwun- 
den, dagegen tritt in diesem Stück und im wesent- 
lichen nur: in diesem Stück eine Krümmung des 
lL,ängsschnittes auf. Diese Tatsache müßte man 
bei einer ausführlichen Behandlung des Gegen 
standes, etwa aus dem Minimumprinzip der Form 
änderungsarbeit mathematisch erschließen, doch er- 
scheint sie auch dureh obertlächliche Betrachtungen 
in dieser Richtung plausibel. Im Stück CD ist nun 
der Krümmungsradius annähernd konstant und an- 
nähernd gleich dem Krümmungsradius des Quer- 
schnittes im ursprünglichen Zustand. Den Grund 
hierfür kann man sich durch folgende kleine Über- 
lezung klarmachen. 











Nehmen wir an. die Deformationsenerzie der 
Stücke AB und EF sei gleich Null und die Defor- 
mationsenergie der Übergangsstücke BC und DE 
sei unabhäneiz von dem Winkel «a, den die Stücke 
AB und EF miteinander einschließen. Die Defor- 
mationsenergie des Stückes U D kann man sieh aus 
zwei Bestandteilen zusammengesetzt denken. 

1. Aus der Aufbiegarbeit Aj. d. h. der Energie, 
die notwendig ist, um die Krümmung des Quer- 
schnittes rückgängig zu machen. Dieser Bestand- 
teil ist offenbar proportional der Länge dieses 
Stückes. 

2. Aus der Umbiegarbeit Ay. d. h. der Energie, 
die notwendig ist. um die Krümmunz im flachen 
Band (bei dem man sieh die Krümmung in der 
Querrichtung bereits rückgängig gemacht zu denken 
hat) in der Längsrichtung zu bewirken. Bedeutet 
ft; den Krümmungsradius des Querschnittes im 
ursprünglichen Zustand und ZAj den Krümmungs- 
radius im Längsschnitt von ED, z und m zwei 
ortsunabhängige Größen. und bedeutet g den Winkel, 
den die Tangente in einem beliebiren Punkt. der 
die Mittellinie des Stückes ED mit der Mittellinie 
des Stückes AB einschließt. so ergibt sich für das 
3jorenelement im U D 


d Ss — It d G- 


und es ergeben sich für A; und Ay; Ausdrücke von 
der Form 
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Dann folet aus 
A ! Ay -> Min. 


daß Ay konstant sein muß. ferner dureh Diffe- 
rentiation nach Ajr. dab Ay= Ay sein muß, Stellt 
man nun die Ausdrücke für die Aufbierarbeit und 
die Umbiegarbeit durch Doppelintegrale dar, so 
eeht der eine aus dem anderen hervor, indem man 
/?; und Zj; vertauscht. Aus der Gleichheit von 4) 
und „Ij; kann man somit auf die Gleichheit von Rj 
und Zr schließen. 

Prag. Br. F 


unk. 380 


Ein ebenes Diagramm für ternäre Ge- 
mische. lin ebenes Diagramm zur Darstellung 
der Zustandskurven ternärer Gemische. Anwendung 
auf die Entropievermehrung beim idealen ternären 
(Gemisch. Konstruktion dieses Diagrammes. Be- 
stimmung des Höchstwertes. Beziehung zwischen 
intropievermehrung und theoretischer Trennungs- 
arbeit bei Gasgemischen Benutzung des Dia- 
erammes zur Ermittlung der Trennungsarbeit. 

An anderer Stelle!) wurde ein ebenes Diacramm 
zur Darstellung der Gleichgewiehtsisobaren und 
Gleichzrewiehtsisothermen von ternären Gemischen 
entwickelt und seine Anwendung am Flüssiekeits- 
Dampfgemisch Sauerstofi-Stickstoff-Argon gezeigt. 
Im nachfoleenden wird die Verwendung eines 
solehen Diarrammes zur Darstellung von Zustands- 
kurven ternärer Gemische am Beispiel der Entropie- 
vermehrungz im idealen ternären Gemisch erläutert. 

l;s seien a, b, ec drei Stoffe, die sämtlich ent- 
weder im dampfförmigen, flüssigen oder festen Zu- 
stand sind. e sei die Entropievermehrung je Mol 
des Gemisches. Sind ferner 24. 2b» Ze die Mol- 
bzw. Volumenkonzentrationen. so ist 

N 1.985 (2. In 2a 


„Inz»—+ 2.1Inz.). (1). 


und bei Einführung der Briggschen Logarithmen 


Ü 4,5655 (2, 12 24 + ale zu + 2.122.) (02). 
wobei 
Aa n } 1 (3) 
vilt. 


X\uf Grund von (3) können wir die e als Ordi 
naten in einem ebenen Diarramm mit zwei Ordi 
natenachsen (darstellen. die den Abstand 1 von- 
einander haben. und dessen Abszissen die drei 
Konzentrationen 2,. 2, und 2. sind. Zur Verein- 
fachunzr wählen wir als Ordinaten ‚-.-.. 
1.5655 

Wir setzen nun in (2) 2, konstant und schreiben 
daher 


N 1.5695 (A — zu le zu 


ze leg 2.) . (4). 


Für A# wählen wir der Reihe nach 2,=0,0.1,0,2 
usw, und erhalten damit aus (4) die in Abb. 1 dar- 
eestellte Kurvensehar. Diese ist nach unten dureh 
die gestrichelt gezeichnete Kurve VÜ1 begrenzt. 
die die Entropievermehrung e,. des binären Ge- 

1) V, Fischer: Ann. d. 
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misches der Stoffe b und e bestimmt. >ie ist dureh 


die Gleichung 


TyR 4,565 (zulezu + 2elg£2) -» » » Ö) 
eeeeben. Die weitere Begrenzung der Kurvenschar 


erfolgt durch das Ordinatenstück I A, dessen Länge 
nach (5) gleich —1g0,5 ist, entsprechend dem 
Höchstwert der Entropievermehrung beim binären 
(Gemisch. und durch die gestrichelt gezeichnete Ein 
hüllende der Kurvenschar. 

Ziehen wir die Gerade O0 A, dann wird BU=DE, 
so daß die Ordinatenstücke. die durch die Ein 
hüllende und die Gerade 0 A begrenzt sind. eben- 
falls die Entropievermehrung beim binären Gemisch 
ergeben. 

Wie aus Gl. (2) folgt, liegen die Schnittpunkte 
C, G und F je zweier Kurven der Schar mit dem 
Berührungspunkt E der Einhüllenden auf einer 
eemeinsamen Ordinate. Daraus ergibt sich eine 
einfache Konstruktion der Kurvenschar, sobald 
die ersten Kurven aus (2) berechnet sind. 

Aus dem Diagramm der Abb. 1 lesen wir z. B. mit 
2 0.2, 2» = 0.4 und 2.0.4 für den zugehörigen 
Punkt F die Entropievermehrung e = 4,5695 0,458 

209 ab. Denselben Wert e erhalten wir auch 
für 2,=0,4, »=0,2 und 2. =0,4. 


Für den Höchstwert der Entropievermehrung 
beim ternären Gemisch gilt 
Ve Np 
oe \ O8 h 
0 wand |—) =0.... 0. 
0 2a’: 0 Ah’za 


Damit folgt aus (1) und (3) 


A I) “C ’ 15) 


und 

u + re 
Das gleiche muß sich auch für den Höchstwert der 
Einhüllenden ergeben, wie nachfolgend gezeigt wird. 
Bezeiehnen wir die Abszissen der Einhüllenden mit 
2, und zy’. die Ordinaten der Geraden 0 4 mit «, 
und beachten wir, daß 1A4=1e2 ist. so lautet die 
Gleichung dieser Geraden 


Rp v. a lg . . . ö . . . . (9). 


Daher eilt für die Gleiehung der Einhüllenden 


1 


e 1.085 (2a In > t 2a In 2a’ 2» In u) (10). 
Daraus erhalten wir für den Höchstwert von e 
a2 N 
Öc 
y R — In .) + In By ” In zu _—— (0) (11 ). 
( “A pe 
Mithin wird 
1 MR 1 
x .) 
o # . > (12) 
“(1 _ 
und 
4) 1 
a # — PER oy\ 
<a 3 . <) 3) . . . . . . (13). 


Führen wir (13) in (10) ein. so ergibt sich 


e = 4.5655 le 3 
in Übereinstimmung mit (8). 


Während die Entropievermehrung je Mol des 
idealen Gemisches für sämtliche Stoffe gleich wird. 
ist dies für die Entropievermehrung je Massen- 
einheit nicht der Fall. Bezeichnen wir diese mit 
e, das Mol eines Gemisches mit m und die Mole- 


8 IE Ye 





rch 


(5) 
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em 
ren 
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:ulargewichte seiner drei Bestandteile mit n,„, 
77, und m., So wird 


M—Ma2Zat my u» Me 2e -» . (14), 


l 


ınd es gilt für verschiedene Gemische mit den 
2 


Molekulargewichten m, m’, m’’ usw. 


t=en=e'm=e"mM"—= .. (15). 

Das Gemisch bestehe z. B. aus Argon, Sauerstoff 
und Stickstoff, und wir wollen für die Zusammen- 
setzung 24. =0.2, ,=0,4 und z,=0,4 aus dem 
Diagramm Abb. 1 die Entropievermehrung be- 
stimmen. Es ist jetzt m„=40, m,=32 und 
n„=28. Damit folgt aus (14) m=32. Zu der 


u Dr er ae 


Oy4r— + t 
940 m 
| 
036} ! 
| Y 
032} 7 
028 t 





Unter Einführung des mechanischen Wärmeäqui 
valentes A folgt aus (17) für die Trennungsarbeit a 


je Mol 


Ja=Ama 1.9855 T (zaln2a- 2n Inz, | 
(18). 
2e In 2) 
Nach (1) ergibt sich daraus 
TE Ve 
und 
u N 05 £ See; 3 











7) 07 02 0 


[44301] 


Abh. 


obiren Zusammensetzung des Gemisches haben wir 

bereits aus dem Punkt F des Diagrammes € = 2,09 er- 
2.09 u 

— (), 693. 
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Le Ze 


mittelt. Aus (15) erhalten wir nun e = 


Unter der theoretischen Arbeit zur Trennung 
eines Gasgemisches in seine Bestandteile versteht 
man bekanntlich jene Arbeit, die notwendig ist, 
um die einzelnen Bestandteile von ihrem Teildruck 
im Gemisch isothermisch auf den Druck zu ver- 
diehten, unter dem das Gasgemisch steht. 

Bezeichnen wir die Teildrücke in einem ternären 
(sasgemisch, das unter dem Druck p steht, mit Pu. 
?n und pe. den spezifischen Rauminhalt des Ge- 
misches mit ® und seine Teilvolumina mit ”,. ®% 
und ©., sowie die Trennungsarbeit je Masseneinheit 
mit a, dann wird demnach 

a=pv(‘"in F +2 4. Ey = . (16). 

® Da v pn v Pe 
Bezeichnen wir ferner die Gaskonstante des Ge- 
misches mit A, seine absolute Temperatur mit 7 


v ) ! r 

und beachten, daß e = = — 2,, usw. ist, so können 
) ) 

wir (16) auch schreiben 


a=-- RT(z„Inzua+2zlnz»-+ 2eInz.) . (17. 





l. 


Wir können also auch die theoretische Arbeit 
zur Trennung eines ternären Gasgemisches aus 
dem Diagramm Abb. 1 bestimmen. Sei z. B. für 
ein A» 0 N Gemisch die Zusammensetzung 
wieder 2, =0.2, 2, =0.4, z2,=0.4 und die Tempe- 
ratur T == 288’ abs! dann wird nach (20) 


a == 427-288 0,0653 = 8030 mke/ke. 


das sind, falls in einer Stunde I kg des Gemisches 
zerlegt werden soll, 0.0297 PS. 


Zusammenfassung. Aus der Bedingung, dab die 
Summe der drei Konzentrationen eines ternären 
Gemisches zleich 1 ist, ergibt sich die Möglichkeit, 
die Zustandskurven ternärer Gemische in einem 
ebenen, orthogonalen Koordinatensystem darzu- 
stellen. Es wird als Beispiel die Kurvenschar 
ermittelt, durch die die Entropievermehrung beim 
idealen ternären Gemisch in dem neuen Diagramımı 
dargestellt ist. Ferner wird eine Gleichung ab- 
geleitet, die den Zusammenhang zwischen Entropie- 
vermehrung und theoretischer Trennungsarbeit 
ergibt, so daß das erhaltene Diagramm auch die 
letztere abzulesen gestattet. 
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BUCHBESPRECHUNGEN 


FELIX KLEIN, Vorlesungen über die 
h  ypergeometr ische Funktion, eehalten 
a. d. Univ, Göttingen im WS. 1893/94.  Ausgear- 
beitet von Ernst RITTER, herausgeg. u. m. Anm. 
versehen von Otto HAUPT. Prof. d. Math. a. d. 
Univ, Erlaneen. (Bd. 39 d. „Grundlehren der ma- 
themat. Wissenschaften“, herausgeg. v. R. Cou- 
rant.) IX + 344 8. m. 96 Fig. Berlin 1933, Verlag 
von Julius Springer, Preis geb. 23.60 RM. 


Die vorlierende Vorlesung hat Klein im Winter 
1803/94 vor einer zerineen Zahl von Hörern ge- 
halten. War doch auch in Göttingen in jenen Jahren 
lie Zahl der Mathematikstudierenden sehr zurück- 
veeaneen. Die Hörer waren überwiegend damals 
\usländer, unter ihnen zum ersten Male auch zwei 
Stiulentinnen der Mathematik. deren Dissertationen 
übrigens aus dieser Vorlesung erwachsen sind. Für 
den Referenten war es die erste Vorlesung, die er 
bei Klein zehört hat. Um es gleich vorauszu- 
nehmen: Die Absicht des Herausgebers. möglichst 
den persönlichen Kleinschen Stil zu wahren, ist 
sehr gut gelungen. Das soll um so mehr hervor- 
eehoben werden. als der Herausgeber bei Klein 
selbst nie eine Vorlesung gehört hat. Referent sah 
beim Lesen unwillkürlich Klein vor sich auf dem 
Katheder stehend, mit dem Folioblatt in der Hand. 
auf dem der Gedankengang jeder Vorlesung sorgfältig 
‚lisponiert war. Es steigt aber auch die Erinnerung an 
den genialen Assistenten Kleins auf, der die Vor- 
lesung «damals ausgearbeitet hat: Ernst Ritter. 
der noch nieht zwei Jahre später (September 1895) 
im Hafen von New York einem schweren Typhus 
zum Opfer fiel. der ihn auf der Fahrt nach Amerika. 
wohin er als Professor berufen war. befallen hatte, 
Bei einer neuen Auflage der gedruckten Ausgabe 
möchte im Vorwort doch auch Ernst Ritters 
etlacht werden unter Hinweis auf den Nachruf. den 
Klein seinem Schüler „ewidmet hat. (Jahres- 
bericht der Deutschen Mathematiker-Vereinigung, 
Band IV, Seite 52-54: versehentlich ist dort bei 
der Jahresangabe ein Druckfehler stehengeblieben: 
es mub 18095. nieht 1885 heißen.) 

Absicht der Kleinschen Vorlesung war es, eine 
erobe Funktionsgattung, die in der damaligen Zeit 
nicht hinreichend allgemein studiert war, im Gegen- 
satz z. B. zu den elliptischen Funktionen und deren 
Veralleemeinerungen, zur Geltung zu bringen unter 
besonderer Anlehnung an die eroße Arbeit Rie- 
manns über diese Funktionen, Mit der Analyse 
lerRiemannsechen Arbeit beschäftigt sich der erste 
Feil der Vorlesungen. wobei die hypergzeometrische 
Funktion. die Riemannsche P-Funktion. als 
l,ösung, einer Differentialgleichung zweiter Ordnung 
erscheint. Der Ausdruck hypergeometrisch kommt 
zuerst bei Wallis vor für eine Reihe. deren Bil- 
dungsgesetz allgemeiner ist als das der geome- 
trischen Reihe. In dem wertvollen, bis auf die 
neueste Zeit erzänzten Literaturverzeiehnis hätte 
Wallis vielleieht auch mit einem eenauen Zitat 
bedacht werden können !). | 


!) Joh. Wallis: Arithmetiea Infinitornm Opera math. 
Vol. 1 (1690) pag. 466, 4. Absatz. In der vom 19. Juli 1655 
datierten, für die Gesehiehte seiner Forschung sehr inter- 
essanten ausführlichen Widmung der Arithmetiea an den 
Lkandpfarrer W. Ougthred (matheseos peritissimum) 
sagt er, daß er für Reihen, die noeh mehr zusammen 
gesetzt sind als geometrische, geglaubt hat, einen nenen 
Namen einführen zu müssen, um so mehr, da nieht zu 
hoffen sei, aus solehen Reihen dureh „Interpolation“ 
mittlere Glieder als „erklärbare Zahlen“ i(numeri ex 
plieabiles) zu gewinnen; erhält man doch sehon hei der 
gewöhnlichen geometrischen Reihe „unmögliche Zahlen“ 
wie V 2. Den Ausdruck „hypergeometrisch" nennt er 
dort aber nieht. Op. math. Vol. I pag. 359. Dagegen er- 
wähnt er ihn ausdrücklieh in einem Brief an Leibniz 
vom 16. Januar 1690: Series quam dieo hypergeometrieam 
12624 120 ete., Operum mathematicorum Vol. Tertinm 


Die Reihe 


2 ab ata+VDb(b+N1) 
F(a,b:c: x) 1 + cH x? 
l-c 1-2-c(c+]) 


hat Euler zuerst als eine Funktion von 7 be 
trachtet. 

Der im Literaturverzeichnis S. 307 geenanntı 
Band der Opera Omnia, Serie I, Band 16, ist in 
wischen auch erschienen. Klein fabt den Begrift 
etwas weiter und versteht unter einer hyper 
eeometrischen Funktion 


Creed r)’-F(a.b:c:r). 


wo a und zunächst reelle Zahlen sind. In der 
eeschiehtliehen Entwieklung wird Pfaff genannt, 
der bedeutende Helmstedter Mathematiker, bei dem 
Gauß mehrere Monate gewohnt und bei dem er 
auch promoviert hat. Bei der Nennung von 
Pfaff kann die wortgetreue Fassung der jetzigen 
Druckausgabe für jüngere historisch nicht genügen!l 
orientierte Leser vielleicht irreführen: statt zu 
saeen „Ende des vorigen Jahrhunderts“, was für 
1893 paßte, möchte man jetzt lieber sagen: „Ende 
oder riehtiger aus der Wende des 18. Jahrhun 
derts“; ist doch Pfaff, der 1810 nach der Auf- 
hebung der Universität Helmstedt als Ordinarius 
nach Halle gekommen ist. dort 1825 gestorben, 
Wesen der von Klein angeregten Frage, die 
näheren Beziehungen der Arbeit Pfaffs zu deı 
späteren großen Arbeit von Gaub zu untersuchen. 
könnte auf Schlesinzer verwiesen werden: 
Über Gauß’ Arbeiten zur Funktionen-Theorie. Mate 
rialien für eine wissenschaftliche Biographie von 
Gauß 1912. S. 35. Freilich fehlt auch heute noch 
die in Aussicht gestellte genauere Darstellung des 
wissenschaftlichen Verhältnisses der beiden Männer. 

Der zweite Teil behandelt die durch den Quo- 
tienten zweier Lösungen der Differentialgleichung 
zweiter Ordnung  zelieferte konforme Abbildung. 
Dieser Quotient genügt einer Gleichung dritter 
Ordnung. Es ergeben sich Beziehungen zu der Rieca- 
tischen Differentialgleichung zu den Bessel 
schen Funktionen und anderen. Eingeschoben wirt 
ein Ausblick auf Lies Auffassung. Zum Studium 
des Quotienten dient eine Darstellung auf der kom 
plexen Kugel. und so ergibt sieh als Hilfsmittel die 
sphärische Trigonometrie; sie wird durch die Namen 
auß. Möbius, Study charakterisiert. Dieser 
Absehnitt wendet sich dann zu den Ergebnissen von 
F. Scehilling. dureh die die Formeln der 
sphärischen Trironometrie für komplexe Werte der 
Veränderlichen geometrisch erläutert werden. In 
der Vorlesung selbst konnte Klein nur aul 
die bevorstehenden Untersuchungsergebnisse Schil- 
lines hinweisen, die damals noch nicht abge 
schlossen waren. 

Wenn mit Rückicht auf den zur Verfügung 
stehenden Raum hier nur stichwortartig Gruppen 
theorie. Algebra mit besonderer Rücksicht aufs 
Ikosaeder, Nichteuklidische Geometrie. Gedanken- 
kreis von Pieeard-Vessiot aufgezählt werden. 
so mub das eenügen, einen Einbliek in den reichen 
(edankeninhalt der Klein schen Vorlesung zu geben, 
(die wie jede Kleinsche Vorlesung immer von dem 
Bestreben «eleitet ist. die verschiedenen Gebiete 
der Mathematik miteinander zu verbinden. Mit 
Recht sagt der Herausgeber im Vorwort. daß diese 
Vorlesung eine treffliche Ergänzung und Weiter- 
führung dessen bringt. was der Studierende mitt- 


(1699) pag. 694. An die gleiche Reihe der Fakultäten 
knüpft Enler in der nach seinem Tode veröffentliehten 
Petersburger Akademieabhandlung an: De termino gene 
rali serierum hypergeometriearum. Op. omnia Ser. 1, 16, 
S, 14, 


- 
Pa 


u 


EW 


:ch. 


nt« 
in 
rifi 


ler 
nt, 
em 


'On 
"en 
nl 
ZU 
für 
le 
11 
uf- 
jus 
el, 
lie 
leı 
BR. 
3: 
te 
on 
ch 
les 
ET. 
LO- 
ne 


Bin. 
A - 


rl 
111 
In 

lie 
el 
IST 
ON 
er 


In 
ul 


Fü 
n 


n- 
N. 
N 
N, 
m 
te 
jt 

se 


T- 


el 
el 
ie 

Ib, 


2" Kr 


“ 


— 
5 “ 


sand 14, Heft 4 
August 1934 


Nachrichten 2355 





rer Semester an Geometrie und Funktionentheorie 
ennen gelernt hat. Sehr wertvoll sind auch die 
‚3 Seiten Anmerkungen, durch die der Herausgeber 
. a. die seit 1893 erzielten Fortschritte kenn- 
ejehnet, 

Leipzig, z. Zt. Frankfurt a. M. 


W. Lorey. 


H. KAYSER, Prof. der Physik an der Universität 
Bonn, und H. KONEN, !’rof. der Physik an der Uni- 
versität Bonn, Handbuch der Speectro- 
scopie. T. Bd. 3. u. letzte Liefg. Leipzig 1934, 
Hirzel Verlag. Umfang des 7. Bandes (enthält 
3 Lieferungen) IV + 1473 Seiten. Preis brosch. 
32,20 M., geb. 141,— M. 


Mit der vorliegenden dritten Lieferung ist der 
Band des Handbuchs der Speetroscopie von 
nayser undKonen abgeschlosssen. Der Gesamt- 
und umfaßt 57 Elemente in alphabetischer Reihen- 
'olee von Argon bis Niobium. Bei der ungeheuren 
"ülle der spektroskopischen Publikationen der letz- 
ten Jahre haben die Verfasser den Anschluß des 
behandelten Stoffes an die Jetztzeit durch ein An- 
näherungsverfahren zu erreichen versucht. Die erste 
Lieferung des 7. Bandes vom Jahre 1923 ist durch 
lie inzwischen erschienene erste Lieferung des 
Sn, Bandes auf den Stand des Jahres 1932 gebracht. 
Die zweite Lieferung dieses 8. Bandes ist noch für 
1934 in Aussicht gestellt, so dab dann in den beiden 
Bänden 7 und 8 die genannten 57 Elemente so nahe 
wie möglich an die Gegenwart herangeführt sind. 
\uech in dem neuen Band ist die bekannte und be- 
währte Form beibehalten worden. Die spektroskopi- 
schen Daten sind so vollständige angeführt, daß ge- 
ade auf dieser restlosen Umfassung eines umfang- 
reichen Teilgebiets der Physik die hohe Bedeutung 
les Handbucehs für die Forschung beruht. 


J. 


W. BAUMANN und R. MECKE, Das ultra- 
'ote Sonnenspektrum von A 10000 bis 
, 7600 A. E. Leipzig 1934, Johann Ambrosius Barth. 
56 8. m. 6 Taf. Preis 12 M. 


In sechs ausgezeichnet wiederzegebenen Tafeln 
ist das ultrarote Sonnenspektrum zwischen 7574 und 
809 A.E. mit einer Dispersion von 0.44 A. E./mm 
(largestellt. Diese Aufnahmen bilden eine erwünschte 
Krweiterung des bekannten Rowlandschen 
Sonnenatlasses, der nur bis 7331 A.E. reicht. In 
len beigefügten Tabellen. die sich z. T. auf eigene, 
7. T. auf amerikanische Messungen stützen, sind 
3297 Linien zwischen 10000 und 7600 A.E. ent- 
halten, von denen eine große Zahl durch Banden- 
analyse wie auch durch direkte Absorptionsversuche 
dem Wasserdampf und dem Sauerstoff mit Sicher- 
heit zugeschrieben werden können. Diese genaue 
l’urehmusterunge dürfte besonders dem  Astro- 
physiker willkommen sein, dessen Interesse durch 
len kürzlich erfolgten Nachweis des Methans und 
\mmoniaks in einzelnen Planetenatmosphären ge- 
rade auf dieses Gebiet der Absorptionsspektroskopie 
velenkt wird. 

Berlin. 
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Ferner sind bei der 
bücher 
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Schriftleitung 
(ausführliche 


folgende 
eingegangen Besprechung 


vorbehalten): 


Dr.-Ing. Dr. phil. R. POZDENA, Meter und 
Kilogramm. Entstehung und Sicherung (des 
internationalen metrischen Maßsystems.  Mathem. 
Phys. Bibl.. Reihe 1, Bd. 76. 45 8. m. 15 Fig. 
leipzig und Berlin 1934, Verlag B. G. Teubner. 
Preis kart. 1.20 M. 


Dr. Dr.-Ing. GEORG SCHEFFERS, o. 
Teehn. Hochsch. Berlin. Wie findet und 
zeichnet man Gradnetze von Land 
und Sternkarten? Mathem.-Phys.  Bibl.. 
Reihe 1. Bd. 85/86. 98 S. m. 12 Taf. u. 27 Abb. 


Prof. a. ıl 


leipzig und Berlin 1934, Verlag B. G. Teubner. 
Preis kart. 2,40 M. 
Dr.-Ing. FERDINAND WINTER \VDlIl. Das 


Spiel der 30 bunten Würfel. 128 SS. m. 
31 Fig. u. 27 Tab. Leipzig und Berlin 1934, Ver 
lax B. G. Teubner. Preis kart. 3,60 M. " 

R. BECKER, o. Professor a. d. Techn. Hochseh. 
Berlin, Theorie der Elektrizität. Band |. 
Kinführung in die Maxwellsche Theorie der Elektri- 
zität. mit einem einleitenden Abschnitte über (as 
Rechnen mit Vektorgrößen in der Physik. 10. Aufl. 
262 S. m. 59 Fig. Leipzig und Berlin 1933. Verlag 
B. G. Teubner. Preis geb. 14.50 M. 


C, JUEL, Professor emeritus an der Technischen 
Hochschule Kopenhagen, Vorlesungen über 
’ro ] ektive Geometrie, mit besonderer Be 
rücksiehtieung der v. Staudtschen Imaginärtheorie 
Band XLI. Die Grundlehren der Mathematischen 
Wissenschaften in Einzeldarstellungen, 282 S. m. 


87 Fig. Berlin 1934, Verlag Julius Springer. Preis 
eeb. 2250 M. 
MARIO MEDICI, Prof. a. d. Technischen Hoch 


schule Padua, Le Machine Idrauliche. 
Biblioteca DellIngegenere Bd. 3 (Herausg. Prof. 6. 
Albenea). 520 S. m. 325 Fire. Torino 1934 \1l. 
Unione Tipogratico-Editriee Torinese. 


F. ODQVIST, Dr. phil. 
'Jastieitetsteori meld 
ear. 80 8. Stockholm 1934. 
erafiska Anstalts Förlag. 


Dozent. Stockholm. 
Tillämpnin 
(eneralstabens Lito 


Dr.-Ine. © K. FRÖHLICH, Druckver 

teilune im Baugrunde mit besonderer Be 
rücksichtigung der plastischen Erscheinungen. 
183 S. m. 68 Textabb. Wien 1934, Verlag ‚Julius 
Springer. Preis kart. 15 M. 


Abhandlungen der Internationalen Vereinigung 
für Brückenbau und Hochbau, Band 2. 1933-1991 
Herauszereben vom Generalsekretariat in Zürich. 
495 8. 1934. Kommissionsverlag A.G. Gebr. Lee 
mann & Co.. Zürich. Preis kart. 34.80 Schw. Franken. 


NACHRICHTEN 


IV. Internationaler Kongreß für angewandte 
Mechanik. 


Vom 3. bis 9. Juli fand in Cambridge der vierte 
Internationale Kongreß für angewandte Mechanik 
statt, der sieh, wie die früheren Kongresse, eines 
sehr zahlreichen Besuches der Fachgenossen aus 
aller Herren Länder erfreute. Nachdem der inter- 
nationale Kongrebauschuß bereits am 
des 3. Juli unter dem Vorsitz des Präsidenten, Prof. 


Vormittag 


C. E.1Ingelis, zu einer Sitzung zusammengetreten 
war, erfolgte die feierliche Eröffnung des Kon 
eresses durch den Vizekanzler der Universität Cam- 
bridge, Herrn J. F. Cameron (M.A, Master of 
Conville u. Caius College), am Nachmittage im 
Senatshause der Universität. Die sachlichen Ver 
handlunzen des Kongresses begannen am folgenden 
Tage. Mittwoch, 4. Juli: sie wurden in Vor- und 
Nachmittagssitzungen (mit Ausnahme des Sonntags) 
bis Montage. 9. Juli, durchgeführt. Für die all 
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gemeinen Vorträge stand je eine volle, für die 
Sektionsvorträge einschließlich Diskussion eine 
halbe Stunde zur Verfügung. Wie schon in Stock- 
holm, so war auch diesmal in den Nachmittags- 
stunden nach Schluß der Sitzungen Gelegenheit zur 
Fortführung der Diskussion eeboten, soweit ein be- 
sonders lebhaftes Interesse an einzelnen Vorträgen 
dies erwünscht erscheinen ließ. Auf diese Weise 
wurde die Einhaltune des Programms ermöglicht. 
ohne dab die Diskussionsmöglichkeiten über Gebühr 
eingeschränkt wurden. 

Die allgemeinen Vorträge, «deren Themen im 
letzten Heft mitgeteilt wurden, erregten naturgemäß 
das Hauptinteresse, da sie einen Überbliek über die 
Fortschritte auf den gegenwärtig wichtigsten For- 
schungsgebieten gaben. Zur Mechanik der 
Flüssigkeiten eab Kärmän einen Überblick über 
die Entwicklung der Vorstellungen und über den 
gegenwärtigen Stand der Theorie «der turbulenten 
Strömung, zu der er selbst so wesentliches beizetraren 
hat. Wagner sprach über das Gleiten von Körpern 
auf der Wasserobertläche, wobei er sich gleichfalls 
auf eigene, frühere und neuere Arbeiten stützen 
konnte. Zur Schwingungstheorie sprach Den 
Hartog in einem außerordentlich anregenden, 
dureh Modellvorführung ergänzten Vortrag, der an 
sehr interessanten Beispielen einen Überbliek über 
die Fülle der Probleme auf diesem praktisch so 
wichtigen Gebiete der Ingenieurmechanik bot. 
Die Klastizitäts- und Festiekeitslehre kam in 2 der 
allgemeinen Vorträge zur Geltung. Neben einem 
anregenden Vortrag von Caquot über grundsätz- 
lieh wiehtige Fragen der Elastizitätstheorie fesselte 
der Vortrag von Taylor (Cambridge). der einen 
neuen und  vielversprechenden Ansatz zur Er- 
klärung der Festirkeitseigenschaften der Kristalle 
vah, Der praktischen Ingenieurtätiekeit am 
nächsten stand wohl der Vortrag von Schmidt 
Danzig). der über theoretische und experimentelle 
Untersuchungen zur Frage des Wärmeüberganges 
von festen Körpern an flüssige und gasförmige 


Stoffe berichtete. Über Rechenmaschinen uni 
mathematische Instrumente sprach Bush (M.ILT. 
Cambridge U.S.A.). Vom deutschen Standpunkt 


aus muß auf diesen Vortrag@ besonders hingewiesen 
werden, weil es sich hier um ein Gebiet handelt, 
las mehr Aufmerksamkeit verdient, als ihm in 
Deutschland gerenwärtie zuteil wird. Es sei er- 
wähnt, daß u. a. eine Maschine zur Lösung von 
linearen Gleiehungen mit 11 Unbekannten vor 
eeführt wurde, 

Über die Einzelvorträge in den Sektionen zu 
berichten, ist schon deshalb kaum möglich 
weil der Referent von den gleichzeitige statt 
findenden Vorträgen nur jeweils einen hören 
konnte, leider wird es nieht möglich sein. wie in 
früheren Jahren einen Kongreßbericht zu drucken. 
in dem alle Vorträge ausführlich wiedergegeben 
wurden. Der Verhandlungsberieht wird nur die 
alleemeinen Vorträge in extenso enthalten, von den 
Sektionsvorträgen nur kurze Auszüge. 

Der internationale Kongreb-Ausschuß hielt seine 
Schlußsitzung am Sonnabend, dem 7. Juli. An Stelle 
der Herren vv Karmän und v. Mises, die als 
Vertreter der UNS.A. und der Türkei dem Aus 
schuß weiterhin angehören werden, wurden die 
Herren Grammel und Trefftz als deutsche 
Vertreter zugewählt. ferner Herr Federhofer 
für Österreich. Die Schlußsitzung (des Kongresses 
fand am Montag, dem 9. Juli, nachmittags, statt. 

Es braueht wohl nicht besonders gesagt zu 
werden. daß die Teilnehmer des Kongresses (ie 
freundliehsten Erinnerungen an die Cambridger 
Taeune mitzenommen haben. Neben dem wissen 


schaftlichen Teil. um den sieh die Herren des Orga 
nisationsausschusses besonders bemüht haben, ver 
dient der gesellschaftliche Teil, in dem die traditionell: 
englische Gastlichkeit zur Geltung kam, eine dank 
bare Erwähnune. Den Herren und Damen des 
Orranisationsausschusses, dem Mavor (der >taili 
Cambridge und den Masters and Fellows von 
Christ's College, St. John’s College, Trinity College, 
Kine’s Colleee und Pembroke College ist deı 
dauernde Dank der Teilnehmer sicher. Nicht zu 
vereessen: Herr P. ©. Fitzgerald M.A.. der einen 
iiberaus reizvollen Vortra® über die Geschichte uni! 
die Gebäude von Cambridge hielt und sich uner 
müdlich zeigte, den ausländischen Gästen Ver 
eaneenheit und Gegenwart von Cambridge nahe 
zubringen. Trefftz, 


Gesellschaft für angewandte Mathematik und 
Mechanik. 

Für die Tagung in Pyrmont sind die folgenden 

Vorträge angemeldet worden: 

Freitag, den 14. September, vormittags 9 Uhr. 

1.6, Sehulz. Berlin. Umkehrung von Integral 
transformationen. 

>, H. Sehaefer, Jena. Angenäherte Berechnung 
des kleinsten EKizenwertes zusammengesetzteı 
Systeme. 

3. L. Collatz. Berlin. Verallzemeinerung des 
DifferenzenverfahrensfürDifferentialgleichungen. 

I. W. Cauer, Göttingen. Topologische Dualitäts 
sätze und Reziprozitätstheoreme der Schaltungs- 
theorie. 

5. R.Höger,. Prag. Schwingungen elastischer Seile. 

6. W. Quade, Karlsruhe. Schwingungsvorgänge 
in Systemen mit zwei Freiheitsgraden. 


Nachmittags 16 Uhr. 


1. H. Neuber, München. Theorie der räumlichen 
Kerbwirkung. 

2, H. Fromm. Berlin. Zerreißversuche an Flach- 
stäben mit ebener Formänderung. 

3. H. Ebner, Berlin. Spannungszustand dureh 
Drillune in dünnwandigen Kastenträgern bei 
verhinderter Endverwölbung. 

4.K. Ludwig. Hannover. Kinwirkung von 
Bodenschwinzungen auf Pfeiler. 


Sonnabend, den 15. September, vormittags Y Uhr. 
Miteliederversammlung. 
Anschließend: 

1. W., Müller. Prag. Über einige nieht-stationäre 
Bewerungzen von reibenden Flüssigkeiten. 

> W. Tollmien. Göttingen. Ein Instabilitäts- 
kriterium für laminare Geschwindigkeitsver- 
teiluneen. 

3. H. Gebelein, Stuttgart. Zustandekommen 
der turbulenten Austauschgröße. 

I. C. Weber, Dresden. Wellenerregung in einem 
Kanale dureh eine bewegte Wand. 

5.0, Weber, Dresden. Theoretische Ansätze 
für die Schulersehen Umschlagwellen. 

6.6. Vogelpohl, Berlin. Über die Anwendung 
des Impulssatzes der Strömungslehre auf den 
Strahlantrieb (Raketenwirkung). 


Rechenmaschinen. 

Wir weisen unsere Leser darauf hin. daß in deı 
Internationalen Büroausstellung. die vom 7. bis 
16. September in Berlin stattfindet, Gelegenheit 
vereben ist, die neuesten deutschen und aus 
ländischen Rechenmascehinen kennenzulernen. 
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